Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,144,158 C→A 85.7% 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,1580CA85.7% 573.8 / 95.6 231intergenic (+253/+19)yghA/exbDputative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base C (20/13);  new base A (122/76);  total (142/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 233 total aligned reads displayed.

TGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGG  >  CP009273/3144010‑3144304
                                                                                                                                                    |                                                                                                                                                  
      aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAg                                                                                                                                             >  1:1081011/1‑150 (MQ=255)
      aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAg                                                                                                                                             >  2:3487956/1‑150 (MQ=255)
          gATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGACCGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCaa                                                                                                                                         >  2:760393/1‑150 (MQ=255)
           aTTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAAt                                                                                                                                        >  2:3445267/1‑150 (MQ=255)
             tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTc                                                                                                                                      <  2:2059944/150‑1 (MQ=255)
              ttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCACCATGTTGCACTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCt                                                                                                                                     >  2:3531848/1‑150 (MQ=255)
              ttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCt                                                                                                                                     >  1:74223/1‑150 (MQ=255)
               tggaGGCCTGATACGCGGAGCGCATCAGGCATTGAGCACCTCATGCCGGCTGCGTCGTGCACGCCTTCTCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCta                                                                                                                                    <  2:4105966/147‑1 (MQ=25)
               tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCta                                                                                                                                    >  2:1824614/1‑150 (MQ=255)
                gTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCTACATGTGGCATTTGATTTTGTAGGGCTGATGAGCGGAGAGCATCAGGCCATTCtag                                                                                                                                   >  2:2120751/1‑149 (MQ=37)
                 tAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCACTTCtact                                                                                                                                  >  1:2105229/1‑150 (MQ=255)
                 tAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtact                                                                                                                                  <  2:2786607/150‑1 (MQ=255)
                  aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactt                                                                                                                                 >  1:3561071/1‑150 (MQ=255)
                  aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactt                                                                                                                                 >  2:3625074/1‑150 (MQ=255)
                  aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAGTCtactt                                                                                                                                 >  2:2920250/1‑150 (MQ=255)
                   ggCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttt                                                                                                                                <  2:1082457/150‑1 (MQ=255)
                    gCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttta                                                                                                                               <  2:914065/150‑1 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCACTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttac                                                                                                                              >  2:3223369/1‑150 (MQ=255)
                      cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttact                                                                                                                             >  1:408619/1‑150 (MQ=255)
                       tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCtactttactt                                                                                                                            >  1:2029177/1‑150 (MQ=255)
                       tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCtactttactt                                                                                                                            >  1:3423576/1‑150 (MQ=255)
                       tGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttactt                                                                                                                            <  2:1735528/150‑1 (MQ=255)
                         aTAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCg                                                                                                                          <  1:703200/150‑1 (MQ=255)
                            aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                       <  2:3920961/150‑1 (MQ=255)
                              cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTg                                                                                                                     >  2:1581628/1‑150 (MQ=255)
                                tAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGc                                                                                                                   >  2:516834/1‑150 (MQ=255)
                                 aGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCg                                                                                                                  >  1:3600312/1‑150 (MQ=255)
                                 aGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCg                                                                                                                  >  1:3810727/1‑150 (MQ=255)
                                     cATCAGGCATTGAGCACCTACTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGttt                                                                                                              <  1:622377/150‑1 (MQ=255)
                                          ggCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCg                                                                                                         >  2:2969567/1‑150 (MQ=255)
                                             aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGACCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCg                                                                                                      <  1:1896787/150‑1 (MQ=255)
                                              ttGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACCAACGCACGAAGATTCAACATGTTGCAATAGATGTTGTAGGCCTGATGAGCGTAGCGCATCAGGCCATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGa                                                                                                     >  2:3277421/1‑150 (MQ=255)
                                               tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGAc                                                                                                    >  2:3814617/1‑150 (MQ=255)
                                               tGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATGCAACATGGTGCCATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAGTCTACTTTACTTCGCTGTGGCGGTGTCTTCGCCGAc                                                                                                    >  2:3814511/1‑150 (MQ=37)
                                                  gCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAg                                                                                                 <  2:649537/150‑1 (MQ=255)
                                                     ccTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGAcc                                                                                              >  1:3214897/1‑150 (MQ=255)
                                                      cTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCt                                                                                             >  2:2978944/1‑150 (MQ=255)
                                                      cTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGAAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCTCCGACCAGACCt                                                                                             >  1:2115422/1‑150 (MQ=255)
                                                        aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTAt                                                                                           >  2:3926672/1‑150 (MQ=255)
                                                         aTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATc                                                                                          <  1:2575490/150‑1 (MQ=255)
                                                         aTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATc                                                                                          >  1:2293345/1‑150 (MQ=255)
                                                          tGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCt                                                                                         <  2:2546337/150‑1 (MQ=255)
                                                              ggATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAg                                                                                     >  2:198331/1‑150 (MQ=255)
                                                               gATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAgg                                                                                    <  1:55320/150‑1 (MQ=255)
                                                                aTGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGt                                                                                   >  2:4325928/1‑150 (MQ=255)
                                                                   cgtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                                >  2:1357742/1‑150 (MQ=255)
                                                                     tcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAccc                                                                              <  1:1012029/150‑1 (MQ=255)
                                                                     tcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAccc                                                                              <  1:232218/150‑1 (MQ=255)
                                                                      cgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCg                                                                             <  2:487060/150‑1 (MQ=255)
                                                                        tGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGcc                                                                           >  2:2786227/1‑150 (MQ=255)
                                                                           aCGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGa                                                                        >  2:1865377/1‑150 (MQ=255)
                                                                            cGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTATGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTAGTTCGCTTTGGCGGTTTCTTAGCCGACCAGACCTATATTCAGGTAACCCGCCAGAt                                                                       >  1:27749/1‑150 (MQ=25)
                                                                            cGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGTCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAt                                                                       >  1:27708/1‑150 (MQ=255)
                                                                              ccTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGc                                                                     <  1:943272/150‑1 (MQ=255)
                                                                                ttATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAg                                                                   >  1:3335481/1‑150 (MQ=255)
                                                                                  aTCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCg                                                                 >  1:2921165/1‑150 (MQ=255)
                                                                                   tCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGt                                                                >  1:417487/1‑150 (MQ=255)
                                                                                     cGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTAt                                                              >  2:3042086/1‑150 (MQ=255)
                                                                                      ggCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATc                                                             <  2:3636367/150‑1 (MQ=255)
                                                                                         caaCAAACGCACGAATATTCAACATGTTGCAATAGATTTTGTAGGCCTGTTGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAt                                                          <  1:1828201/148‑1 (MQ=255)
                                                                                            cAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCAGATGAGCGTAGAGAATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTAc                                                       >  2:1524284/1‑150 (MQ=255)
                                                                                             aaaCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTAcc                                                      >  1:1122096/1‑150 (MQ=255)
                                                                                              aaCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCt                                                     >  1:1483452/1‑150 (MQ=255)
                                                                                                   aCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatc                                                >  2:1139507/1‑150 (MQ=255)
                                                                                                   aCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatc                                                >  1:2529532/1‑150 (MQ=255)
                                                                                                   aCGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatc                                                >  2:3336424/1‑150 (MQ=255)
                                                                                                    cGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtcatca                                               <  1:4289543/150‑1 (MQ=255)
                                                                                                     gAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCaa                                              >  1:3251587/1‑150 (MQ=255)
                                                                                                           tCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtctc                                        <  2:1396455/150‑1 (MQ=255)
                                                                                                            cAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCg                                       <  1:725832/150‑1 (MQ=255)
                                                                                                              aCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTa                                     <  1:1993147/150‑1 (MQ=255)
                                                                                                              aCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTa                                     <  1:316612/150‑1 (MQ=255)
                                                                                                              aCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTa                                     <  1:222591/150‑1 (MQ=255)
                                                                                                                aTGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAAt                                   <  2:3689293/150‑1 (MQ=255)
                                                                                                                 tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                  <  1:4262140/150‑1 (MQ=255)
                                                                                                                 tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                  >  2:3094442/1‑150 (MQ=255)
                                                                                                                 tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                  <  1:3628953/150‑1 (MQ=255)
                                                                                                                  gTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCg                                 >  2:2537968/1‑150 (MQ=255)
                                                                                                                  gTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCg                                 >  1:2571524/1‑150 (MQ=255)
                                                                                                                        aTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtt                           >  1:3957794/1‑150 (MQ=255)
                                                                                                                        aTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtt                           <  2:4169752/150‑1 (MQ=255)
                                                                                                                          aGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                         <  1:2642226/150‑1 (MQ=255)
                                                                                                                          aGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                         >  2:2691662/1‑150 (MQ=255)
                                                                                                                          aGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                         >  2:2009369/1‑150 (MQ=255)
                                                                                                                            aTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACGTCGCTTTGGCGGTTTATTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTACTCGCCGGTTgtgt                       >  2:1467448/1‑150 (MQ=255)
                                                                                                                              tttGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTcc                     <  1:2162598/150‑1 (MQ=255)
                                                                                                                              tttGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTcc                     <  2:3999213/150‑1 (MQ=255)
                                                                                                                                 gTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCt                  >  1:1644392/1‑150 (MQ=255)
                                                                                                                                   aGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCg                >  2:1896564/1‑150 (MQ=255)
                                                                                                                                    ggCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCgg               >  1:916662/1‑150 (MQ=255)
                                                                                                                                    ggCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCgg               >  2:3140234/1‑150 (MQ=255)
                                                                                                                                      ccTGATGTGCGTAGAGCATAAGGAAATTCTACTTTACTTCGCTTTGGCGGTTTATTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCTTATCAATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaa             <  1:1595627/150‑1 (MQ=255)
                                                                                                                                      ccTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaa             >  2:1242809/1‑150 (MQ=255)
                                                                                                                                       cTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCCGCGTATCCATTACCTTCATCAAAGTCTCGTAATCGACGGTTTTGTCCGCTCGgaag            >  2:4251692/1‑150 (MQ=255)
                                                                                                                                        tgatgaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaaga           >  2:3665121/1‑150 (MQ=255)
                                                                                                                                          atgaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGaagaag         >  2:2397453/1‑150 (MQ=255)
                                                                                                                                            gaGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAt       <  2:2192180/150‑1 (MQ=255)
                                                                                                                                             aGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtg      <  2:1946466/150‑1 (MQ=255)
                                                                                                                                               cGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggt    >  1:457988/1‑150 (MQ=255)
                                                                                                                                                 tAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggtgg  <  2:1556292/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGG  >  CP009273/3144010‑3144304

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