Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,984,342 G→T 45.6% intergenic (‑216/‑171) hemC ← / → cyaA hydroxymethylbilane synthase/adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,984,3420GT45.6% 239.2 / inf 333intergenic (‑216/‑171)hemC/cyaAhydroxymethylbilane synthase/adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (83/98);  new base T (53/98);  total (137/196)
Fisher's exact test for biased strand distribution p-value = 5.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 340 total aligned reads displayed.

GGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAC  >  CP009273/3984193‑3984488
                                                                                                                                                     |                                                                                                                                                  
  tgAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATgt                                                                                                                                                  <  1:4540825/150‑1 (MQ=255)
    aaaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATgtgt                                                                                                                                                >  2:1921970/1‑150 (MQ=255)
    aaaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGgtgt                                                                                                                                                >  2:2016077/1‑150 (MQ=255)
     aaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGtt                                                                                                                                               >  2:4527437/1‑150 (MQ=255)
     aaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGtt                                                                                                                                               <  1:3022665/150‑1 (MQ=255)
     aaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGtt                                                                                                                                               <  2:2727506/150‑1 (MQ=255)
       cGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTaa                                                                                                                                             <  2:4391835/150‑1 (MQ=255)
       cGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTaa                                                                                                                                             <  1:1805264/150‑1 (MQ=255)
        gTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTaaa                                                                                                                                            >  1:565759/1‑150 (MQ=255)
         tAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAAt                                                                                                                                           <  2:4280450/150‑1 (MQ=255)
             cGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTTGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGAt                                                                                                                                       <  2:756982/150‑1 (MQ=255)
              gCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATc                                                                                                                                      <  2:4170132/150‑1 (MQ=255)
               ccTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCa                                                                                                                                     <  1:2513370/150‑1 (MQ=255)
                cTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCAc                                                                                                                                    <  1:3727294/150‑1 (MQ=255)
                 tGCGCAAACCATAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACg                                                                                                                                   <  1:4324655/150‑1 (MQ=255)
                    gcAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGttt                                                                                                                                <  1:1119333/150‑1 (MQ=255)
                    gcAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTATCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACgagt                                                                                                                                >  1:1834936/1‑147 (MQ=255)
                      aaaCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTa                                                                                                                              >  2:321623/1‑150 (MQ=255)
                      aaaCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTa                                                                                                                              >  2:321542/1‑150 (MQ=255)
                      aaaCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTa                                                                                                                              <  1:2792661/150‑1 (MQ=255)
                        aCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGa                                                                                                                            <  2:42998/150‑1 (MQ=255)
                         ccGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGAc                                                                                                                           <  1:2308018/150‑1 (MQ=255)
                          cGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGAcc                                                                                                                          >  1:1514002/1‑150 (MQ=255)
                           gTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCa                                                                                                                         <  2:3620142/150‑1 (MQ=255)
                           gTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCa                                                                                                                         >  1:710044/1‑150 (MQ=255)
                           gTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCa                                                                                                                         <  2:3358807/150‑1 (MQ=255)
                               aaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAtttt                                                                                                                     <  2:2244235/150‑1 (MQ=255)
                               aaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAtttt                                                                                                                     <  2:2242944/150‑1 (MQ=255)
                                  gAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtc                                                                                                                  >  1:1941269/1‑150 (MQ=255)
                                    ggTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtcgt                                                                                                                >  1:2587842/1‑150 (MQ=255)
                                       cTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTcgtcgt                                                                                                             >  2:2344424/1‑150 (MQ=255)
                                          gCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGaa                                                                                                          <  2:3721630/150‑1 (MQ=255)
                                          gCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGaa                                                                                                          <  2:1457120/150‑1 (MQ=255)
                                          gCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGaa                                                                                                          <  2:2758803/150‑1 (MQ=255)
                                             gTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACt                                                                                                       >  2:3518689/1‑150 (MQ=255)
                                              tGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTa                                                                                                      <  1:897461/150‑1 (MQ=255)
                                               ggATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaa                                                                                                     >  1:66855/1‑150 (MQ=255)
                                               ggATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaa                                                                                                     <  1:526947/150‑1 (MQ=255)
                                                        cAGGCTTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAg                                                                                            <  2:3060328/150‑1 (MQ=255)
                                                        cAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAg                                                                                            <  2:3961128/150‑1 (MQ=255)
                                                        cAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAg                                                                                            >  2:1115787/1‑150 (MQ=255)
                                                           gCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGgcg                                                                                         >  2:1305630/1‑150 (MQ=255)
                                                               ttCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGaa                                                                                     <  2:886323/150‑1 (MQ=255)
                                                                tCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGaaa                                                                                    <  1:3926845/150‑1 (MQ=255)
                                                                   cGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGt                                                                                 >  2:2331693/1‑150 (MQ=255)
                                                                     ccGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTgg                                                                               <  1:4317561/150‑1 (MQ=255)
                                                                     ccGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTgg                                                                               <  1:3249890/150‑1 (MQ=255)
                                                                       gTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTa                                                                             <  2:3577098/150‑1 (MQ=255)
                                                                         tGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAAc                                                                           <  2:258109/150‑1 (MQ=255)
                                                                         tGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAAc                                                                           <  2:4354675/150‑1 (MQ=255)
                                                                             ataaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGtt                                                                       <  1:3165179/150‑1 (MQ=255)
                                                                                   aaTTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCttt                                                                 <  1:593957/150‑1 (MQ=255)
                                                                                      ttACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGAc                                                              >  1:3547939/1‑150 (MQ=255)
                                                                                         cAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATa                                                           <  1:1049666/150‑1 (MQ=255)
                                                                                         cAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATa                                                           >  1:2102597/1‑150 (MQ=255)
                                                                                         cAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATa                                                           >  1:4088169/1‑150 (MQ=255)
                                                                                           gagaATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTGACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACg                                                         <  2:2848849/150‑1 (MQ=255)
                                                                                               aTAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTTAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAat                                                     <  2:409818/150‑1 (MQ=255)
                                                                                                 aaaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGACACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                                   <  1:2937907/150‑1 (MQ=255)
                                                                                                 aaaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                                   <  1:1701816/150‑1 (MQ=255)
                                                                                                 aaaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                                   <  1:876704/150‑1 (MQ=255)
                                                                                                     ggTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCg                                               <  1:4193234/150‑1 (MQ=255)
                                                                                                      gTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGa                                              <  2:2076858/150‑1 (MQ=255)
                                                                                                      gTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGa                                              >  1:1078700/1‑150 (MQ=255)
                                                                                                       tGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGaa                                             >  2:154405/1‑150 (MQ=255)
                                                                                                       tGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGaa                                             >  1:4379767/1‑150 (MQ=255)
                                                                                                         cTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATg                                           >  1:2132210/1‑150 (MQ=255)
                                                                                                            cacTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcg                                        <  1:2264994/150‑1 (MQ=255)
                                                                                                             acTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcgc                                       <  2:266024/150‑1 (MQ=255)
                                                                                                               ttGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                                     <  2:1361136/150‑1 (MQ=255)
                                                                                                                    tgtaGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAc                                <  2:2223414/150‑1 (MQ=255)
                                                                                                                    tgtaGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAc                                >  1:1251119/1‑150 (MQ=255)
                                                                                                                    tgtaGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAc                                >  1:1250873/1‑150 (MQ=255)
                                                                                                                     gtaGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAcc                               <  1:3274225/150‑1 (MQ=255)
                                                                                                                      taGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCg                              >  1:1422281/1‑150 (MQ=255)
                                                                                                                      taGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCg                              >  1:3689671/1‑150 (MQ=255)
                                                                                                                       aGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGt                             <  1:3577546/150‑1 (MQ=255)
                                                                                                                           cATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGag                         >  2:2928094/1‑149 (MQ=255)
                                                                                                                            aTCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATg                        <  2:3881246/150‑1 (MQ=255)
                                                                                                                            aTCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATg                        >  1:3384653/1‑150 (MQ=255)
                                                                                                                                ttctttACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTgg                    >  1:2189615/1‑150 (MQ=255)
                                                                                                                                ttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTgg                    <  2:867312/150‑1 (MQ=255)
                                                                                                                                ttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTgg                    <  1:3467416/150‑1 (MQ=255)
                                                                                                                                ttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTgg                    <  2:4290108/150‑1 (MQ=255)
                                                                                                                                     tACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGaa               >  1:1553264/1‑150 (MQ=255)
                                                                                                                                       cGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATc             <  1:4368739/150‑1 (MQ=255)
                                                                                                                                        ggTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATca            <  1:1656210/150‑1 (MQ=255)
                                                                                                                                          tCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATcaca          <  2:2537061/150‑1 (MQ=255)
                                                                                                                                            aaTCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGt        <  2:710044/150‑1 (MQ=255)
                                                                                                                                            aaTCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGt        <  2:1172995/150‑1 (MQ=255)
                                                                                                                                              tCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa      <  2:2029541/150‑1 (MQ=255)
                                                                                                                                              tCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa      >  2:2990646/1‑150 (MQ=255)
                                                                                                                                              tCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa      >  1:265800/1‑150 (MQ=255)
                                                                                                                                              tCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa      >  2:534743/1‑150 (MQ=255)
                                                                                                                                              tCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa      >  2:535290/1‑150 (MQ=255)
                                                                                                                                              tCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa      <  1:2045850/150‑1 (MQ=255)
                                                                                                                                                aGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATg    >  2:4418226/1‑150 (MQ=255)
                                                                                                                                                  cAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAc  <  1:1491728/150‑1 (MQ=255)
                                                                                                                                                  cAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAc  <  2:4365165/150‑1 (MQ=255)
                                                                                                                                                  cAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAc  <  2:3070922/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAC  >  CP009273/3984193‑3984488

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: