Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,430,948 T→G 8.2% A156A (GCA→GCC tsaC ← threonylcarbamoyl‑AMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,430,9480TG 8.2% 382.7 / 24.7 134A156A (GCA→GCCtsaCthreonylcarbamoyl‑AMP synthase
Reads supporting (aligned to +/- strand):  ref base T (59/64);  new base G (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 = 9.18e-01

Only 100 of 137 total aligned reads displayed.

GTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCC  >  NC_000913/3430800‑3431091
                                                                                                                                                    |                                                                                                                                               
gtgGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATtgtg                                                                                                                                                >  2:1094565/1‑150 (MQ=255)
 tgGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTgc                                                                                                                                               <  1:1207158/150‑1 (MQ=255)
   ggCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTgcgc                                                                                                                                             <  2:1782230/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:553936/1‑150 (MQ=255)
                ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                >  1:2203275/1‑150 (MQ=255)
                ccAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTAc                                                                                                                                >  1:901805/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)
                                tttCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGc                                                                                                                <  2:1454769/150‑1 (MQ=255)
                                  tCCATTATGTTACCCCTGTCGAAACAGTTCACCCTTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCaa                                                                                                              <  1:2368349/150‑1 (MQ=255)
                                   ccATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAAt                                                                                                             >  1:1547860/1‑150 (MQ=255)
                                   ccATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAAt                                                                                                             >  2:1247664/1‑150 (MQ=255)
                                   ccATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAAt                                                                                                             >  2:809706/1‑150 (MQ=255)
                                    cATTATGTTACACCTGTCGAAACAGTTCACCAGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTATCACCAGGCACAACCGGCAACGCCGCGCCAAATAGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGCCAAta                                                                                                            >  2:630048/1‑149 (MQ=255)
                                     aTTATGTTACCCCTGTCGAAACAGTTAACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGCCACAACCGGGAACGCCGAGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAAAcc                                                                                                           >  2:1073171/1‑150 (MQ=255)
                                         tGTTACCCCTGTCGAAAGAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTc                                                                                                       >  2:2156029/1‑150 (MQ=255)
                                            tACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCGGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGAAAAGGTGCAAATCAACTCAAg                                                                                                    >  2:425583/1‑150 (MQ=255)
                                              ccccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGtt                                                                                                  >  1:80753/1‑150 (MQ=255)
                                               cccTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTg                                                                                                 <  1:1956451/150‑1 (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:523473/150‑1 (MQ=255)
                                                   gTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGCTGGCAATCCACTCAAGTTGGCAc                                                                                             >  1:2399027/1‑150 (MQ=255)
                                                    tCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAAGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGAc                                                                                                                        <  2:1729218/122‑1 (MQ=255)
                                                     cGAAACAGTTCCCCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTGTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTg                                                                                           >  2:2460484/1‑150 (MQ=255)
                                                     cGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTg                                                                                           >  2:164099/1‑150 (MQ=255)
                                                     cGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTg                                                                                           >  1:1556383/1‑150 (MQ=255)
                                                       aaaCAGTTCCCCCGTCAGGGCATCGCGGCTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACCACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCCCTCAAGTTGGCACTGGt                                                                                         >  2:2152383/1‑150 (MQ=255)
                                                             ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   >  1:1263064/1‑150 (MQ=255)
                                                             ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   <  2:1962109/150‑1 (MQ=255)
                                                             ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAAc                                                                                   >  1:2292878/1‑150 (MQ=255)
                                                               cACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCa                                                                                 >  1:1390120/1‑150 (MQ=255)
                                                                aCCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAg                                                                                <  1:2031334/150‑1 (MQ=255)
                                                                aCCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAg                                                                                <  2:1306487/150‑1 (MQ=255)
                                                                   cGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCTAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCgg                                                                             <  1:669285/150‑1 (MQ=255)
                                                                   cGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCgg                                                                             <  1:414735/150‑1 (MQ=255)
                                                                    gTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTTGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGt                                                                            >  1:2260035/1‑150 (MQ=255)
                                                                       aGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTa                                                                         >  1:1793213/1‑150 (MQ=255)
                                                                        gggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAc                                                                        >  2:1837370/1‑150 (MQ=255)
                                                                        gggCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTAc                                                                        >  1:1727690/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)
                                                                                  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)
                                                                                      tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcac                                                          <  1:337856/150‑1 (MQ=255)
                                                                                      tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcac                                                          >  2:2373580/1‑150 (MQ=255)
                                                                                       cTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGTGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                         >  2:2486184/1‑150 (MQ=255)
                                                                                         gAAGGATTTAAACGCCCCCCGGTATCACCAGGCACAACCGGGACCGCCGGCACAAATTGTGCCCGAAATACGTCAACTGTACGACAAGGTGGCAAACCACTCAAGTGGGCACTGGTAGAAAGATGCGGTTTACCATAAGCCTCCCACaaa                                                       >  2:343247/1‑150 (MQ=255)
                                                                                         gAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACaaa                                                       >  1:670187/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)
                                                                                               tttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAacc                                                 <  2:665612/150‑1 (MQ=255)
                                                                                                  aaaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTTGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAATCCTGGCACAAAGCAaccacc                                              <  2:2260035/150‑1 (MQ=255)
                                                                                                        ccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                        >  1:2209636/1‑150 (MQ=255)
                                                                                                        ccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                        <  2:901805/150‑1 (MQ=255)
                                                                                                            ccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGt                                    <  2:80753/150‑1 (MQ=255)
                                                                                                                 tCACCAGGCACAACCGTGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                               <  1:2486184/150‑1 (MQ=255)
                                                                                                                 tCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                               <  1:35483/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:332714/150‑1 (MQ=255)
                                                                                                                     cAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTc                           <  2:989471/150‑1 (MQ=255)
                                                                                                                      aGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                          <  1:318299/150‑1 (MQ=255)
                                                                                                                       ggCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                         >  1:1005180/1‑150 (MQ=255)
                                                                                                                       ggCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGt                         >  1:887867/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:714370/1‑150 (MQ=255)
                                                                                                                             aCCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCa                   >  1:2004043/1‑150 (MQ=255)
                                                                                                                                  gaagttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGa              >  2:1519614/7‑150 (MQ=255)
                                                                                                                                    agttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTTGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAAt            >  2:1480078/5‑150 (MQ=255)
                                                                                                                                    agttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCGGATGGTCGGTGACTCGTACAGCAAGCGAAt            >  1:1801008/5‑150 (MQ=255)
                                                                                                                                     gttcgcgCACAATTTGGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATc           >  2:247425/4‑150 (MQ=255)
                                                                                                                                     cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATc           >  1:2218121/1‑150 (MQ=255)
                                                                                                                                          cgcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGc      <  1:538919/150‑1 (MQ=255)
                                                                                                                                          cgcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGc      <  1:424084/150‑1 (MQ=255)
                                                                                                                                          cgcCAAATTGTGCGCGAACTTCGTCTACTGGTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCAGTGACTCGTACAGCAAGCGAATCAAAGc      <  2:2053991/150‑1 (MQ=255)
                                                                                                                                              aaaTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc  >  1:1820053/1‑150 (MQ=255)
                                                                                                                                                    tgCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATc           <  1:1133572/135‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTGGGCTATCGGATTACCAAAAACAGCATAGGTTTCCATTATGTTACCCCTGTCGAAACAGTTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCC  >  NC_000913/3430800‑3431091

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