Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,637,458 G→T 20.9% S14Y (TCT→TAT)  quuQ ← Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,637,4580GT20.9% 595.9 / 105.1 229S14Y (TCT→TAT) quuQQin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base G (79/102);  new base T (25/23);  total (104/125)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

AACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATC  >  CP009273/1637311‑1637604
                                                                                                                                                   |                                                                                                                                                  
 aCTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGgaga                                                                                                                                                   >  2:1783000/1‑148 (MQ=255)
 aCTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGgaga                                                                                                                                                   >  2:2758268/1‑148 (MQ=255)
  cTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAt                                                                                                                                                  <  1:78120/148‑1 (MQ=255)
    gAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAttt                                                                                                                                                >  1:3347341/1‑148 (MQ=255)
     aGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCACGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAtttt                                                                                                                                               <  2:490311/148‑1 (MQ=255)
      gggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTg                                                                                                                                              <  1:1896401/148‑1 (MQ=255)
       ggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgg                                                                                                                                             >  2:1257988/1‑148 (MQ=255)
          ccTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                                                                                                          <  2:3046016/148‑1 (MQ=255)
          ccTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                                                                                                          <  1:2334689/148‑1 (MQ=255)
          ccTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAGTTTGGGGa                                                                                                                                          >  1:1403781/1‑148 (MQ=255)
           cTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAg                                                                                                                                         <  2:2879681/148‑1 (MQ=255)
            tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCCATGTTGGAGATTTTGGGGAGt                                                                                                                                        >  2:2126954/1‑148 (MQ=255)
            tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGt                                                                                                                                        >  1:160162/1‑148 (MQ=255)
            tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGt                                                                                                                                        >  2:2126535/1‑148 (MQ=255)
            tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGt                                                                                                                                        >  1:881552/1‑148 (MQ=255)
             gAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTa                                                                                                                                       <  2:763417/148‑1 (MQ=255)
             gAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTa                                                                                                                                       >  2:1983214/1‑148 (MQ=255)
               acacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAat                                                                                                                                     <  1:575589/148‑1 (MQ=255)
               acacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAat                                                                                                                                     <  1:3074629/148‑1 (MQ=255)
                cacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAata                                                                                                                                    <  1:2625670/148‑1 (MQ=255)
                cacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAata                                                                                                                                    >  2:3087990/1‑148 (MQ=255)
                 acCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAatat                                                                                                                                   <  1:412088/148‑1 (MQ=255)
                  ccAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtt                                                                                                                                  <  1:1859616/148‑1 (MQ=255)
                   cAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAttt                                                                                                                                 <  1:479664/148‑1 (MQ=255)
                    aaCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtttt                                                                                                                                <  1:3056465/148‑1 (MQ=255)
                     aCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTg                                                                                                                               >  2:1189687/1‑148 (MQ=255)
                     aCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTg                                                                                                                               <  1:1283807/148‑1 (MQ=255)
                      cTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTgg                                                                                                                              >  1:2543591/1‑148 (MQ=255)
                         ttGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTaa                                                                                                                           <  2:1951730/148‑1 (MQ=255)
                             cAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTTCACAATCCCCTGCGCTGCCTTTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                                                                                                                       <  2:749499/148‑1 (MQ=255)
                             cAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                                                                                                                       <  2:592504/148‑1 (MQ=255)
                             cAGAAATAAGGCCAACCCAAGTGGTGCTTTCGAATGCACCATCCCGTGCGCTGCCATTACATCGGTAATTGTTATACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGTATATTTTGGGGAGTAATATTTTGGTAAGGCt                                                                                                                       <  1:2953776/148‑1 (MQ=255)
                                aaaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCa                                                                                                                    >  1:1020919/1‑148 (MQ=255)
                                aaaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCa                                                                                                                    >  1:2264956/1‑148 (MQ=255)
                                  aTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAg                                                                                                                  >  1:1574627/1‑148 (MQ=255)
                                  aTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAg                                                                                                                  >  1:1198618/1‑148 (MQ=255)
                                  aTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAg                                                                                                                  >  1:131330/1‑148 (MQ=255)
                                         cAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATg                                                                                                           <  1:769861/148‑1 (MQ=255)
                                          aaCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCTCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGc                                                                                                          <  2:2708163/148‑1 (MQ=255)
                                          aaCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATCTTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGc                                                                                                          >  2:3296169/1‑148 (MQ=255)
                                          aaCCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGc                                                                                                          <  1:607455/148‑1 (MQ=255)
                                              cAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGt                                                                                                      <  1:2050748/148‑1 (MQ=255)
                                                aGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTaa                                                                                                                           >  1:3112444/1‑125 (MQ=255)
                                                aGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTaa                                                                                                                           <  2:1361589/125‑1 (MQ=255)
                                                 gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct                                                                                                   <  2:1868722/148‑1 (MQ=255)
                                                 gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct                                                                                                   <  2:459090/148‑1 (MQ=255)
                                                   ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc                                                                                                 >  1:737835/1‑148 (MQ=255)
                                                    gtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCa                                                                                                >  1:3194818/1‑148 (MQ=255)
                                                    gtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCa                                                                                                >  2:1971621/1‑148 (MQ=255)
                                                          tCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACg                                                                                          <  2:3204220/148‑1 (MQ=255)
                                                               tGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGt                                                                                     <  1:3037875/148‑1 (MQ=255)
                                                                  aCCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACg                                                                                  <  2:141292/148‑1 (MQ=255)
                                                                      tCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCaa                                                                              >  1:1467750/1‑148 (MQ=255)
                                                                      tCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCaa                                                                              <  1:688693/148‑1 (MQ=255)
                                                                       ccccTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAAt                                                                             >  1:3409976/1‑148 (MQ=255)
                                                                          cTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGc                                                                          >  2:1123543/1‑148 (MQ=255)
                                                                           tGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGcc                                                                         <  2:2841226/148‑1 (MQ=255)
                                                                                   cATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGcgcgcg                                                                 <  1:624525/148‑1 (MQ=255)
                                                                                   cATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGcgcgcg                                                                 >  2:1635001/1‑148 (MQ=255)
                                                                                     ttACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgat                                                               <  1:2402174/148‑1 (MQ=255)
                                                                                     ttACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgat                                                               <  1:2403287/148‑1 (MQ=255)
                                                                                       aCATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatcg                                                             >  2:3318188/1‑148 (MQ=255)
                                                                                        cATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGcgatcga                                                            <  1:1757423/148‑1 (MQ=255)
                                                                                          tCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATa                                                          >  1:1049301/1‑148 (MQ=255)
                                                                                           ccGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATaa                                                         >  2:2280637/1‑148 (MQ=255)
                                                                                             gTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATaaaa                                                       <  1:869878/148‑1 (MQ=255)
                                                                                              tAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAAc                                                      >  1:1136801/1‑148 (MQ=255)
                                                                                                 ttGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGaa                                                   <  2:1165676/148‑1 (MQ=255)
                                                                                                  tGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGaaa                                                  >  2:67814/1‑148 (MQ=255)
                                                                                                  tGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGaaa                                                  <  2:2449257/148‑1 (MQ=255)
                                                                                                   gTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAat                                                 <  2:2860612/148‑1 (MQ=255)
                                                                                                      aaaCAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATc                                              >  2:424522/1‑148 (MQ=255)
                                                                                                        acaacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAg                                            <  2:881552/148‑1 (MQ=255)
                                                                                                         caacCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGc                                           <  1:788334/148‑1 (MQ=255)
                                                                                                           acCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTc                                         >  2:3435738/1‑148 (MQ=255)
                                                                                                           acCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTc                                         <  2:2320872/148‑1 (MQ=255)
                                                                                                            ccGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTcc                                        <  1:1966902/148‑1 (MQ=255)
                                                                                                                 gCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTg                                   <  1:137556/148‑1 (MQ=255)
                                                                                                                  ccTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTgg                                  >  1:2823215/1‑148 (MQ=255)
                                                                                                                   cTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTggg                                 <  1:2724664/148‑1 (MQ=255)
                                                                                                                      tCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGc                              <  2:1306691/148‑1 (MQ=255)
                                                                                                                      tCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGc                              <  2:3128565/148‑1 (MQ=255)
                                                                                                                             aGCGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACtt                       <  1:2142005/148‑1 (MQ=255)
                                                                                                                               cGTCATCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTAt                     >  2:104762/1‑148 (MQ=255)
                                                                                                                               cGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTAt                     <  2:1929837/148‑1 (MQ=255)
                                                                                                                               cGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTAt                     <  2:2970964/148‑1 (MQ=255)
                                                                                                                                gTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATc                    <  1:2093289/148‑1 (MQ=255)
                                                                                                                                 tCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCt                   <  2:3114349/148‑1 (MQ=255)
                                                                                                                                 tCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCt                   <  2:1850190/148‑1 (MQ=255)
                                                                                                                                  cATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCtt                  >  2:1590522/1‑148 (MQ=255)
                                                                                                                                  cATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCtt                  <  2:1698796/148‑1 (MQ=255)
                                                                                                                                   aTCGCTCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCACGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTc                 >  2:1170802/1‑148 (MQ=255)
                                                                                                                                     cGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCaa               <  1:134212/148‑1 (MQ=255)
                                                                                                                                      gCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCaaa              >  2:3496148/1‑148 (MQ=255)
                                                                                                                                      gCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCaaa              <  2:2085418/148‑1 (MQ=255)
                                                                                                                                        tCAATTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATg            <  2:3441696/148‑1 (MQ=255)
                                                                                                                                           aTTTTGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCa         <  1:1189624/148‑1 (MQ=255)
                                                                                                                                            ttttGGATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCAc        >  2:526215/1‑148 (MQ=255)
                                                                                                                                             tttGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACt       >  1:2660891/1‑148 (MQ=255)
                                                                                                                                                ggATATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTa    >  2:1442258/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGCCAGTACGCCAATTGCCAGCGCGCGATCGATAAAACGAAATATCAGCTCCAGTTGGGAGCCATACTTATCTTCAAATGCCACTGTATC  >  CP009273/1637311‑1637604

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: ‑