Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,501,358 G→T 7.4% R6S (CGC→AGC)  fryC ← putative enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,501,3580GT 7.4% 1007.1 / 35.2 270R6S (CGC→AGC) fryCputative enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base G (101/149);  new base T (8/12);  total (109/161)
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.98e-01

Only 100 of 278 total aligned reads displayed.

CACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGT  >  CP009273/2501209‑2501507
                                                                                                                                                     |                                                                                                                                                     
cacaCCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCg                                                                                                                                                       >  2:838289/1‑150 (MQ=255)
cacaCCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCg                                                                                                                                                       <  1:2415881/150‑1 (MQ=255)
cacaCCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCg                                                                                                                                                       <  1:2415731/150‑1 (MQ=255)
    ccGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTCCGGGAGTTCACCCCAGAAGCAGGTTTTAGAGGCCTCCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttt                                                                                                                                                   >  1:2931308/1‑150 (MQ=255)
      gCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttttt                                                                                                                                                 >  2:1674054/1‑150 (MQ=255)
      gCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttttt                                                                                                                                                 >  2:32691/1‑150 (MQ=255)
      gCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGtttttt                                                                                                                                                 >  2:33011/1‑150 (MQ=255)
       cccATAATTAAGGTCGGCGCCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTa                                                                                                                                                >  1:3337859/1‑150 (MQ=255)
        ccATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTaa                                                                                                                                               <  2:2636301/150‑1 (MQ=255)
                aaGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATa                                                                                                                                       <  1:3584467/150‑1 (MQ=255)
                aaGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATa                                                                                                                                       <  1:1313575/150‑1 (MQ=255)
                 aGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATaa                                                                                                                                      >  2:168299/1‑150 (MQ=255)
                 aGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATaa                                                                                                                                      <  1:657738/150‑1 (MQ=255)
                 aGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATaa                                                                                                                                      <  2:2603808/150‑1 (MQ=255)
                  ggTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAAc                                                                                                                                     >  1:1860552/1‑150 (MQ=255)
                  ggTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAAc                                                                                                                                     >  1:1860448/1‑150 (MQ=255)
                      ggCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAAtt                                                                                                                                 <  2:3361288/150‑1 (MQ=255)
                          ccATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATcc                                                                                                                             <  1:172501/150‑1 (MQ=255)
                            aTGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCt                                                                                                                           >  1:2252258/1‑150 (MQ=255)
                            aTGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCt                                                                                                                           >  1:3826882/1‑150 (MQ=255)
                              gCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGt                                                                                                                         <  2:1392583/150‑1 (MQ=255)
                               cGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTc                                                                                                                        <  2:4231449/150‑1 (MQ=255)
                                  gAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGaa                                                                                                                     <  2:2871210/150‑1 (MQ=255)
                                   aaaTCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAat                                                                                                                    <  1:3582850/150‑1 (MQ=255)
                                    aaTCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAata                                                                                                                   <  1:4386342/150‑1 (MQ=255)
                                      tCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACTACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATAtt                                                                                                                 <  1:541811/150‑1 (MQ=255)
                                       cccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                                >  1:390370/1‑150 (MQ=255)
                                       cccTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTa                                                                                                                <  2:3843679/150‑1 (MQ=255)
                                          tGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACATCGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtg                                                                                                             <  2:1024914/150‑1 (MQ=255)
                                          tGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtg                                                                                                             >  1:2931808/1‑150 (MQ=255)
                                            aCATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgtt                                                                                                           >  1:1977503/1‑150 (MQ=255)
                                             cATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACtgttg                                                                                                          <  1:896608/150‑1 (MQ=255)
                                               tCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTgttgtt                                                                                                        >  2:1560720/1‑150 (MQ=255)
                                                 acaTGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCa                                                                                                      >  2:3555913/1‑150 (MQ=255)
                                                  caTGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCag                                                                                                     >  2:1812423/1‑150 (MQ=255)
                                                       tgcGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCa                                                                                                <  2:2722686/150‑1 (MQ=255)
                                                        gcGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCaa                                                                                               <  1:610168/150‑1 (MQ=255)
                                                              gTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCAtct                                                                                         <  1:3861728/150‑1 (MQ=255)
                                                               ttCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCAtctc                                                                                        <  2:419719/150‑1 (MQ=255)
                                                               ttCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCAtctc                                                                                        <  2:419028/150‑1 (MQ=255)
                                                                  aCCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTc                                                                                     <  1:1844841/150‑1 (MQ=255)
                                                                   ccccAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCg                                                                                    <  1:732788/150‑1 (MQ=255)
                                                                   ccccAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCg                                                                                    <  1:1547266/150‑1 (MQ=255)
                                                                        gAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTc                                                                               >  1:3800117/1‑150 (MQ=255)
                                                                        gAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTc                                                                               <  1:3376186/150‑1 (MQ=255)
                                                                        gAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTc                                                                               >  1:964639/1‑150 (MQ=255)
                                                                        gAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTc                                                                               <  2:1292901/150‑1 (MQ=255)
                                                                           gCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCttt                                                                            >  2:773860/1‑150 (MQ=255)
                                                                            cTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTtg                                                                           >  2:3455148/1‑150 (MQ=255)
                                                                                      agGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCg                                                                 >  2:1204836/1‑150 (MQ=255)
                                                                                         ccTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAg                                                              <  1:4307248/150‑1 (MQ=255)
                                                                                            gCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGt                                                           >  2:4322773/1‑150 (MQ=255)
                                                                                            gCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGt                                                           >  2:419959/1‑150 (MQ=255)
                                                                                              ggATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGttt                                                         <  2:3826085/150‑1 (MQ=255)
                                                                                               gATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGtttt                                                        <  1:766959/150‑1 (MQ=255)
                                                                                                aTTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGttttt                                                       <  1:1096907/150‑1 (MQ=255)
                                                                                                 ttCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTa                                                      <  2:2086096/150‑1 (MQ=255)
                                                                                                 ttCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTa                                                      <  2:2085384/150‑1 (MQ=255)
                                                                                                 ttCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTa                                                      >  1:812166/1‑150 (MQ=255)
                                                                                                 ttCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTa                                                      <  2:1737477/150‑1 (MQ=255)
                                                                                                    ttAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAAtt                                                   <  2:1026788/150‑1 (MQ=255)
                                                                                                      aaCTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGc                                                 >  2:872174/1‑150 (MQ=255)
                                                                                                      aaCTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGc                                                 >  1:331265/1‑150 (MQ=255)
                                                                                                       aCTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCg                                                <  2:1625317/150‑1 (MQ=255)
                                                                                                        ctgccgctgccgCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGt                                               <  2:2138729/150‑1 (MQ=255)
                                                                                                            cgctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTg                                           <  2:3762604/150‑1 (MQ=255)
                                                                                                             gctgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGc                                          <  1:3785448/150‑1 (MQ=255)
                                                                                                               tgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTCATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTATTTGCGGCCTGcga                                        >  2:4452647/1‑148 (MQ=255)
                                                                                                               tgccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCaa                                        >  1:1526906/1‑150 (MQ=255)
                                                                                                                gccgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCaaa                                       <  1:2389325/150‑1 (MQ=255)
                                                                                                                  cgcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCACAGt                                     >  1:3510983/1‑150 (MQ=255)
                                                                                                                   gcACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTg                                    <  2:2026027/150‑1 (MQ=255)
                                                                                                                     aCCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAt                                  <  2:1658722/150‑1 (MQ=255)
                                                                                                                     aCCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAt                                  >  2:1000807/1‑150 (MQ=255)
                                                                                                                       cGGATGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGAttt                                <  2:321393/150‑1 (MQ=255)
                                                                                                                         gATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCa                              <  2:1829791/150‑1 (MQ=255)
                                                                                                                           tGCGCCAGGCACAACGGTTGCACTGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATa                            <  1:1290082/150‑1 (MQ=255)
                                                                                                                           tGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATa                            >  1:1231345/1‑150 (MQ=255)
                                                                                                                           tGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATa                            >  1:118462/1‑150 (MQ=255)
                                                                                                                                cAGGCACAACGTTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCGTTGATGATGCCCGCATCGTTTTTAATTGCGACCTGCAAAGTGATTTCATAAACGt                       >  1:1889780/1‑150 (MQ=255)
                                                                                                                                cAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGt                       <  2:4448414/150‑1 (MQ=255)
                                                                                                                                    cacaACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCGCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTcgcg                   <  1:1024683/150‑1 (MQ=255)
                                                                                                                                     acaACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGt                  >  1:1508572/1‑150 (MQ=255)
                                                                                                                                      caACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTg                 >  2:3030677/1‑150 (MQ=255)
                                                                                                                                      caACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTg                 <  1:4158608/150‑1 (MQ=255)
                                                                                                                                        aCGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAt               >  1:2001452/1‑150 (MQ=255)
                                                                                                                                        aCGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAt               >  1:3109457/1‑150 (MQ=255)
                                                                                                                                        aCGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAt               <  1:446539/150‑1 (MQ=255)
                                                                                                                                         cGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGAtt              <  2:1962625/150‑1 (MQ=255)
                                                                                                                                           gTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCg            >  2:1306210/1‑150 (MQ=255)
                                                                                                                                           gTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCg            <  1:1449008/150‑1 (MQ=255)
                                                                                                                                            ttGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGa           >  2:2281797/1‑150 (MQ=255)
                                                                                                                                              gCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGaaa         <  2:4487195/150‑1 (MQ=255)
                                                                                                                                                aCTGCGTTTTTTAATGGCCATAACATTTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcg       <  2:3022517/150‑1 (MQ=255)
                                                                                                                                                aCTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcg       >  1:2186717/1‑150 (MQ=255)
                                                                                                                                                 cTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcgt      >  1:1603345/1‑150 (MQ=255)
                                                                                                                                                  tGCTTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcgtt     >  1:2285020/1‑150 (MQ=255)
                                                                                                                                                   gCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcgttc    <  2:65827/150‑1 (MQ=255)
                                                                                                                                                     gTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcgttcgt  <  2:4145532/150‑1 (MQ=255)
                                                                                                                                                     gTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAAcgttcgt  <  1:3455158/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACACCGCCCATAATTAAGGTCGGCACCATGCGTGAAATCCCTGACATCACATGCTGCGGGAGTTCACCCCAGAAGCTGGTTTTAGAGGCCTGCGGATTCTTAACTGCCGCTGCCGCACCGGATGCGCCAGGCACAACGGTTGCACTGCGTTTTTTAATGGCCATAACAATTATCCCTGTCGAATATTACTGTTGTTCAGAGGCAATCATCTCTTCGATTTCTTTGATGATGCCCGCAGCGTTTTTAATTGCGTCCTGCAAAGTGATTTCATAAACGTCGCGTGATTCGAAACGTTCGT  >  CP009273/2501209‑2501507

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: