Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 80,123 C→T 93.0% G25S (GGC→AGC)  leuL ← leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927380,1230CT93.0% 935.2 / 45.3 270G25S (GGC→AGC) leuLleu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (7/12);  new base T (125/126);  total (132/138)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 272 total aligned reads displayed.

TCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCA  >  CP009273/79976‑80268
                                                                                                                                                   |                                                                                                                                                 
tCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAtgctgc                                                                                                                                                 >  2:3802866/1‑150 (MQ=255)
     cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACt                                                                                                                                            >  2:3087684/1‑150 (MQ=255)
     cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACt                                                                                                                                            >  2:4239743/1‑150 (MQ=255)
      aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTc                                                                                                                                           >  1:2885678/1‑150 (MQ=255)
      aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTc                                                                                                                                           <  2:2789901/150‑1 (MQ=255)
          tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCt                                                                                                                                       >  1:3822737/1‑150 (MQ=255)
          tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCt                                                                                                                                       <  1:4054883/150‑1 (MQ=255)
          tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCt                                                                                                                                       >  2:2611978/1‑150 (MQ=255)
           ggTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTa                                                                                                                                      >  1:3799136/1‑150 (MQ=255)
              aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCACTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcg                                                                                                                                   <  2:1664180/150‑1 (MQ=255)
              aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcg                                                                                                                                   <  1:3462984/150‑1 (MQ=255)
              aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcg                                                                                                                                   >  2:149321/1‑150 (MQ=255)
                cgcaaaTCATGACGTGCTGGCTCATGGTGTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTCATGCTGGATGCTGCTCACTCGGCTCCcgcg                                                                                                                                 <  2:2079770/147‑1 (MQ=255)
                  aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCa                                                                                                                               >  1:1900027/1‑150 (MQ=255)
                  aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCa                                                                                                                               <  1:3125209/150‑1 (MQ=255)
                   aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCaa                                                                                                                              >  1:4197270/1‑150 (MQ=255)
                   aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCaa                                                                                                                              >  1:1606170/1‑150 (MQ=255)
                     ataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAaag                                                                                                                            <  1:3880100/150‑1 (MQ=255)
                      taatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAaaga                                                                                                                           <  2:2267783/150‑1 (MQ=255)
                          gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGAt                                                                                                                       >  1:5122577/1‑150 (MQ=255)
                           aCTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATg                                                                                                                      >  2:2742132/1‑150 (MQ=255)
                            cgtgcttgctGATGGTTTGGGTGCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAACAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGTAGATGc                                                                                                                     <  2:4518073/148‑1 (MQ=255)
                             ttGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCg                                                                                                                    >  2:1112075/1‑150 (MQ=255)
                             ttGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCg                                                                                                                    >  1:129794/1‑150 (MQ=255)
                              tGCTGGCTCTTGTTTTGTGTCCTTGTCTCTTTTAGAGCGCCTCGATTCTGGGATAAAAAAACCCTCGCAATGGCGCGTGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGGTGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                                                                                                                   <  2:581834/150‑1 (MQ=255)
                              tGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                                                                                                                   >  1:3182899/1‑150 (MQ=255)
                                cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCATCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAAAGATGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                                                                 >  2:2294450/1‑150 (MQ=255)
                                 tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                                                                                                                <  1:3286786/150‑1 (MQ=255)
                                      cagggTTTGGGTCTTGGCCCTTTTTCTTGCGCCCCGCTCCGGGAATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagt                                                                                                           <  2:3572810/147‑1 (MQ=255)
                                          gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagta                                                                                                       <  2:1749968/150‑1 (MQ=255)
                                               ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCg                                                                                                  >  1:4967960/1‑150 (MQ=255)
                                                 tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGAt                                                                                                >  1:1392220/1‑150 (MQ=255)
                                                  ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATa                                                                                               >  2:3254453/1‑150 (MQ=255)
                                                     tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                                                                            <  1:342037/150‑1 (MQ=255)
                                                     tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                                                                            <  2:31252/150‑1 (MQ=255)
                                                      gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                                                                           <  2:285440/150‑1 (MQ=255)
                                                       tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCg                                                                                          <  1:2128514/150‑1 (MQ=255)
                                                          ctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                                                                       >  2:1228139/1‑150 (MQ=255)
                                                           ttttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACg                                                                                      >  2:3266606/1‑150 (MQ=255)
                                                             ttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                                                                    >  2:2970226/1‑150 (MQ=255)
                                                             ttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                                                                    >  2:2119097/1‑150 (MQ=255)
                                                              tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                                                                   >  1:2462818/1‑150 (MQ=255)
                                                              tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                                                                   >  2:5112719/1‑150 (MQ=255)
                                                               agagCGCCTCGCTTCGGGCATAAAAAATCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                                                                  >  1:4805142/1‑150 (MQ=255)
                                                                 agCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                                                                <  2:2571718/150‑1 (MQ=255)
                                                                    gcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTc                                                                             <  2:1204871/150‑1 (MQ=255)
                                                                     ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                                                                            <  1:3899399/150‑1 (MQ=255)
                                                                     ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                                                                            <  1:472138/150‑1 (MQ=255)
                                                                      cTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAt                                                                           >  2:3326559/1‑150 (MQ=255)
                                                                      cTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAt                                                                           >  2:4486992/1‑150 (MQ=255)
                                                                      cTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAt                                                                           <  1:1890622/150‑1 (MQ=255)
                                                                        cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                                                                         >  2:3817489/1‑150 (MQ=255)
                                                                         gCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTaa                                                                        >  2:4882520/1‑150 (MQ=255)
                                                                             cGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCa                                                                    <  1:2909304/150‑1 (MQ=255)
                                                                             cGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCa                                                                    >  1:3572033/1‑150 (MQ=255)
                                                                              gggCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                                                                   >  1:1484265/1‑150 (MQ=255)
                                                                               ggCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGc                                                                  >  1:3711643/1‑150 (MQ=255)
                                                                                  aTAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTcc                                                               >  2:3838127/1‑150 (MQ=255)
                                                                                    aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg                                                             <  1:1062545/150‑1 (MQ=255)
                                                                                       aaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg                                                          >  1:4532951/1‑150 (MQ=255)
                                                                                          aCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAAt                                                       <  2:602198/150‑1 (MQ=255)
                                                                                          aCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAAt                                                       <  2:3151180/150‑1 (MQ=255)
                                                                                              gcgcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGa                                                   <  2:779300/150‑1 (MQ=255)
                                                                                               cgcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGat                                                  <  1:2690309/150‑1 (MQ=255)
                                                                                               cgcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGat                                                  <  2:2385904/150‑1 (MQ=255)
                                                                                                gccaTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGata                                                 <  2:1835353/147‑1 (MQ=255)
                                                                                                 cAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGatat                                                >  2:4114049/1‑150 (MQ=255)
                                                                                                  aaTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATg                                               <  1:2487906/150‑1 (MQ=255)
                                                                                                  aaTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATg                                               <  1:944267/150‑1 (MQ=255)
                                                                                                   agggCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGc                                              <  1:4342341/148‑1 (MQ=255)
                                                                                                    tGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCt                                             <  2:4676040/150‑1 (MQ=255)
                                                                                                     ggCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCtt                                            <  1:3617022/150‑1 (MQ=255)
                                                                                                       cgcgGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCtttt                                          <  2:4725122/150‑1 (MQ=255)
                                                                                                        gcgGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTa                                         >  1:2238046/1‑150 (MQ=255)
                                                                                                           ggTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTagag                                      >  1:1342225/1‑150 (MQ=255)
                                                                                                             ttttttGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGtt                                    <  1:2507915/150‑1 (MQ=255)
                                                                                                              tttttGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTa                                   >  1:3994949/1‑150 (MQ=255)
                                                                                                                tttgtttgACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACt                                 >  2:4147757/1‑150 (MQ=255)
                                                                                                                   gtttgACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGa                              <  1:1815770/150‑1 (MQ=255)
                                                                                                                     ttgACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATa                            <  2:1367829/150‑1 (MQ=255)
                                                                                                                     ttgACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATa                            <  2:195738/150‑1 (MQ=255)
                                                                                                                       gACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACa                          >  1:4421613/1‑150 (MQ=255)
                                                                                                                         cTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaa                        >  1:196178/1‑150 (MQ=255)
                                                                                                                          tGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaaa                       >  2:2194917/1‑150 (MQ=255)
                                                                                                                            cGTGCTGGCTTAATGCTGGTTGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAAc                     <  2:4047983/150‑1 (MQ=255)
                                                                                                                             gTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACg                    <  2:621841/150‑1 (MQ=255)
                                                                                                                                cgggCTTAATGCTGGATGCTGCTCACTCGTCTACCGCTCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGAt                 <  2:3381236/148‑1 (MQ=255)
                                                                                                                                 tGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATg                >  2:2308854/1‑150 (MQ=255)
                                                                                                                                  ggCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGt               >  1:4647115/1‑150 (MQ=255)
                                                                                                                                   gCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTc              >  1:3022928/1‑150 (MQ=255)
                                                                                                                                   gCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTc              >  1:1762919/1‑150 (MQ=255)
                                                                                                                                   gCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTc              >  2:252699/1‑150 (MQ=255)
                                                                                                                                    cTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCa             <  1:3258122/150‑1 (MQ=255)
                                                                                                                                    cTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCa             <  2:2489305/150‑1 (MQ=255)
                                                                                                                                      tAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAc           <  1:1928306/150‑1 (MQ=255)
                                                                                                                                            tGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAc     <  2:924800/150‑1 (MQ=255)
                                                                                                                                             ggATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACg    >  1:764064/1‑150 (MQ=255)
                                                                                                                                              gATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGc   >  1:4725085/1‑150 (MQ=255)
                                                                                                                                               aTGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCa  <  1:2338247/150‑1 (MQ=255)
                                                                                                                                               aTGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCa  <  1:108688/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCA  >  CP009273/79976‑80268

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: