Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,144,158 C→A 90.8% intergenic (+253/+19) yghA → / ← exbD putative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,144,1580CA90.8% 732.9 / 70.6 260intergenic (+253/+19)yghA/exbDputative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base C (21/3);  new base A (144/92);  total (165/95)
Fisher's exact test for biased strand distribution p-value = 1.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 260 total aligned reads displayed.

GCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG  >  CP009273/3144014‑3144303
                                                                                                                                                |                                                                                                                                                 
  aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAg                                                                                                                                            >  1:889795/1‑150 (MQ=255)
     aGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCa                                                                                                                                         >  1:4400081/1‑150 (MQ=255)
          ttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCt                                                                                                                                    <  2:2646873/150‑1 (MQ=255)
             tAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtact                                                                                                                                 <  2:1718459/150‑1 (MQ=255)
               ggCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttt                                                                                                                               <  1:677219/150‑1 (MQ=255)
               ggCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttt                                                                                                                               <  2:1895697/150‑1 (MQ=255)
                gCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttta                                                                                                                              <  1:3193494/150‑1 (MQ=255)
                 ccTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttac                                                                                                                             >  1:1551082/1‑150 (MQ=255)
                 ccTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttac                                                                                                                             >  2:1535206/1‑150 (MQ=255)
                 ccTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGAGTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttac                                                                                                                             >  2:4383373/1‑150 (MQ=255)
                   tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttactt                                                                                                                           >  1:2451234/1‑150 (MQ=255)
                    gATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTc                                                                                                                          <  2:838867/150‑1 (MQ=255)
                     aTAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCg                                                                                                                         >  1:1699416/1‑150 (MQ=255)
                      tAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGc                                                                                                                        >  2:1440257/1‑150 (MQ=255)
                      tAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGc                                                                                                                        >  1:2396243/1‑150 (MQ=255)
                       aaGCGTAGCGCATCAGGCATTGAGCCCCTAATGCCGGATGCGTCGTGAACGCCTCATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCt                                                                                                                       >  2:2888400/1‑150 (MQ=255)
                        aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                      >  2:3763403/1‑150 (MQ=255)
                           gTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCTACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGcgtggg                                                                                                                   >  2:4137002/1‑145 (MQ=255)
                             aGCGCATCAGTCAATGAGCAACTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACTCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGATTTGGCg                                                                                                                 <  2:4335786/150‑1 (MQ=25)
                             aGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAGTCTACTTTACTTCGCTTTGGCg                                                                                                                 >  2:1991734/1‑150 (MQ=255)
                                 cATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGttt                                                                                                             >  1:1006180/1‑150 (MQ=255)
                                   tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttct                                                                                                           <  2:64321/150‑1 (MQ=255)
                                   tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttct                                                                                                           <  2:1823541/150‑1 (MQ=255)
                                       gCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGc                                                                                                       >  1:20436/1‑150 (MQ=255)
                                        cATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGcc                                                                                                      <  1:4511944/150‑1 (MQ=255)
                                         aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCg                                                                                                     <  2:1140325/150‑1 (MQ=255)
                                          ttGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGa                                                                                                    <  2:4277320/150‑1 (MQ=255)
                                           tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAc                                                                                                   <  1:2646782/150‑1 (MQ=255)
                                           tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAc                                                                                                   >  1:783910/1‑150 (MQ=255)
                                           tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAc                                                                                                   >  1:1362414/1‑150 (MQ=255)
                                           tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAc                                                                                                   >  1:3682357/1‑150 (MQ=255)
                                           tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTGGGCGGTTTCTTCGCCGAc                                                                                                   >  2:317121/1‑150 (MQ=255)
                                               cACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGa                                                                                               >  1:3018369/1‑150 (MQ=255)
                                               cACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGa                                                                                               >  1:4048676/1‑150 (MQ=255)
                                                aCCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGAc                                                                                              <  1:1907114/150‑1 (MQ=255)
                                                  cTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCt                                                                                            <  1:2843166/150‑1 (MQ=255)
                                                  cTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCt                                                                                            <  1:2542736/150‑1 (MQ=255)
                                                    aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTAt                                                                                          >  1:753048/1‑150 (MQ=255)
                                                      tGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCt                                                                                        <  1:2993378/150‑1 (MQ=255)
                                                        ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTc                                                                                      >  2:4092784/1‑150 (MQ=255)
                                                        ccGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTc                                                                                      <  1:3551482/150‑1 (MQ=255)
                                                            aTGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGt                                                                                  >  1:2029348/1‑150 (MQ=255)
                                                               cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                               >  1:1082511/1‑150 (MQ=255)
                                                               cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                               >  2:2699378/1‑150 (MQ=255)
                                                               cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                               >  1:3712559/1‑150 (MQ=255)
                                                                  cgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCg                                                                            >  2:4551290/1‑150 (MQ=255)
                                                                   gtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGc                                                                           >  1:653008/1‑150 (MQ=255)
                                                                    tGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGcc                                                                          >  1:2141631/1‑150 (MQ=255)
                                                                      aaCGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTg                                                                        >  1:2336129/1‑150 (MQ=255)
                                                                       aCGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGa                                                                       >  1:3066973/1‑150 (MQ=255)
                                                                        cGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAt                                                                      <  1:3320909/150‑1 (MQ=255)
                                                                        cGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAt                                                                      <  1:1818069/150‑1 (MQ=255)
                                                                        cGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAt                                                                      <  1:3320929/150‑1 (MQ=255)
                                                                         gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATg                                                                     >  2:4231308/1‑150 (MQ=255)
                                                                           cTTATCCGGCCTACAAACGCACGAAGATTCAACCTGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCGTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCa                                                                   >  2:3125459/1‑150 (MQ=255)
                                                                           cTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCa                                                                   >  2:808183/1‑150 (MQ=255)
                                                                            ttATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAg                                                                  <  1:3451139/150‑1 (MQ=255)
                                                                            ttATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAg                                                                  <  2:244128/150‑1 (MQ=255)
                                                                                ccGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTa                                                              <  1:1374898/150‑1 (MQ=255)
                                                                                 cGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTAt                                                             >  2:1195007/1‑150 (MQ=255)
                                                                                  ggCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATc                                                            >  1:834425/1‑150 (MQ=255)
                                                                                  ggCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATc                                                            >  1:3230132/1‑150 (MQ=255)
                                                                                  ggCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATc                                                            >  2:2336693/1‑150 (MQ=255)
                                                                                     cTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAt                                                         >  1:2819336/1‑150 (MQ=255)
                                                                                     cTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAt                                                         >  2:3427764/1‑150 (MQ=255)
                                                                                        cAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTAc                                                      >  2:4605746/1‑150 (MQ=255)
                                                                                         aaaCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTAcc                                                     >  1:3984686/1‑150 (MQ=255)
                                                                                            cGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtc                                                  >  2:382380/1‑150 (MQ=255)
                                                                                              cACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcat                                                >  1:4477832/1‑150 (MQ=255)
                                                                                               aCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatc                                               >  1:3778232/1‑150 (MQ=255)
                                                                                               aCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCACTTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatc                                               >  2:2474190/1‑150 (MQ=255)
                                                                                                 gAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTATTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGTTGCAGCGTATCCATTACATTCAGCaa                                             >  2:2254610/1‑150 (MQ=38)
                                                                                                   aGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACg                                           >  1:4081466/1‑150 (MQ=255)
                                                                                                    gATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACgg                                          >  2:3525497/1‑149 (MQ=255)
                                                                                                     aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                         >  2:2993482/1‑150 (MQ=255)
                                                                                                     aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                         >  2:4614536/1‑150 (MQ=255)
                                                                                                        cAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCg                                      >  2:4521932/1‑150 (MQ=255)
                                                                                                         aaCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGt                                     >  1:1481898/1‑150 (MQ=255)
                                                                                                         aaCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGt                                     >  2:3492801/1‑150 (MQ=255)
                                                                                                           cATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTaa                                   >  2:2999689/1‑150 (MQ=255)
                                                                                                             tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                 <  2:2789983/150‑1 (MQ=255)
                                                                                                             tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                 >  1:1177001/1‑150 (MQ=255)
                                                                                                              gTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCg                                <  1:3921886/150‑1 (MQ=255)
                                                                                                              gTTGAAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGAAAATCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCTCCTGATGCATCGTATCCATTACCTTCATCAACGTCTCGTAATCg                                >  2:933768/1‑150 (MQ=255)
                                                                                                                tGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGAc                              <  1:4136937/150‑1 (MQ=255)
                                                                                                                tGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGAc                              <  2:3439889/150‑1 (MQ=255)
                                                                                                                 gCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACg                             <  1:1977935/150‑1 (MQ=255)
                                                                                                                   aaTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGt                           >  1:125346/1‑150 (MQ=255)
                                                                                                                    aTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtt                          >  1:160433/1‑150 (MQ=255)
                                                                                                                    aTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtt                          >  2:2533543/1‑150 (MQ=255)
                                                                                                                      aGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                        >  2:3561121/1‑150 (MQ=255)
                                                                                                                        aTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGt                      >  1:982997/1‑150 (MQ=255)
                                                                                                                        aTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGt                      >  1:524964/1‑150 (MQ=255)
                                                                                                                        aTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGt                      >  1:1991975/1‑150 (MQ=255)
                                                                                                                             gTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCt                 <  2:668463/150‑1 (MQ=255)
                                                                                                                                  ccTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaa            >  1:2841723/1‑150 (MQ=255)
                                                                                                                                  ccTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaa            >  2:2526776/1‑150 (MQ=255)
                                                                                                                                     gatgaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaagaa         >  2:569864/1‑150 (MQ=255)
                                                                                                                                         aGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtg     <  1:2724149/150‑1 (MQ=255)
                                                                                                                                            gTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggtg  >  1:2740342/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG  >  CP009273/3144014‑3144303

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: