Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,317,805 C→A 94.5% D156Y (GAT→TAT)  folP ← 7,8‑dihydropteroate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,317,8050CA94.5% 628.0 / 25.3 217D156Y (GAT→TAT) folP7,8‑dihydropteroate synthase
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base A (98/106);  total (105/112)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 217 total aligned reads displayed.

AGGCGCGCCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTT  >  CP009273/3317660‑3317949
                                                                                                                                                 |                                                                                                                                                
aGGCGCGCCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGt                                                                                                                                                >  1:3146008/1‑148 (MQ=255)
  gcgcgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCa                                                                                                                                              >  1:1420153/1‑148 (MQ=255)
  gcgcgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCa                                                                                                                                              >  2:1683189/1‑148 (MQ=255)
  gcgcgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCa                                                                                                                                              >  2:3084935/1‑148 (MQ=255)
   cgcgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCAt                                                                                                                                             <  1:1190080/148‑1 (MQ=255)
   cgcgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCAt                                                                                                                                             <  1:2941618/148‑1 (MQ=255)
     cgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATAc                                                                                                                                           <  2:1819129/148‑1 (MQ=255)
     cgcCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATAc                                                                                                                                           >  2:1950679/1‑148 (MQ=255)
       ccAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACtt                                                                                                                                         <  1:3032216/148‑1 (MQ=255)
         aGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCg                                                                                                                                       >  2:1848839/1‑148 (MQ=255)
          gTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCgg                                                                                                                                      <  2:361196/148‑1 (MQ=255)
           tAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGa                                                                                                                                     >  2:662520/1‑148 (MQ=255)
            aaTGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAg                                                                                                                                    <  2:2576635/148‑1 (MQ=255)
            aaTGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAg                                                                                                                                    >  1:1985420/1‑148 (MQ=255)
              tGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCt                                                                                                                                  >  1:585246/1‑148 (MQ=255)
                aaTAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTc                                                                                                                                <  2:1712844/148‑1 (MQ=255)
                 aTAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTcc                                                                                                                               <  1:886319/148‑1 (MQ=255)
                  tAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCt                                                                                                                              >  2:3177815/1‑148 (MQ=255)
                   aGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTg                                                                                                                             >  2:2651804/1‑148 (MQ=255)
                   aGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTg                                                                                                                             >  1:836645/1‑148 (MQ=255)
                     ttATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCa                                                                                                                           <  2:1899109/148‑1 (MQ=255)
                     ttATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCa                                                                                                                           <  2:3389454/148‑1 (MQ=255)
                     ttATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGAAATAGTCATACTTCGGAGCTTCCTGCa                                                                                                                           <  1:1067558/148‑1 (MQ=255)
                      tATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCAt                                                                                                                          >  1:259919/1‑148 (MQ=255)
                       aTGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATg                                                                                                                         <  2:2956607/148‑1 (MQ=255)
                        tGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATgg                                                                                                                        <  2:3343015/148‑1 (MQ=255)
                         gggAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGt                                                                                                                       <  2:818342/148‑1 (MQ=255)
                           gagagaTTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGttt                                                                                                                     <  2:1280695/148‑1 (MQ=255)
                           gagagaTTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGttt                                                                                                                     <  2:536944/148‑1 (MQ=255)
                             gagaTTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGttttt                                                                                                                   <  2:109125/148‑1 (MQ=255)
                                aTTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGa                                                                                                                <  2:153901/148‑1 (MQ=255)
                                aTTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGa                                                                                                                <  1:1577444/148‑1 (MQ=255)
                                 tttttACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGAt                                                                                                               <  2:39334/148‑1 (MQ=255)
                                  ttttACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGAtt                                                                                                              >  2:1962730/1‑148 (MQ=255)
                                  ttttACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGAtt                                                                                                              >  2:2112165/1‑148 (MQ=255)
                                   tttACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGAttt                                                                                                             <  2:915686/148‑1 (MQ=255)
                                      accgaaaccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTccc                                                                                                          <  2:21711/148‑1 (MQ=255)
                                       ccgaaaccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCt                                                                                                         >  2:105336/1‑148 (MQ=255)
                                       ccgaaaccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCt                                                                                                         >  1:1737661/1‑148 (MQ=255)
                                         gaaaccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGc                                                                                                       >  1:285897/1‑148 (MQ=255)
                                         gaaaccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGc                                                                                                       <  2:676308/148‑1 (MQ=255)
                                           aaccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCat                                                                                                     <  1:151454/148‑1 (MQ=255)
                                             ccgaaTCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCatat                                                                                                   <  1:1607406/148‑1 (MQ=255)
                                                    cGGGGTCGAGCAACAATTTCTCTTTTGCGAGACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAg                                                                                            <  2:948633/148‑1 (MQ=255)
                                                      gggTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGAc                                                                                          >  2:649100/1‑148 (MQ=255)
                                                       ggTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACa                                                                                         <  2:2842288/148‑1 (MQ=255)
                                                                acaaTTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaa                                                                                >  1:3300189/1‑148 (MQ=255)
                                                                acaaTTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaa                                                                                <  1:3403805/148‑1 (MQ=255)
                                                                acaaTTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaa                                                                                <  1:1666152/148‑1 (MQ=255)
                                                                 caaTTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaaa                                                                               >  2:970446/1‑148 (MQ=255)
                                                                 caaTTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaaa                                                                               <  2:1872288/148‑1 (MQ=255)
                                                                 caaTTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaaa                                                                               >  1:1605845/1‑148 (MQ=255)
                                                                    tttCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaaaccg                                                                            >  1:812818/1‑148 (MQ=255)
                                                                     ttCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaaaccgg                                                                           <  1:946958/148‑1 (MQ=255)
                                                                      tctcTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACaaaccggtaaaccggt                                                                          <  2:1892869/148‑1 (MQ=255)
                                                                       ctctTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGtt                                                                         >  2:2259158/1‑148 (MQ=255)
                                                                           tttGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTg                                                                     >  1:87735/1‑148 (MQ=255)
                                                                               cGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAgc                                                                 <  1:927637/148‑1 (MQ=255)
                                                                                 aTACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTCGCGATTCACTTCTGCAAAGACATAGTCATCCTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAgccg                                                               >  2:3031784/1‑148 (MQ=255)
                                                                                      cGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTcc                                                          <  1:2554775/148‑1 (MQ=255)
                                                                                                  aCGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGt                                              <  2:3220398/148‑1 (MQ=255)
                                                                                                   cGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGtt                                             >  2:3495192/1‑148 (MQ=255)
                                                                                                      gCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTACACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCgg                                          >  2:3126961/1‑148 (MQ=255)
                                                                                                        tATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGaa                                        <  1:2484916/148‑1 (MQ=255)
                                                                                                        tATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGaa                                        <  1:3166278/148‑1 (MQ=255)
                                                                                                         aTTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGaaa                                       <  1:484940/148‑1 (MQ=255)
                                                                                                         aTTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGaaa                                       >  2:2597340/1‑148 (MQ=255)
                                                                                                         aTTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGaaa                                       <  2:2073367/148‑1 (MQ=255)
                                                                                                          tttGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAg                                      >  1:2416397/1‑148 (MQ=255)
                                                                                                          tttGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAg                                      >  1:2416488/1‑148 (MQ=255)
                                                                                                          tttGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTACACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAg                                      >  1:960163/1‑148 (MQ=255)
                                                                                                          tttGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCGGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAg                                      <  2:3151335/148‑1 (MQ=255)
                                                                                                           ttGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAgg                                     <  2:1420153/148‑1 (MQ=255)
                                                                                                           ttGCTCAATAAAGTAGCGATTCACTTCTGAAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAgg                                     <  2:1990708/148‑1 (MQ=255)
                                                                                                                cAATAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCg                                >  2:237408/1‑148 (MQ=255)
                                                                                                                  aTAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGa                              <  2:1074722/148‑1 (MQ=255)
                                                                                                                   tAAAGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGat                             <  1:1383244/148‑1 (MQ=255)
                                                                                                                    aaaGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGata                            >  2:854967/1‑148 (MQ=255)
                                                                                                                      aGTAGCGATTCTCTTCTTCAAAGACATTGTCATACTTCTGAGGTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTACAGAGCGCCAGGTTCGGAAAGGGAGCGGATATc                          <  1:2725616/148‑1 (MQ=255)
                                                                                                                      aGTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATc                          <  1:2922315/148‑1 (MQ=255)
                                                                                                                       gTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCa                         <  1:1848864/148‑1 (MQ=255)
                                                                                                                       gTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCa                         >  2:2889202/1‑148 (MQ=255)
                                                                                                                       gTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCa                         >  1:580197/1‑148 (MQ=255)
                                                                                                                       gTAGCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCa                         >  1:2885350/1‑148 (MQ=255)
                                                                                                                          gCGATTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATta                      >  1:2040869/1‑148 (MQ=255)
                                                                                                                             aTTCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATtaata                   <  1:631179/148‑1 (MQ=255)
                                                                                                                              ttCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATtaataa                  >  1:28234/1‑148 (MQ=255)
                                                                                                                              ttCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATtaataa                  >  2:356405/1‑148 (MQ=255)
                                                                                                                              ttCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATtaataa                  >  2:2291182/1‑148 (MQ=255)
                                                                                                                              ttCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATtaataa                  >  2:232384/1‑148 (MQ=255)
                                                                                                                               tCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAATGGGAGCGGATATCATTaataat                 >  2:74838/1‑148 (MQ=255)
                                                                                                                               tCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTaataat                 >  2:629498/1‑148 (MQ=255)
                                                                                                                               tCACTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTaataat                 >  1:1312353/1‑148 (MQ=255)
                                                                                                                                 aCTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTACACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAAtgt               >  2:3251238/1‑148 (MQ=255)
                                                                                                                                  cTTCTGCAAAGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAAtgtg              >  2:1050880/1‑148 (MQ=255)
                                                                                                                                          aaGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCaa      >  2:24632/1‑148 (MQ=255)
                                                                                                                                          aaGACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCaa      >  2:2675900/1‑148 (MQ=255)
                                                                                                                                            gACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACt    >  1:3067717/1‑148 (MQ=255)
                                                                                                                                            gACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACt    <  1:728407/148‑1 (MQ=255)
                                                                                                                                            gACATAGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACt    <  1:100664/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGCGCGCCAGTAATGAATAGTTATGGGAGAGATTTTTACCGAAACCGAATCCGGGGTCGAGCAACAATTTCTCTTTTGCGATACCCGCCTGCTCGCAACGTGCTATTTGCTCAATAAAGTAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTT  >  CP009273/3317660‑3317949

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