Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,526,794 C→T 100% intergenic (+55/‑329) cysK → / → ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,526,7940CT98.2% 722.3 / 2.6 219intergenic (+55/‑329)cysK/ptsHcysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (117/98);  total (119/100)
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.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 219 total aligned reads displayed.

TATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCC  >  CP009273/2526649‑2526938
                                                                                                                                                 |                                                                                                                                                
 atTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTg                                                                                                                                               >  1:1578809/1‑148 (MQ=255)
    gtgGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCtt                                                                                                                                            <  2:656384/148‑1 (MQ=255)
      ggTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCa                                                                                                                                          >  2:2391353/1‑148 (MQ=255)
        ttATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCaaa                                                                                                                                        >  1:4020/1‑148 (MQ=255)
         tATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAAc                                                                                                                                       <  2:2965561/148‑1 (MQ=255)
               aCCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGc                                                                                                                                 >  2:502649/1‑148 (MQ=255)
                ccATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGcc                                                                                                                                >  1:2400385/1‑148 (MQ=255)
                    catcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCcctc                                                                                                                            <  2:2872763/148‑1 (MQ=255)
                         ggTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGc                                                                                                                       >  1:2794044/1‑148 (MQ=255)
                         ggTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGc                                                                                                                       >  1:2011132/1‑148 (MQ=255)
                             aGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCCTTg                                                                                                                   >  2:2659773/1‑148 (MQ=255)
                             aGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTg                                                                                                                   >  2:54044/1‑148 (MQ=255)
                             aGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTg                                                                                                                   >  1:3120324/1‑148 (MQ=255)
                              gCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAACTGCGTAAACCAGCCCCACCTTAGCTGCAAAATTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGa                                                                                                                  <  2:3347979/148‑1 (MQ=255)
                              gCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGa                                                                                                                  >  1:631979/1‑148 (MQ=255)
                               cGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGAt                                                                                                                 >  2:1147135/1‑148 (MQ=255)
                                 ttATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTc                                                                                                               >  2:405951/1‑148 (MQ=255)
                                  tATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCa                                                                                                              >  2:391630/1‑148 (MQ=255)
                                    tttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGc                                                                                                            <  2:1114375/148‑1 (MQ=255)
                                    tttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGc                                                                                                            >  2:1325208/1‑148 (MQ=255)
                                     ttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGcc                                                                                                           <  2:1850696/148‑1 (MQ=255)
                                     ttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGcc                                                                                                           <  2:1165841/148‑1 (MQ=255)
                                        aGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGt                                                                                                        <  1:1478358/148‑1 (MQ=255)
                                            ccGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGa                                                                                                    >  2:149323/1‑148 (MQ=255)
                                             cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGaa                                                                                                   <  1:745833/148‑1 (MQ=255)
                                             cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGaa                                                                                                   >  1:1326269/1‑148 (MQ=255)
                                             cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGaa                                                                                                   >  2:2086384/1‑148 (MQ=255)
                                                  tGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGt                                                                                              <  1:412478/148‑1 (MQ=255)
                                                  tGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGt                                                                                              >  2:3228227/1‑148 (MQ=255)
                                                   gTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                                             >  1:558908/1‑148 (MQ=255)
                                                   gTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTa                                                                                             <  1:1169843/148‑1 (MQ=255)
                                                     ttGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTAtt                                                                                           <  2:3060828/148‑1 (MQ=255)
                                                        ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTaa                                                                                        <  1:1318204/148‑1 (MQ=255)
                                                        ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTaa                                                                                        <  2:340094/148‑1 (MQ=255)
                                                          gATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAAcc                                                                                      <  1:2708216/148‑1 (MQ=255)
                                                            tctcTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAg                                                                                    >  1:3127599/1‑148 (MQ=255)
                                                              tctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGAc                                                                                  <  1:713637/148‑1 (MQ=255)
                                                               ctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACt                                                                                 >  2:72229/1‑148 (MQ=255)
                                                               ctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACt                                                                                 >  2:72292/1‑148 (MQ=255)
                                                                ttCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTa                                                                                >  2:1441863/1‑148 (MQ=255)
                                                                 tCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTaa                                                                               <  2:3484416/148‑1 (MQ=255)
                                                                    cTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAatta                                                                            >  2:871419/1‑148 (MQ=255)
                                                                    cTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTGTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAatta                                                                            >  1:3159268/1‑148 (MQ=255)
                                                                     tGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAattat                                                                           <  1:3077539/148‑1 (MQ=255)
                                                                      gagaAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAattatt                                                                          <  2:1408597/148‑1 (MQ=255)
                                                                       agaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAttt                                                                         >  2:1024714/1‑148 (MQ=255)
                                                                       agaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTCTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAttt                                                                         >  2:2777819/1‑148 (MQ=255)
                                                                        gaaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTAtttt                                                                        >  1:1234440/1‑148 (MQ=255)
                                                                         aaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTg                                                                       <  1:84841/148‑1 (MQ=255)
                                                                         aaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTg                                                                       <  1:899967/148‑1 (MQ=255)
                                                                         aaagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTg                                                                       >  2:230640/1‑148 (MQ=255)
                                                                          aagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTCTTTTGa                                                                      >  2:1731892/1‑148 (MQ=255)
                                                                          aagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGa                                                                      <  1:1617276/148‑1 (MQ=255)
                                                                            gaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATg                                                                    <  2:2864744/148‑1 (MQ=255)
                                                                              aTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATgcg                                                                  <  1:2394333/148‑1 (MQ=255)
                                                                                tGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGcgcg                                                                >  2:2451689/1‑148 (MQ=255)
                                                                                  cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGaa                                                              >  1:1164829/1‑148 (MQ=255)
                                                                                   aaCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGaaa                                                             <  2:1258165/148‑1 (MQ=255)
                                                                                    aCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAAt                                                            <  2:2171949/148‑1 (MQ=255)
                                                                                       gTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTaa                                                         >  1:2921458/1‑148 (MQ=255)
                                                                                        tAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAAt                                                        <  1:2633659/148‑1 (MQ=255)
                                                                                            gCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGtt                                                    >  1:3197447/1‑148 (MQ=255)
                                                                                              cAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTAc                                                  >  2:2218265/1‑148 (MQ=255)
                                                                                                gCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAg                                                >  1:3474292/1‑148 (MQ=255)
                                                                                                     ttAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaaa                                           <  2:2618279/148‑1 (MQ=255)
                                                                                                     ttAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaaa                                           >  2:2446646/1‑148 (MQ=255)
                                                                                                     ttAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGaaaa                                           >  2:685640/1‑148 (MQ=255)
                                                                                                           aTGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAGGCCaaa                                     >  2:597971/1‑148 (MQ=255)
                                                                                                                tAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGa                                <  2:428539/148‑1 (MQ=255)
                                                                                                                 aaaaaaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGaa                               >  1:1056514/1‑148 (MQ=255)
                                                                                                                  aaaaaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt                              <  2:3267037/148‑1 (MQ=255)
                                                                                                                  aaaaaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAAt                              <  1:2788806/148‑1 (MQ=255)
                                                                                                                   aaaaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATc                             >  2:1252536/1‑148 (MQ=255)
                                                                                                                    aaaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCg                            >  1:2852580/1‑148 (MQ=255)
                                                                                                                    aaaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCg                            >  2:1594897/1‑148 (MQ=255)
                                                                                                                     aaGCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGa                           >  1:539972/1‑148 (MQ=255)
                                                                                                                       gCACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAtt                         >  2:3140147/1‑148 (MQ=255)
                                                                                                                        cACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAttt                        >  1:961753/1‑148 (MQ=255)
                                                                                                                        cACCTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGAttt                        <  2:1745129/148‑1 (MQ=255)
                                                                                                                          ccTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTa                      <  1:391630/148‑1 (MQ=255)
                                                                                                                          ccTTTTTAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTAGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTa                      >  2:1939655/1‑148 (MQ=255)
                                                                                                                            tttttAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATg                    <  1:1383300/148‑1 (MQ=255)
                                                                                                                            tttttAGGTGCTTTTGTGTGGTCTTATTCAAACTTTCGCGCCTCCTGGCATTGATTCAGCCTGTCGTAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTATTCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTAtt                    <  2:2487617/148‑2 (MQ=255)
                                                                                                                             ttttAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGa                   >  2:1543286/1‑148 (MQ=255)
                                                                                                                              tttAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAt                  <  2:2787188/148‑1 (MQ=255)
                                                                                                                               ttAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAtt                 <  1:2451081/148‑1 (MQ=255)
                                                                                                                                tttgTGCTTTTTTGTGGTCTGATTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTAGGAACTGGTATTTAACAAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAttt                <  1:697585/145‑1 (MQ=255)
                                                                                                                                tAGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGAttt                <  1:186367/148‑1 (MQ=255)
                                                                                                                                 aGGTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTg               <  1:1913463/148‑1 (MQ=255)
                                                                                                                                   gTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGt             >  2:2012090/1‑148 (MQ=255)
                                                                                                                                   gTGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGt             >  1:2883782/1‑148 (MQ=255)
                                                                                                                                    tGCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGtt            >  1:2556070/1‑148 (MQ=255)
                                                                                                                                     gCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTc           >  2:36643/1‑148 (MQ=255)
                                                                                                                                     gCTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTc           >  1:1500326/1‑148 (MQ=255)
                                                                                                                                      cTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCa          >  1:1029983/1‑148 (MQ=255)
                                                                                                                                      cTTTTTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCa          >  1:1030044/1‑148 (MQ=255)
                                                                                                                                        tttttGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAt        >  2:527209/1‑148 (MQ=255)
                                                                                                                                        tttttGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAt        <  2:3117905/148‑1 (MQ=255)
                                                                                                                                           ttGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAAttct     <  1:502649/148‑1 (MQ=255)
                                                                                                                                              tgGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTcc  >  2:190939/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TATTGTGGTTATTCTACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTAATTATTTTGATGCGCGAAATTAATCGTTACAGGAAAAGCCAAAGCTGAATCGATTTTATGATTTGGTTCAATTCTTCC  >  CP009273/2526649‑2526938

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: