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,8000GT96.9% 387.6 / 4.3 127E297* (GAA→TAA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (52/71);  total (54/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 127 total aligned reads displayed.

GCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTG  >  CP009273/279655‑279943
                                                                                                                                                 |                                                                                                                                               
gcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaa                                                                                                                                               <  2:1563292/148‑1 (MQ=255)
gcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaa                                                                                                                                               <  2:1563299/148‑1 (MQ=255)
gcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaa                                                                                                                                               <  2:3130368/148‑1 (MQ=255)
 cTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaa                                                                                                                                              <  1:2524859/148‑1 (MQ=255)
 cTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaa                                                                                                                                              <  1:2524860/148‑1 (MQ=255)
 cTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaa                                                                                                                                              <  1:2526469/148‑1 (MQ=255)
  tGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaaa                                                                                                                                             <  1:1232581/148‑1 (MQ=255)
      tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                         <  2:1915160/148‑1 (MQ=255)
       cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGa                                                                                                                                        <  1:3325367/148‑1 (MQ=255)
       cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGa                                                                                                                                        <  1:3208521/148‑1 (MQ=255)
         ggCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATg                                                                                                                                      <  2:1643514/148‑1 (MQ=255)
         ggCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATg                                                                                                                                      >  2:2706577/1‑148 (MQ=255)
            gcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTg                                                                                                                                   <  1:2161441/148‑1 (MQ=255)
             cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                  <  1:2841512/148‑1 (MQ=255)
               gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                <  1:1926408/148‑1 (MQ=255)
               gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                <  2:2993661/148‑1 (MQ=255)
               gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                <  1:1431759/148‑1 (MQ=255)
                 cGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCa                                                                                                                              <  1:2789103/148‑1 (MQ=255)
                      tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACgcgg                                                                                                                         <  1:859295/148‑1 (MQ=255)
                      tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACgcgg                                                                                                                         <  1:981877/148‑1 (MQ=255)
                        cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCATCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggcg                                                                                                                       <  2:1317823/148‑1 (MQ=255)
                        cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggcg                                                                                                                       <  1:138919/148‑1 (MQ=255)
                         gcgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGt                                                                                                                      <  2:3130990/148‑1 (MQ=255)
                          cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                     >  1:2741139/1‑148 (MQ=255)
                          cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                     <  1:427333/148‑1 (MQ=255)
                           gcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTc                                                                                                                    <  2:1497848/148‑1 (MQ=255)
                             tGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcgg                                                                                                                  <  2:252565/148‑1 (MQ=255)
                                      ccGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                         >  2:3363296/1‑148 (MQ=255)
                                      ccGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                         >  2:3363294/1‑148 (MQ=255)
                                         ggCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCaa                                                                                                      <  1:1816923/148‑1 (MQ=255)
                                              gCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAc                                                                                                 >  1:1830574/1‑148 (MQ=255)
                                              gCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAc                                                                                                 >  1:2804014/1‑148 (MQ=255)
                                                ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                               <  2:1073049/148‑1 (MQ=255)
                                                 gtgtGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActgc                                                                                              >  1:1972745/1‑148 (MQ=255)
                                                    tgGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACtgctgc                                                                                           >  1:3205092/1‑148 (MQ=255)
                                                      gCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcac                                                                                         <  2:359931/148‑1 (MQ=255)
                                                      gCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcac                                                                                         <  2:2128577/148‑1 (MQ=255)
                                                      gCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcac                                                                                         <  1:2041013/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  1:2278625/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  2:2408079/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  2:645279/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  1:1071772/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  2:2408130/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     >  2:3360287/1‑148 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  1:1071550/148‑1 (MQ=255)
                                                             aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                  >  1:1556292/1‑148 (MQ=255)
                                                                     ccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGcc                                                                          <  2:158392/148‑1 (MQ=255)
                                                                       aGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCa                                                                        <  2:2052313/148‑1 (MQ=255)
                                                                                gcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCggg                                                               <  1:1757082/148‑1 (MQ=255)
                                                                                 cgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGc                                                              <  2:963586/148‑1 (MQ=255)
                                                                                       cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                        >  1:609934/1‑148 (MQ=255)
                                                                                       cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                        <  2:2779261/148‑1 (MQ=255)
                                                                                       cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                        >  1:2372820/1‑148 (MQ=255)
                                                                                         gCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGAt                                                      >  1:1590004/1‑148 (MQ=255)
                                                                                          cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                     <  2:600345/148‑1 (MQ=255)
                                                                                          cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                     <  1:1023389/148‑1 (MQ=255)
                                                                                          cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAAGTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                     >  2:774789/1‑148 (MQ=255)
                                                                                                 gATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACg                                              <  1:3174245/148‑1 (MQ=255)
                                                                                                 gATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACg                                              <  1:3174270/148‑1 (MQ=255)
                                                                                                  aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                             <  2:1814004/148‑1 (MQ=255)
                                                                                                    aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                           >  1:1092308/1‑148 (MQ=255)
                                                                                                      ccGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAg                                         >  2:801840/1‑148 (MQ=255)
                                                                                                            cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                   >  2:561949/1‑148 (MQ=255)
                                                                                                            cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                   >  1:762285/1‑148 (MQ=255)
                                                                                                             aTCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGa                                  >  2:2898031/1‑148 (MQ=255)
                                                                                                               caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGTTGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                <  2:1409127/148‑1 (MQ=255)
                                                                                                               caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                <  2:3124939/148‑1 (MQ=255)
                                                                                                                accacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgc                               >  1:1391642/1‑148 (MQ=255)
                                                                                                                  cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                             <  2:2717558/148‑1 (MQ=255)
                                                                                                                  cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                             <  2:2717632/148‑1 (MQ=255)
                                                                                                                    cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                           <  2:1830574/148‑1 (MQ=255)
                                                                                                                    cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                           <  1:1607365/148‑1 (MQ=255)
                                                                                                                      cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                         >  1:3110161/1‑148 (MQ=255)
                                                                                                                      cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                         >  1:2334077/1‑148 (MQ=255)
                                                                                                                      cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                         >  2:3022812/1‑148 (MQ=255)
                                                                                                                         gATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAcc                      >  2:665398/1‑148 (MQ=255)
                                                                                                                         gATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAcc                      >  2:1067703/1‑148 (MQ=255)
                                                                                                                         gATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAcc                      >  2:665799/1‑148 (MQ=255)
                                                                                                                            atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                   >  1:1698607/1‑148 (MQ=255)
                                                                                                                             tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                  >  1:713015/1‑148 (MQ=255)
                                                                                                                             tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                  >  1:713200/1‑148 (MQ=255)
                                                                                                                             tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                  <  2:1216663/148‑1 (MQ=255)
                                                                                                                              ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                 >  1:3464529/1‑148 (MQ=255)
                                                                                                                              ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                 >  1:3464482/1‑148 (MQ=255)
                                                                                                                              ccTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                 >  2:1048007/1‑148 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                <  2:897448/148‑1 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                >  2:2425154/1‑148 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                <  1:872997/148‑1 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATGTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                <  1:873242/148‑1 (MQ=255)
                                                                                                                                tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt               >  1:3282588/1‑148 (MQ=255)
                                                                                                                                tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt               >  1:3282592/1‑148 (MQ=255)
                                                                                                                                  tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc             <  2:3171471/148‑1 (MQ=255)
                                                                                                                                   ccGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGcc            >  2:1314707/1‑148 (MQ=255)
                                                                                                                                       tAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGt        >  1:3118594/1‑148 (MQ=255)
                                                                                                                                        aaaGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTc       >  2:610456/1‑148 (MQ=255)
                                                                                                                                          aGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg     <  1:1471709/148‑1 (MQ=255)
                                                                                                                                          aGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg     <  2:1972745/148‑1 (MQ=255)
                                                                                                                                          aGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg     >  1:411351/1‑148 (MQ=255)
                                                                                                                                             cATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTg  >  1:3075/1‑148 (MQ=255)
                                                                                                                                             cATCTAAAACGCGATGATGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTg  >  1:2810479/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTG  >  CP009273/279655‑279943

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: