Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 80,123 C→T 100% 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,1230CT96.4% 757.0 / 11.1 219G25S (GGC→AGC) leuLleu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base T (117/94);  total (119/100)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 220 total aligned reads displayed.

TCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGACGCA  >  CP009273/79976‑80268
                                                                                                                                                   |                                                                                                                                                 
  gcgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTc                                                                                                                                               <  2:3995026/150‑1 (MQ=255)
  gcgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTc                                                                                                                                               <  2:2625318/150‑1 (MQ=255)
   cgcTATGTGGTATCGAAAATAATGCCTTGCTGGCTCATGGTGTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCa                                                                                                                                              <  1:4218698/150‑1 (MQ=255)
   cgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCa                                                                                                                                              <  1:2568529/150‑1 (MQ=255)
    gcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCAc                                                                                                                                             >  2:3746052/1‑150 (MQ=255)
     cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACt                                                                                                                                            <  1:2510912/150‑1 (MQ=255)
      aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAATCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTc                                                                                                                                           >  2:780759/1‑150 (MQ=255)
           ggTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTa                                                                                                                                      >  1:2018803/1‑150 (MQ=255)
              aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcg                                                                                                                                   >  2:3841882/1‑150 (MQ=255)
              aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcg                                                                                                                                   >  1:4265062/1‑150 (MQ=255)
              aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcg                                                                                                                                   >  1:1773152/1‑150 (MQ=255)
                cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCACTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcgcg                                                                                                                                 <  2:3665734/150‑1 (MQ=255)
                cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcgcg                                                                                                                                 >  2:4599564/1‑150 (MQ=255)
                cGAAAATAATGACTTGCTGGCTCAGGGTTTGGGTCCTTGTCTCTTGTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACcgcg                                                                                                                                 >  1:2073264/1‑150 (MQ=255)
                 gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCgcgc                                                                                                                                >  2:2113928/1‑150 (MQ=255)
                  aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCa                                                                                                                               >  2:4215978/1‑150 (MQ=255)
                     ataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAaag                                                                                                                            >  2:2424511/1‑150 (MQ=255)
                     ataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAaag                                                                                                                            >  1:1339950/1‑150 (MQ=255)
                      taatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAaaga                                                                                                                           <  1:3269638/150‑1 (MQ=255)
                        atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAaagaag                                                                                                                         <  1:371732/150‑1 (MQ=255)
                             ttGCTGGCTCATTGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGTGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTTACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCg                                                                                                                    <  2:731388/150‑1 (MQ=255)
                             ttGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCg                                                                                                                    >  2:1522885/1‑150 (MQ=255)
                                cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                                                                 >  1:4472766/1‑150 (MQ=255)
                                 tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                                                                                                                <  1:2664067/150‑1 (MQ=255)
                                    cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                                                                             >  1:3104550/1‑150 (MQ=255)
                                     tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag                                                                                                            >  2:4263900/1‑150 (MQ=255)
                                      cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagt                                                                                                           <  2:3632068/150‑1 (MQ=255)
                                        tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtag                                                                                                         <  2:662386/150‑1 (MQ=255)
                                        tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtag                                                                                                         >  1:420674/1‑150 (MQ=255)
                                            ttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGa                                                                                                     <  1:182480/150‑1 (MQ=255)
                                             tGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                                                                    <  2:110365/150‑1 (MQ=255)
                                                gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCACTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                                                                 <  1:935696/150‑1 (MQ=255)
                                                 tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGAt                                                                                                >  1:2313626/1‑150 (MQ=255)
                                                 tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGAt                                                                                                <  2:1773152/150‑1 (MQ=255)
                                                 tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGAt                                                                                                <  2:4265062/150‑1 (MQ=255)
                                                   cTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaa                                                                                              <  1:441592/150‑1 (MQ=255)
                                                     tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                                                                            >  1:1615490/1‑150 (MQ=255)
                                                      gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                                                                           <  2:3079839/150‑1 (MQ=255)
                                                        ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGa                                                                                         <  1:1464877/150‑1 (MQ=255)
                                                          ctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                                                                       <  1:2496132/150‑1 (MQ=255)
                                                          ctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                                                                       <  1:2504579/150‑1 (MQ=255)
                                                                   cgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGt                                                                              >  2:2418498/1‑150 (MQ=255)
                                                                     ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                                                                            >  2:2993090/1‑150 (MQ=255)
                                                                       tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAtt                                                                          >  1:2398663/1‑150 (MQ=255)
                                                                        cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                                                                         >  2:1396849/1‑150 (MQ=255)
                                                                        cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                                                                         >  1:482250/1‑150 (MQ=255)
                                                                        cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTa                                                                         >  1:2126123/1‑150 (MQ=255)
                                                                        cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGAGGGTTTATTGTTTGACTGCGTGCAGGCTTAATGCTGGATGCTGCTCACACGACTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGACCGATGTGAGTCATTa                                                                         >  1:166035/1‑150 (MQ=255)
                                                                         gCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTaa                                                                        >  1:3926617/1‑150 (MQ=255)
                                                                         gCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTaa                                                                        >  1:3763115/1‑150 (MQ=255)
                                                                            tCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAACGCGAACGATGTGAGTCATTAAATc                                                                     >  2:3561407/1‑150 (MQ=255)
                                                                            tCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATc                                                                     >  1:3364001/1‑150 (MQ=255)
                                                                            tCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATc                                                                     >  1:2954596/1‑150 (MQ=255)
                                                                             cGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCa                                                                    >  2:3620080/1‑150 (MQ=255)
                                                                                gCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCt                                                                 <  1:591414/150‑1 (MQ=255)
                                                                                gCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCt                                                                 <  2:3129863/150‑1 (MQ=255)
                                                                                   tAAAAATACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGAGTGCTGGCTTAATTCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATTTGAGTCATTAAATCAGCTCCa                                                              >  2:2023744/1‑150 (MQ=255)
                                                                                   tAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa                                                              >  2:2151858/1‑150 (MQ=255)
                                                                                    aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg                                                             >  1:740865/1‑150 (MQ=255)
                                                                                    aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg                                                             >  2:1246412/1‑150 (MQ=255)
                                                                                         aaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa                                                        >  2:453834/1‑150 (MQ=255)
                                                                                         aaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa                                                        >  1:536743/1‑150 (MQ=255)
                                                                                         aaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa                                                        >  1:1558817/1‑150 (MQ=255)
                                                                                             cgcgCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCg                                                    >  1:2670066/1‑150 (MQ=255)
                                                                                             cgcgCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCg                                                    <  1:4620989/150‑1 (MQ=255)
                                                                                               cgcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGat                                                  >  1:298830/1‑150 (MQ=255)
                                                                                                 cAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGatat                                                >  2:303993/1‑150 (MQ=255)
                                                                                                   aTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGc                                              <  1:2168318/150‑1 (MQ=255)
                                                                                                   aTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTT‑‑‑TAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCttt                                           >  1:2619286/1‑150 (MQ=255)
                                                                                                    tGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCt                                             <  1:463649/150‑1 (MQ=255)
                                                                                                     ggCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCtt                                            >  1:4088545/1‑150 (MQ=255)
                                                                                                     ggCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCtt                                            >  2:4426512/1‑150 (MQ=255)
                                                                                                          gggTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTaga                                       <  2:4408003/150‑1 (MQ=255)
                                                                                                          gggTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTaga                                       <  1:4587680/150‑1 (MQ=255)
                                                                                                                   gtttgACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGa                              >  1:3682819/1‑150 (MQ=255)
                                                                                                                      tgACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATAc                           >  2:536336/1‑150 (MQ=255)
                                                                                                                      tgACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATAc                           >  2:1143566/1‑150 (MQ=255)
                                                                                                                       gACTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACa                          >  1:1721355/1‑150 (MQ=255)
                                                                                                                         cTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaa                        >  1:2578478/1‑150 (MQ=255)
                                                                                                                         cTGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaa                        <  1:3468332/150‑1 (MQ=255)
                                                                                                                         cTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaa                        <  1:126046/150‑1 (MQ=255)
                                                                                                                          tGCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaaa                       >  1:1148639/1‑150 (MQ=255)
                                                                                                                           gCGTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACaaaaa                      <  1:419708/150‑1 (MQ=255)
                                                                                                                             gTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACg                    >  2:1295601/1‑150 (MQ=255)
                                                                                                                             gTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACg                    <  1:98319/150‑1 (MQ=255)
                                                                                                                             gTGCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACg                    <  2:3304963/150‑1 (MQ=255)
                                                                                                                               gCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGa                  >  1:415472/1‑150 (MQ=255)
                                                                                                                               gCTGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGa                  <  1:3797683/150‑1 (MQ=255)
                                                                                                                                 tGGCTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTATATCAGCTCCAGATGAATGCGATATGATTTTAGAGTTACTGGATACAAAAACGGATg                <  2:3287950/150‑1 (MQ=255)
                                                                                                                                    cTTAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCa             <  2:351780/150‑1 (MQ=255)
                                                                                                                                      tAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAc           <  2:2779929/150‑1 (MQ=255)
                                                                                                                                      tAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAc           >  1:2339058/1‑150 (MQ=255)
                                                                                                                                      tAATGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAc           <  2:1872930/150‑1 (MQ=255)
                                                                                                                                       aaTGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAcc          >  1:2908891/1‑150 (MQ=255)
                                                                                                                                       aaTGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAcc          >  1:2887388/1‑150 (MQ=255)
                                                                                                                                       aaTGCTGGATGCTGCTCACTCGTCTACAGCGATAAGAATTTGCGTTTAGTAGTTGTAGACCGATTTAGCGATCGATGTTAGTAATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGTTACATAAACGGTTGTCAAcc          <  2:1150089/150‑1 (MQ=255)
                                                                                                                                        aTGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAAccc         >  2:3238024/1‑150 (MQ=255)
                                                                                                                                           cTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGACTGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGa      >  2:2939703/1‑150 (MQ=255)
                                                                                                                                           cTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGa      >  2:4033112/1‑150 (MQ=255)
                                                                                                                                            tGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATGCTTTTAGAGTTACTGGATACAAAAACGGATGTCAACCCTGAc     <  2:2444416/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: