Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,931,363 G→T 9.2% R152L (CGT→CTT)  yebK → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,931,3630GT 9.2% 858.0 / 49.6 239R152L (CGT→CTT) yebKputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (130/87);  new base T (12/10);  total (142/97)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 245 total aligned reads displayed.

AATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAGCTGGCACGCGAAAACG  >  CP009273/1931214‑1931512
                                                                                                                                                      |                                                                                                                                                     
aaTGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCg                                                                                                                                                       >  1:2205700/1‑150 (MQ=255)
    gCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGtttt                                                                                                                                                   <  2:1282632/150‑1 (MQ=255)
     cAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTa                                                                                                                                                  >  1:4172852/1‑150 (MQ=255)
      aaCGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTaa                                                                                                                                                 >  1:182095/1‑150 (MQ=255)
       aCGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAAt                                                                                                                                                <  2:2896960/150‑1 (MQ=255)
         gCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGt                                                                                                                                              <  2:4442828/150‑1 (MQ=255)
         gCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGt                                                                                                                                              <  2:227933/150‑1 (MQ=255)
          cTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGtt                                                                                                                                             <  1:3075312/150‑1 (MQ=255)
          cTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGtt                                                                                                                                             >  1:1219122/1‑150 (MQ=255)
            tGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTcc                                                                                                                                           <  1:889556/150‑1 (MQ=255)
             gATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTTATGTTCCg                                                                                                                                          >  2:2567471/1‑150 (MQ=255)
             gATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCg                                                                                                                                          <  2:2240222/150‑1 (MQ=255)
              aTCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCgg                                                                                                                                         >  1:734848/1‑150 (MQ=255)
               tCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCggt                                                                                                                                        <  2:3800498/150‑1 (MQ=255)
               tCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCggt                                                                                                                                        <  2:3799553/150‑1 (MQ=255)
                cATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCATTTTTCGGATTAGTCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCggtg                                                                                                                                       >  1:4688048/1‑150 (MQ=255)
                cATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATAGCCTTTTTAGGATTTGGCTCTTCAGCCGCCGTTGCCCAAGATGCGATGAATAATTTATTTCGTTTTTATTTTGCggtg                                                                                                                                       >  1:242986/1‑150 (MQ=255)
                   gTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGgtggtg                                                                                                                                    <  1:3309771/150‑1 (MQ=255)
                     acGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTa                                                                                                                                  >  2:2000255/2‑150 (MQ=255)
                         cATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTcc                                                                                                                              >  2:4687221/1‑150 (MQ=255)
                         cATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTcc                                                                                                                              >  2:239230/1‑150 (MQ=255)
                          aTTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCg                                                                                                                             >  2:4663890/1‑150 (MQ=255)
                          aTTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCg                                                                                                                             >  2:4388616/1‑150 (MQ=255)
                           ttCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCga                                                                                                                            >  1:557819/1‑150 (MQ=255)
                            tCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgat                                                                                                                           >  2:1658867/1‑150 (MQ=255)
                            tCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgat                                                                                                                           >  2:1101479/1‑150 (MQ=255)
                             cACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCgatg                                                                                                                          <  1:2673882/150‑1 (MQ=255)
                              aCTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgatga                                                                                                                         <  2:4473691/150‑1 (MQ=255)
                               cTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgatga                                                                                                                         >  1:915611/1‑150 (MQ=255)
                               cTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCgatgat                                                                                                                        >  1:3760627/1‑150 (MQ=255)
                               cTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCgatgat                                                                                                                        >  2:3462580/1‑150 (MQ=255)
                                 ggATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGatat                                                                                                                      <  1:2106624/150‑1 (MQ=255)
                                  gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATc                                                                                                                     <  1:1945284/150‑1 (MQ=255)
                                  gATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATc                                                                                                                     >  1:1360537/1‑150 (MQ=255)
                                   aTAAATCTGCCATCAACCGCGCCGTCGACTTGTTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCGATGATATCg                                                                                                                    >  2:3668809/1‑150 (MQ=255)
                                   aTAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCg                                                                                                                    <  1:2700269/150‑1 (MQ=255)
                                   aTAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCg                                                                                                                    >  2:3317831/1‑150 (MQ=255)
                                   aTAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCg                                                                                                                    <  1:4039646/150‑1 (MQ=255)
                                    tAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGt                                                                                                                   >  1:1005963/1‑150 (MQ=255)
                                     aaaTCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGtg                                                                                                                  >  1:3164053/1‑150 (MQ=255)
                                     aaaTCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGtg                                                                                                                  >  1:795039/1‑150 (MQ=255)
                                       aTCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGtgct                                                                                                                >  1:3705544/1‑150 (MQ=255)
                                          tGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCa                                                                                                             <  2:176322/150‑1 (MQ=255)
                                          tGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCa                                                                                                             >  2:829368/1‑150 (MQ=255)
                                          tGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCa                                                                                                             >  2:1876179/1‑150 (MQ=255)
                                            ccATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAAc                                                                                                           <  2:267997/150‑1 (MQ=255)
                                            ccATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAAc                                                                                                           <  2:617867/150‑1 (MQ=255)
                                            ccATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGACTAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACg                                                                                                          >  2:4362452/1‑150 (MQ=255)
                                                cAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCAt                                                                                                       >  2:1204793/1‑150 (MQ=255)
                                                  aCCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGa                                                                                                     >  2:9284/1‑150 (MQ=255)
                                                       gcCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGt                                                                                                >  2:3216084/1‑150 (MQ=255)
                                                         cGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTAt                                                                                              >  1:1696777/1‑150 (MQ=255)
                                                         cGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTAt                                                                                              >  1:3452787/1‑150 (MQ=255)
                                                          gTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATTAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATg                                                                                             >  2:1410896/1‑150 (MQ=255)
                                                            cGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGaa                                                                                           >  2:2956812/1‑150 (MQ=255)
                                                             gACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCAGGAGTTGTATGAAt                                                                                          >  2:2855454/1‑150 (MQ=255)
                                                              aCTTGCTCACTCAGGCAAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAAt                                                                                          >  1:368887/1‑150 (MQ=255)
                                                              aCTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAAtt                                                                                         >  1:2542081/1‑150 (MQ=255)
                                                                ttGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGt                                                                                       >  1:2879456/1‑150 (MQ=255)
                                                                  gCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAg                                                                                     <  2:3017500/150‑1 (MQ=255)
                                                                  gCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAg                                                                                     <  1:604805/150‑1 (MQ=255)
                                                                   ctcactcaGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGc                                                                                    >  2:1124994/1‑150 (MQ=255)
                                                                    tcactcaGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCg                                                                                   >  1:1323954/1‑150 (MQ=255)
                                                                            gCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGAc                                                                           <  1:4488548/150‑1 (MQ=255)
                                                                             ggcaaaaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGAc                                                                           >  2:2725664/4‑150 (MQ=255)
                                                                                aaaaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACgtg                                                                        <  2:1322114/150‑1 (MQ=255)
                                                                                aaaaaaaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACgtg                                                                        >  1:2206600/1‑150 (MQ=255)
                                                                                     aaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACgtggtggt                                                                   <  2:4754244/150‑1 (MQ=255)
                                                                                     aaTCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACgtggtggt                                                                   <  2:4754233/150‑1 (MQ=255)
                                                                                       tCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGc                                                                 >  1:1112375/1‑150 (MQ=255)
                                                                                        cGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCt                                                                <  2:971104/150‑1 (MQ=255)
                                                                                         gCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTg                                                               >  1:1750565/1‑150 (MQ=255)
                                                                                          ccTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGa                                                              >  2:221976/1‑150 (MQ=255)
                                                                                                 cGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTca                                                       >  2:4090904/1‑150 (MQ=255)
                                                                                                  ggATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTcac                                                      <  2:4660430/150‑1 (MQ=255)
                                                                                                    aTTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCacac                                                    >  2:2125021/1‑150 (MQ=255)
                                                                                                     ttAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACt                                                   >  1:1343471/1‑150 (MQ=255)
                                                                                                      tAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTg                                                  >  1:208062/1‑150 (MQ=255)
                                                                                                            ctTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGgaagaa                                            >  1:2464183/1‑150 (MQ=255)
                                                                                                               cAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACaa                                         <  2:2615158/150‑1 (MQ=255)
                                                                                                                aGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACaaa                                        >  2:335157/1‑150 (MQ=255)
                                                                                                                 gccgccGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACaaaa                                       >  1:21629/1‑150 (MQ=255)
                                                                                                                  ccgccgTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACaaaaa                                      <  2:27858/150‑1 (MQ=255)
                                                                                                                   cgccgTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACaaaaaa                                     >  2:2790911/1‑150 (MQ=255)
                                                                                                                      cgTTGCCCACGATGCGATGAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCt                                  >  1:355217/1‑150 (MQ=255)
                                                                                                                       gTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTg                                 <  2:4217141/150‑1 (MQ=255)
                                                                                                                         tGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGt                               <  2:3550971/150‑1 (MQ=255)
                                                                                                                         tGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGt                               <  2:3408626/150‑1 (MQ=255)
                                                                                                                         tGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGt                               <  1:1483642/150‑1 (MQ=255)
                                                                                                                           cccACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCg                             >  1:1121546/1‑150 (MQ=255)
                                                                                                                                gatgcgatgAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTg                        >  2:2357344/1‑150 (MQ=255)
                                                                                                                                  tgcgatgAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGgc                      <  1:3791688/150‑1 (MQ=255)
                                                                                                                                    cgatgAATAAGTTCTTTCTTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGgcgc                    >  1:153983/1‑150 (MQ=255)
                                                                                                                                    cgatgAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGgcgc                    >  1:3033219/1‑150 (MQ=255)
                                                                                                                                     gatgAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCa                   >  1:1766952/1‑150 (MQ=255)
                                                                                                                                      atgAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAg                  >  1:2050322/1‑150 (MQ=255)
                                                                                                                                            aaGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAGCTGGCa            >  2:55328/1‑150 (MQ=255)
                                                                                                                                              gTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAGCTGGCAcg          <  1:4456364/150‑1 (MQ=255)
                                                                                                                                                    tCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAGCTGGCACGCGaaaa    <  2:536341/150‑1 (MQ=255)
                                                                                                                                                      gTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAGCTGGCACGCGAAAacg  >  2:2842341/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAAAATCGCCTTTTTCGGATTAGGCTCTTCAGCCGCCGTTGCCCACGATGCGATGAATAAGTTCTTTCGTTTTAATGTTCCGGTGGTGTACTCCGATGATATCGTGCTGCAACGCATGAGTTGTATGAATTGTAGCGACGGAGACGTGGTGGTGCTGATTTCTCACACTGGAAGAACAAAAAATCTGGTCGAGCTGGCGCAGCTGGCACGCGAAAACG  >  CP009273/1931214‑1931512

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: