Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,800 G→T 86.3% E297* (GAA→TAA)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,8000GT86.3% 231.2 / 36.8 102E297* (GAA→TAA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (4/10);  new base T (44/44);  total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 102 total aligned reads displayed.

TGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGT  >  CP009273/279657‑279942
                                                                                                                                               |                                                                                                                                              
tGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaaa                                                                                                                                            <  1:3094664/148‑1 (MQ=255)
  ggTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAAcg                                                                                                                                          <  2:1602551/148‑1 (MQ=255)
    tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                        <  2:1304813/148‑1 (MQ=255)
     cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGa                                                                                                                                       <  2:820890/148‑1 (MQ=255)
     cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGa                                                                                                                                       <  2:2837160/148‑1 (MQ=255)
        gcgcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATgg                                                                                                                                    >  2:1149318/1‑148 (MQ=255)
          gcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTg                                                                                                                                  <  1:2872085/148‑1 (MQ=255)
           cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                 <  2:1716000/148‑1 (MQ=255)
           cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGa                                                                                                                                 <  1:1583119/148‑1 (MQ=255)
           cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGa                                                                                                                                 <  1:3142343/148‑1 (MQ=255)
             gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGTCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATc                                                                                                                               <  2:2670840/148‑1 (MQ=255)
             gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATc                                                                                                                               <  2:2670836/148‑1 (MQ=255)
              ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATcc                                                                                                                              <  2:1866559/148‑1 (MQ=255)
               cGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCTATCTAAAACGCGATGGTGATCCa                                                                                                                             <  2:1385377/148‑1 (MQ=255)
                   cTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCAcgcg                                                                                                                         >  2:438653/1‑148 (MQ=255)
                    tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcgg                                                                                                                        >  2:1224248/1‑148 (MQ=255)
                    tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcgg                                                                                                                        >  1:2970160/1‑148 (MQ=255)
                    tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcgg                                                                                                                        >  1:2866376/1‑148 (MQ=255)
                     ccGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggc                                                                                                                       >  1:2009079/1‑148 (MQ=255)
                     ccGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggc                                                                                                                       >  1:2008987/1‑148 (MQ=255)
                      cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggcg                                                                                                                      <  2:1987840/148‑1 (MQ=255)
                       gcgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGt                                                                                                                     <  2:1379053/148‑1 (MQ=255)
                        cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                    <  2:1526682/148‑1 (MQ=255)
                            ggCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggc                                                                                                                >  1:1582074/1‑148 (MQ=255)
                              cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                              <  2:2356255/148‑1 (MQ=255)
                              cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                              <  2:2356635/148‑1 (MQ=255)
                              cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                              <  1:2352651/148‑1 (MQ=255)
                                           ggCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTa                                                                                                 >  1:2100596/1‑148 (MQ=255)
                                              ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                              <  1:2449390/148‑1 (MQ=255)
                                              ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                              <  2:1637847/148‑1 (MQ=255)
                                                     cTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcaca                                                                                       <  1:1805971/148‑1 (MQ=255)
                                                        gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                    <  1:2635568/148‑1 (MQ=255)
                                                           aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                 >  1:1970327/1‑148 (MQ=255)
                                                           aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                 >  1:1969843/1‑148 (MQ=255)
                                                                 cgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCg                                                                           >  2:157254/1‑148 (MQ=255)
                                                                 cgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCg                                                                           >  2:157189/1‑148 (MQ=255)
                                                                  gccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGc                                                                          <  1:2936242/148‑1 (MQ=255)
                                                                      gTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAc                                                                      >  1:2036098/1‑148 (MQ=255)
                                                                       tCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAcg                                                                     <  2:1689833/148‑1 (MQ=255)
                                                                       tCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAcg                                                                     <  2:1689184/148‑1 (MQ=255)
                                                                            gcgcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACgcggcg                                                                <  2:2099245/148‑1 (MQ=255)
                                                                            gcgcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACgcggcg                                                                <  2:2100182/148‑1 (MQ=255)
                                                                            gcgcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACgcggcg                                                                <  2:2101490/148‑1 (MQ=255)
                                                                              gcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCggg                                                              <  1:1118897/148‑1 (MQ=255)
                                                                               cgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGc                                                             <  1:2756625/148‑1 (MQ=255)
                                                                                 cgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTg                                                           >  2:819971/1‑148 (MQ=255)
                                                                                 cgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTg                                                           <  2:1658407/148‑1 (MQ=255)
                                                                                   gTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                         >  1:431183/1‑148 (MQ=255)
                                                                                      aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                      <  1:2507861/148‑1 (MQ=255)
                                                                                      aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                      <  1:3137279/148‑1 (MQ=255)
                                                                                        cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCAGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                    >  2:2311030/1‑148 (MQ=255)
                                                                                          aGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATccc                                                  >  2:912624/1‑148 (MQ=255)
                                                                                                 tAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGtt                                           >  2:1153016/1‑148 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:1018405/148‑1 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:2547674/148‑1 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:1018402/148‑1 (MQ=255)
                                                                                                  aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                          <  1:1018388/148‑1 (MQ=255)
                                                                                                   gCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGa                                         >  2:497781/1‑148 (MQ=255)
                                                                                                    ccGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAg                                        <  1:3210230/148‑1 (MQ=255)
                                                                                                    ccGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAg                                        <  1:438653/148‑1 (MQ=255)
                                                                                                     cgcgGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGc                                       <  2:1666754/148‑1 (MQ=255)
                                                                                                         gCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTg                                   >  2:485566/1‑148 (MQ=255)
                                                                                                         gCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTg                                   >  2:485546/1‑148 (MQ=255)
                                                                                                            tCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAc                                >  1:3234055/1‑148 (MQ=255)
                                                                                                            tCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAc                                >  1:2399170/1‑148 (MQ=255)
                                                                                                             caccacGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                               >  2:806803/1‑148 (MQ=255)
                                                                                                             caccacGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                               >  2:2712940/1‑148 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             >  1:503742/1‑148 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             <  2:1657645/148‑1 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             >  1:503934/1‑148 (MQ=255)
                                                                                                               ccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgcg                             <  2:3295176/148‑1 (MQ=255)
                                                                                                                cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                            >  1:2934967/1‑148 (MQ=255)
                                                                                                                cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                            >  1:2934960/1‑148 (MQ=255)
                                                                                                                cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                            >  1:2934849/1‑148 (MQ=255)
                                                                                                                 acGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCa                           <  1:486844/148‑1 (MQ=255)
                                                                                                                  cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                          <  1:2529584/148‑1 (MQ=255)
                                                                                                                    cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                        >  2:1415131/1‑148 (MQ=255)
                                                                                                                      ggATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTGACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAc                      <  1:1589061/148‑1 (MQ=255)
                                                                                                                          atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                  >  2:2588082/1‑148 (MQ=255)
                                                                                                                          atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                  <  2:1499312/148‑1 (MQ=255)
                                                                                                                          atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                  <  2:1499169/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  1:550625/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  1:1925397/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 >  1:2989664/1‑148 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  1:2631061/148‑1 (MQ=255)
                                                                                                                           tCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                 <  2:3075994/148‑1 (MQ=255)
                                                                                                                            ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCCAGTTACTGCTGCACATTCCGGCCATAGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                >  1:2653051/1‑148 (MQ=255)
                                                                                                                            ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                >  1:1875963/1‑148 (MQ=255)
                                                                                                                              tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt              >  1:3146620/1‑148 (MQ=255)
                                                                                                                              tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt              >  1:1258014/1‑148 (MQ=255)
                                                                                                                               ctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTg             <  2:2808343/148‑1 (MQ=255)
                                                                                                                               ctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTg             >  2:2227609/1‑148 (MQ=255)
                                                                                                                                tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc            >  2:332385/1‑148 (MQ=255)
                                                                                                                                tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc            >  1:1866238/1‑148 (MQ=255)
                                                                                                                                    aTAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcgcg        >  2:3183986/1‑148 (MQ=255)
                                                                                                                                     tAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGt       <  2:2202076/148‑1 (MQ=255)
                                                                                                                                      aaaGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTc      >  2:2514478/1‑148 (MQ=255)
                                                                                                                                       aaGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCt     >  1:524233/1‑148 (MQ=255)
                                                                                                                                         gCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTgg   >  1:2579148/1‑148 (MQ=255)
                                                                                                                                          ccATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGt  >  2:1749722/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGT  >  CP009273/279657‑279942

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