Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,800 G→T 100% 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,8000GT97.2% 446.1 / 4.2 143E297* (GAA→TAA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base T (69/70);  total (70/73)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 143 total aligned reads displayed.

GGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGA  >  CP009273/279658‑279944
                                                                                                                                              |                                                                                                                                                
gggTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAc                                                                                                                                             >  2:3443584/1‑148 (MQ=255)
 ggTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcg                                                                                                                                            >  2:1932735/1‑148 (MQ=255)
 ggTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCAGCTAAAAcg                                                                                                                                            >  2:1932743/1‑148 (MQ=255)
  gTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCCGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgc                                                                                                                                           >  2:1971121/1‑148 (MQ=255)
  gTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgc                                                                                                                                           >  2:1969460/1‑148 (MQ=255)
   tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                          <  1:1370962/148‑1 (MQ=255)
   tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                          >  1:2586839/1‑148 (MQ=255)
   tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                          <  2:1365337/148‑1 (MQ=255)
   tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                          >  2:2633999/1‑148 (MQ=255)
      ggCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATg                                                                                                                                       <  2:2712596/148‑1 (MQ=255)
       gcgcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATgg                                                                                                                                      >  2:3205014/1‑148 (MQ=255)
       gcgcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATgg                                                                                                                                      >  2:3205022/1‑148 (MQ=255)
         gcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTg                                                                                                                                    <  2:1540463/148‑1 (MQ=255)
          cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                   <  1:464320/148‑1 (MQ=255)
          cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                   <  1:2101420/148‑1 (MQ=255)
          cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                   <  1:2007425/148‑1 (MQ=255)
           tGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGAt                                                                                                                                  <  1:1753273/148‑1 (MQ=255)
            gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                 <  2:1283738/148‑1 (MQ=255)
            gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                 >  2:3009698/1‑148 (MQ=255)
            gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                 <  2:693180/148‑1 (MQ=255)
             ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATcc                                                                                                                                <  1:1073983/148‑1 (MQ=255)
             ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATcc                                                                                                                                <  1:650977/148‑1 (MQ=255)
                       cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                      <  2:135843/148‑1 (MQ=255)
                       cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                      >  1:226376/1‑148 (MQ=255)
                       cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                      >  1:226508/1‑148 (MQ=255)
                       cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                      >  1:2160028/1‑148 (MQ=255)
                       cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                      <  1:1540385/148‑1 (MQ=255)
                       cgcTGGCGCCGTCCGGGCAGGCGGAGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                      <  1:1540493/148‑1 (MQ=255)
                            gcgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCACTCCGCGCGCGCGGTCAGCGAGCTGGAAAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGAAGGTGATCCACGCGGCGTTcggcg                                                                                                                 >  1:1540304/1‑148 (MQ=255)
                             cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                                <  1:2115771/148‑1 (MQ=255)
                             cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcggcgg                                                                                                                >  2:2969582/1‑148 (MQ=255)
                               ccGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCt                                                                                                              >  2:20492/1‑148 (MQ=255)
                               ccGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCt                                                                                                              >  2:2217948/1‑148 (MQ=255)
                               ccGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCt                                                                                                              >  2:3353849/1‑148 (MQ=255)
                                cGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTc                                                                                                             <  1:1897899/148‑1 (MQ=255)
                                      ggCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCaa                                                                                                       <  2:3405987/148‑1 (MQ=255)
                                      ggCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCaa                                                                                                       >  2:71502/1‑148 (MQ=255)
                                          ggCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTa                                                                                                   <  1:1748548/148‑1 (MQ=255)
                                           gCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAc                                                                                                  >  1:400757/1‑148 (MQ=255)
                                             ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                                <  2:570993/148‑1 (MQ=255)
                                              gtgtGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActgc                                                                                               >  1:3474064/1‑148 (MQ=255)
                                              gtgtGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActgc                                                                                               >  1:3473993/1‑148 (MQ=255)
                                                gtgGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActgctg                                                                                             >  1:1894808/1‑148 (MQ=255)
                                                          aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                   <  1:3026006/148‑1 (MQ=255)
                                                            cgcccgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGcc                                                                                 <  2:408265/148‑1 (MQ=255)
                                                            cgcccgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGcc                                                                                 <  2:2366423/148‑1 (MQ=255)
                                                            cgcccgccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGcc                                                                                 <  2:408237/148‑1 (MQ=255)
                                                                  ccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGcc                                                                           <  2:750999/148‑1 (MQ=255)
                                                                    aGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCa                                                                         >  2:569177/1‑148 (MQ=255)
                                                                          cgcgcgcgGTCAGCTAGCTGGTTATCCGCGTCATCACCACTCGGGATATACTATCCGATAAATCCATCTAAATCTCGATGGTGATCCACGCGGCGTTCGGCGGCTCCAACAATTTACTGCTGCACATGCCGGCCATCGCCCTCgctgc                                                                   <  1:2184107/148‑1 (MQ=255)
                                                                              cgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGc                                                               >  1:1899844/1‑148 (MQ=255)
                                                                                    cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                         >  1:533252/1‑148 (MQ=255)
                                                                                    cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                         >  1:1274509/1‑148 (MQ=255)
                                                                                     aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                        >  1:3104133/1‑148 (MQ=255)
                                                                                     aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                        >  1:3103189/1‑148 (MQ=255)
                                                                                     aGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGa                                                        >  1:2581033/1‑148 (MQ=255)
                                                                                      gCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGAt                                                       >  1:3555758/1‑148 (MQ=255)
                                                                                        gAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGAGGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATcc                                                     <  2:1165563/148‑1 (MQ=255)
                                                                                          gCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCg                                                   >  1:1594490/1‑148 (MQ=255)
                                                                                              gATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACg                                               <  2:2657999/148‑1 (MQ=255)
                                                                                              gATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACg                                               <  1:1132623/148‑1 (MQ=255)
                                                                                               aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                              >  1:2737484/1‑148 (MQ=255)
                                                                                               aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                              <  1:1971121/148‑1 (MQ=255)
                                                                                               aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                              <  2:226508/148‑1 (MQ=255)
                                                                                               aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                              <  2:1885603/148‑1 (MQ=255)
                                                                                               aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                              <  2:2160028/148‑1 (MQ=255)
                                                                                                tAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGtt                                             >  1:2647611/1‑148 (MQ=255)
                                                                                                tAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGtt                                             >  1:2647349/1‑148 (MQ=255)
                                                                                                  gCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGCTGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGa                                           <  2:670800/148‑1 (MQ=255)
                                                                                                  gCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGa                                           <  1:3534408/148‑1 (MQ=255)
                                                                                                  gCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGa                                           <  2:669433/148‑1 (MQ=255)
                                                                                                     gcgGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCa                                        <  2:889915/148‑1 (MQ=255)
                                                                                                     gcgGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCa                                        <  2:891572/148‑1 (MQ=255)
                                                                                                     gcgGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCa                                        <  2:1620508/148‑1 (MQ=255)
                                                                                                     gcgGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGAGGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCa                                        <  2:668409/148‑1 (MQ=255)
                                                                                                         cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                    >  1:153153/1‑148 (MQ=255)
                                                                                                         cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                    >  1:1331631/1‑148 (MQ=255)
                                                                                                         cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                    >  1:918138/1‑148 (MQ=255)
                                                                                                         cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                    >  1:918147/1‑148 (MQ=255)
                                                                                                            caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                 >  1:2796146/1‑148 (MQ=255)
                                                                                                            caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                 <  1:1608198/148‑1 (MQ=255)
                                                                                                            caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                 >  1:2797941/1‑148 (MQ=255)
                                                                                                             accacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgc                                >  2:928129/1‑148 (MQ=255)
                                                                                                             accacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgc                                >  1:234653/1‑148 (MQ=255)
                                                                                                               cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                              <  2:1681090/148‑1 (MQ=255)
                                                                                                               cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                              >  1:2739763/1‑148 (MQ=255)
                                                                                                               cacGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                              <  1:2981745/148‑1 (MQ=255)
                                                                                                                 cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                            <  1:286216/148‑1 (MQ=255)
                                                                                                                     ggATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAc                        <  1:2187289/148‑1 (MQ=255)
                                                                                                                       atatCCTCCCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCg                      >  1:1464663/1‑148 (MQ=255)
                                                                                                                              ctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTg               <  2:1861069/148‑1 (MQ=255)
                                                                                                                               tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc              <  2:1042187/148‑1 (MQ=255)
                                                                                                                                 cGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcg            <  1:2633999/148‑1 (MQ=255)
                                                                                                                                   aTAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcgcg          >  1:2977562/1‑148 (MQ=255)
                                                                                                                                   aTAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcgcg          >  1:3428197/1‑148 (MQ=255)
                                                                                                                                   aTAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcgcg          >  1:3428258/1‑148 (MQ=255)
                                                                                                                                    tAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGt         >  1:3188234/1‑148 (MQ=255)
                                                                                                                                       aGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg      <  1:418173/148‑1 (MQ=255)
                                                                                                                                         ccATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGt    >  1:3068004/1‑148 (MQ=255)
                                                                                                                                           aTCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGa  <  1:2375255/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGA  >  CP009273/279658‑279944

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