Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,430,947 G→T 8.2% Q157K (CAA→AAA)  tsaC ← threonylcarbamoyl‑AMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,430,9470GT 8.2% 445.4 / 18.9 134Q157K (CAA→AAA) tsaCthreonylcarbamoyl‑AMP synthase
Reads supporting (aligned to +/- strand):  ref base G (59/64);  new base T (7/4);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 136 total aligned reads displayed.

CTGTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCC  >  NC_000913/3430798‑3431091
                                                                                                                                                     |                                                                                                                                                
cTGTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCGAATtg                                                                                                                                                  >  2:669285/1‑150 (MQ=255)
  gtgGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATtgtg                                                                                                                                                >  2:1094565/1‑150 (MQ=255)
     ggCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTgcgc                                                                                                                                             <  2:1782230/150‑1 (MQ=255)
      gCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGcgcg                                                                                                                                            <  1:522858/150‑1 (MQ=255)
       cTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGa                                                                                                                                           <  1:1962181/150‑1 (MQ=255)
        tATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggc                                                                                                                                          <  1:2329572/150‑3 (MQ=255)
        tATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGCTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggc                                                                                                                                          <  1:959898/150‑3 (MQ=255)
         aTCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggcg                                                                                                                                         <  2:8010/150‑4 (MQ=255)
          gcGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCAGCAGGCACAGCCGGGAACGCCGCGCCAAATTGTGCGCGAACt                                                                                                                                        >  2:1635915/2‑150 (MQ=255)
            ggATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACttg                                                                                                                                      >  2:303027/1‑149 (MQ=255)
            ggATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTc                                                                                                                                      >  1:806931/1‑150 (MQ=255)
             gATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCACAATTTGGcgcggcgttcc                                                                                                                                     <  1:1688978/150‑8 (MQ=255)
                  ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                >  1:2203275/1‑150 (MQ=255)
                  ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                >  1:553936/1‑150 (MQ=255)
                  ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                >  1:901805/1‑150 (MQ=255)
                  ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                >  1:955979/1‑150 (MQ=255)
                             aTAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAg                                                                                                                     >  1:2145877/1‑150 (MQ=255)
                             aTAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCATCTACTGTTCGACAAg                                                                                                                     <  2:2310970/150‑1 (MQ=255)
                              tAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAgg                                                                                                                    <  1:255604/150‑1 (MQ=255)
                               aGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGt                                                                                                                   <  2:456223/150‑1 (MQ=255)
                                ggTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTg                                                                                                                  <  1:2258091/150‑1 (MQ=255)
                                  tttCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGc                                                                                                                <  2:1454769/150‑1 (MQ=255)
                                     ccATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAAt                                                                                                             >  2:1247664/1‑150 (MQ=255)
                                     ccATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAAt                                                                                                             >  2:809706/1‑150 (MQ=255)
                                       aTTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATcc                                                                                                           >  1:1509273/1‑150 (MQ=255)
                                           tGTTACCCCTGTCGAAAGAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTc                                                                                                       >  2:2156029/1‑150 (MQ=255)
                                              tACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCGGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGAAAAGGTGCAAATCAACTCAAg                                                                                                    >  2:425583/1‑150 (MQ=255)
                                                ccccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGtt                                                                                                  >  1:332714/1‑150 (MQ=255)
                                                 cccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTg                                                                                                 <  2:1539534/150‑1 (MQ=255)
                                                  ccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTgg                                                                                                <  2:1145628/150‑1 (MQ=255)
                                                  ccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTgg                                                                                                >  1:329434/1‑150 (MQ=255)
                                                   cTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGc                                                                                               >  1:1638762/1‑150 (MQ=255)
                                                   cTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGc                                                                                               <  1:523473/150‑1 (MQ=255)
                                                     gTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGCTGGCAATCCACTCAAGTTGGCAc                                                                                             >  1:2399027/1‑150 (MQ=255)
                                                      tCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAAGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGAc                                                                                                                        <  2:1729218/122‑1 (MQ=255)
                                                      tCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAAGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGAc                                                                                                                        >  1:1729218/1‑122 (MQ=255)
                                                       cGAAACAGTTCCCCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTGTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTg                                                                                           >  2:2460484/1‑150 (MQ=255)
                                                       cGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTg                                                                                           >  1:1556383/1‑150 (MQ=255)
                                                       cGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTg                                                                                           >  2:164099/1‑150 (MQ=255)
                                                         aaaCAGTTCCCCCGTCAGGGCATCGCGGCTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACCACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCCCTCAAGTTGGCACTGGt                                                                                         >  2:2152383/1‑150 (MQ=255)
                                                               ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   >  1:2292878/1‑150 (MQ=255)
                                                               ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   >  1:1263064/1‑150 (MQ=255)
                                                               ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   >  1:283827/1‑150 (MQ=255)
                                                               ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   <  2:1962109/150‑1 (MQ=255)
                                                                  aCCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAg                                                                                <  2:1306487/150‑1 (MQ=255)
                                                                   cccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGc                                                                               <  1:1153843/150‑1 (MQ=255)
                                                                     cGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCTAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCgg                                                                             <  1:669285/150‑1 (MQ=255)
                                                                     cGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCgg                                                                             <  1:414735/150‑1 (MQ=255)
                                                                      gTCAGGGCATCGCGGATTTCTGAAGGGTTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGt                                                                            <  1:392990/150‑1 (MQ=255)
                                                                          gggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAc                                                                        >  2:1837370/1‑150 (MQ=255)
                                                                          gcgCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCCGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACAAGCGGTTTAc                                                                        >  2:2067153/3‑150 (MQ=255)
                                                                           ggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAcc                                                                       <  1:312378/150‑1 (MQ=255)
                                                                            gCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCa                                                                      >  2:2309125/1‑150 (MQ=255)
                                                                            gCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCa                                                                      <  1:617727/150‑1 (MQ=255)
                                                                                    gATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTg                                                              >  2:2089308/1‑150 (MQ=255)
                                                                                     aTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTgg                                                             <  1:1094565/150‑1 (MQ=255)
                                                                                      tttCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGc                                                            <  1:1467028/150‑1 (MQ=255)
                                                                                      tttCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGc                                                            <  2:1918074/150‑1 (MQ=255)
                                                                                       ttCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGca                                                           <  1:1295732/150‑1 (MQ=255)
                                                                                        tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcac                                                          >  2:2373580/1‑150 (MQ=255)
                                                                                        tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcac                                                          >  2:463841/1‑150 (MQ=255)
                                                                                        tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcac                                                          <  1:337856/150‑1 (MQ=255)
                                                                                         cTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGTGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                         >  2:2486184/1‑150 (MQ=255)
                                                                                           gAAGGATTTAAACGCCCCCCGGTATCACCAGGCACAACCGGGACCGCCGGCACAAATTGTGCCCGAAATACGTCAACTGTACGACAAGGTGGCAAACCACTCAAGTGGGCACTGGTAGAAAGATGCGGTTTACCATAAGCCTCCCACaaa                                                       >  2:343247/1‑150 (MQ=255)
                                                                                             aGGATATAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGc                                                     <  1:1459663/150‑1 (MQ=255)
                                                                                               gATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCaa                                                   >  2:358679/1‑150 (MQ=255)
                                                                                                aTTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAac                                                  >  2:1744067/1‑150 (MQ=255)
                                                                                                 tttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGcc                                                                >  1:2466011/1‑135 (MQ=255)
                                                                                                 tttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGcc                                                                <  2:2466011/135‑1 (MQ=255)
                                                                                                 tttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAacc                                                 <  2:665612/150‑1 (MQ=255)
                                                                                                    aaaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTTGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAATCCTGGCACAAAGCAaccacc                                              <  2:2260035/150‑1 (MQ=255)
                                                                                                    aaaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                              <  1:187815/150‑1 (MQ=255)
                                                                                                    aaaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                              >  2:2012425/1‑150 (MQ=255)
                                                                                                         cccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCCCCAACgg                                         >  1:700569/1‑150 (MQ=255)
                                                                                                          ccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                        <  2:901805/150‑1 (MQ=255)
                                                                                                              ccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGt                                    <  2:80753/150‑1 (MQ=255)
                                                                                                                gTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                  <  1:1965583/150‑1 (MQ=255)
                                                                                                                 tttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGACAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCgg                                 <  1:2192578/150‑1 (MQ=255)
                                                                                                                  ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGt                                <  2:329434/150‑1 (MQ=255)
                                                                                                                   tCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                               <  1:463841/150‑1 (MQ=255)
                                                                                                                     aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGAc                             <  2:643826/150‑1 (MQ=255)
                                                                                                                     aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGAc                             <  2:2145877/150‑1 (MQ=255)
                                                                                                                     aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGAc                             <  1:1247664/150‑1 (MQ=255)
                                                                                                                     aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCAGTGAc                             <  1:1729167/150‑1 (MQ=255)
                                                                                                                       cAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTc                           <  2:989471/150‑1 (MQ=255)
                                                                                                                         ggCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                         >  1:1005180/1‑150 (MQ=255)
                                                                                                                            acaACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACa                      <  2:1811948/150‑1 (MQ=255)
                                                                                                                             caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAGCATAAGCCTGGCACAAAGCAACCACCAACGGATGCTCGGTGACTCGTACAg                     <  2:1181715/150‑1 (MQ=255)
                                                                                                                             caACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAg                     <  1:823853/150‑1 (MQ=255)
                                                                                                                             caACCGGGAACGCCGCGCCAAATTGTGAGCGAACTTCCTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACAATAAGCCTGGCACAAACCAACCACCAACGCATGGTCGGTGACTCGTTCAg                     >  2:1659831/1‑150 (MQ=255)
                                                                                                                               aCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCa                   >  1:2004043/1‑150 (MQ=255)
                                                                                                                               aCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCa                   >  1:714370/1‑150 (MQ=255)
                                                                                                                                  gggAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGc                <  1:809706/150‑1 (MQ=255)
                                                                                                                                   cgaagttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCTGTGACTCGTACAGCAAGCg               >  1:2457486/8‑150 (MQ=255)
                                                                                                                                   cgaagttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCg               >  1:2361898/8‑150 (MQ=255)
                                                                                                                                    gaagttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGa              >  2:1519614/7‑150 (MQ=255)
                                                                                                                                      agttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTTGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAAt            >  2:1480078/5‑150 (MQ=255)
                                                                                                                                      agttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCGGATGGTCGGTGACTCGTACAGCAAGCGAAt            >  1:1801008/5‑150 (MQ=255)
                                                                                                                                       cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATc           >  1:2218121/1‑150 (MQ=255)
                                                                                                                                                aaaTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc  >  1:1820053/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CTGTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCC  >  NC_000913/3430798‑3431091

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