Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,647,706 T→A 8.7% Y119N (TAT→AAT)  cysB → HTH‑type transcriptional regulator CysB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,647,7060TA 8.7% 135.8 / 4.7 46Y119N (TAT→AAT) cysBHTH‑type transcriptional regulator CysB
Reads supporting (aligned to +/- strand):  ref base T (24/17);  new base A (2/2);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

CTGCGCAAGGTCGAGAGCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGC  >  ACB143/2647571‑2647840
                                                                                                                                       |                                                                                                                                      
cTGCGCAAGGTCGAGAGCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGt                                                                                                                                        <  1:90908/136‑1 (MQ=255)
    gcAAGGTCGAGAGCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATcc                                                                                                                                    <  1:49726/136‑1 (MQ=255)
     cAAGGTCGAGAGCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCCCGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGGTATCAGAAGCTTCATCAAGCAGTATCCg                                                                                                                                   >  2:306639/1‑136 (MQ=255)
                gCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGGAAAGCGTTGCCGCCCGTTATTAGCAGCATCATCAAGGATGTTCCGGAGGTGTCCCt                                                                                                                        >  1:338462/1‑136 (MQ=255)
                gCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCCGTATCCGGAGGTGTCCCt                                                                                                                        >  8:131099/1‑136 (MQ=255)
                gCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCt                                                                                                                        >  6:298689/1‑136 (MQ=255)
                gCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCt                                                                                                                        >  7:25601/1‑136 (MQ=255)
                 cATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTg                                                                                                                       <  6:308646/136‑1 (MQ=255)
                  aTCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGc                                                                                                                      <  2:388395/136‑1 (MQ=255)
                      aGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCTCCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACAt                                                                                                                  >  5:181026/1‑136 (MQ=255)
                        cAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGc                                                                                                                <  4:130147/136‑1 (MQ=255)
                              gCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAgg                                                                                                          >  5:114764/1‑136 (MQ=255)
                              gCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAgg                                                                                                          >  5:22516/1‑136 (MQ=255)
                              gCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAgg                                                                                                          >  7:255244/1‑136 (MQ=255)
                                            cAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATcc                                                                                                                                    <  4:7557/96‑1 (MQ=255)
                                            cAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATcc                                                                                                                                    >  3:7557/1‑96 (MQ=255)
                                            cAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATg                                                                                               >  7:112554/1‑133 (MQ=255)
                                            cAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGc                                                                                              <  5:45739/134‑1 (MQ=255)
                                            cAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGAATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATg                                                                                               >  3:233110/1‑133 (MQ=255)
                                            cAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGAATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATg                                                                                               >  3:238408/1‑133 (MQ=255)
                                                   aagaagGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCg                                                                                     <  1:143338/136‑1 (MQ=255)
                                                       agGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGAt                                                                                 >  2:133760/1‑136 (MQ=255)
                                                       agGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCAATGCAGATTGCCGAGAt                                                                                 >  3:267849/1‑136 (MQ=255)
                                                           cACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGAATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCt                                                                             <  4:233110/136‑1 (MQ=255)
                                                           cACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGAATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCt                                                                             <  4:238408/136‑1 (MQ=255)
                                                                   ccATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGacgg                                                                     >  6:40768/1‑136 (MQ=255)
                                                                     aTCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGacggta                                                                   <  6:377791/136‑1 (MQ=255)
                                                                        gCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGacggtacgg                                                                >  6:54839/1‑136 (MQ=255)
                                                                               cccacacccaGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGAttt                                                         >  4:391650/1‑136 (MQ=255)
                                                                               cccacacccaGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGAtt                                                          >  5:328871/1‑135 (MQ=255)
                                                                                  acacccaGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTtcgc                                                      <  1:361593/136‑1 (MQ=255)
                                                                                           cgcgcTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCAc                                             <  5:30990/136‑1 (MQ=255)
                                                                                           cgcgcTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCAc                                             <  1:54895/136‑1 (MQ=255)
                                                                                            gcgcTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATAGCCAcc                                            >  3:206489/1‑136 (MQ=255)
                                                                                             cgcTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCCCCg                                           >  2:235218/1‑136 (MQ=255)
                                                                                                tATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGAGTTAGCCCTCGCCCCCGAgg                                        >  6:230458/1‑136 (MQ=255)
                                                                                                     gTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAgg                                                                                                                               <  7:8211/44‑1 (MQ=255)
                                                                                                     gTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAgg                                                                                                                               >  8:8211/1‑44 (MQ=255)
                                                                                                     gTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTg                                   >  6:313616/1‑136 (MQ=255)
                                                                                                     gTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTg                                   <  3:18944/136‑1 (MQ=255)
                                                                                                       tGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCCCCGAGGCGCTg                                   >  2:37177/1‑134 (MQ=255)
                                                                                                             cgGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGtt                           <  5:40768/136‑1 (MQ=255)
                                                                                                                  aTCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGGTGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGc                       <  4:193292/135‑1 (MQ=255)
                                                                                                                  aTCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGc                       <  1:362665/135‑1 (MQ=255)
                                                                                                                             catcaAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCat            >  1:322976/1‑135 (MQ=255)
                                                                                                                                  aGCAGTATCCGGAGGTGTCCCTGCACATTCATCAGGGTTAGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTAAGGTAGATTTCTCCATCGCCACCGAGGCGATGGCGCTGTTCGGCAACCTGCTAATGATGcc      >  8:275794/1‑136 (MQ=255)
                                                                                                                                    cAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGt    <  2:148027/136‑1 (MQ=255)
                                                                                                                                       tATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGc  >  7:184583/1‑135 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CTGCGCAAGGTCGAGAGCATCAAGCAGATTGCCCAGGAGTTCTCCAACGAGAAGAAGGGCACGCTGTCCATCGCCACCACCCACACCCAGGCGCGCTATGCGTTGCCGCCGGTTATCAGCAGCTTCATCAAGCAGTATCCGGAGGTGTCCCTGCACATGCATCAGGGTTCGCCCATGCAGATTGCCGAGATGGCTGCGGACGGTACGGTCGATTTCGCCATCGCCACCGAGGCGCTGGAGCTGTTCGGCGACCTGATCATGATGCCGTGC  >  ACB143/2647571‑2647840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑