Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0680AC25.7% 218.4 / 20.3 144intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (47/60);  new base C (5/32);  total (52/92)
Fisher's exact test for biased strand distribution p-value = 7.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 147 total aligned reads displayed.

ACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAT  >  CP009273/160918‑161218
                                                                                                                                                       |                                                                                                                                                       
aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                         >  2:1564149/1‑151 (MQ=255)
 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                        >  1:1628889/1‑151 (MQ=255)
  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGTGGGCCGGGGTGAGGGCATc                                                                                                                                                       >  1:3727805/1‑151 (MQ=255)
  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                       <  1:3508743/151‑1 (MQ=255)
  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                       >  1:2419648/1‑151 (MQ=255)
    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                     <  1:4360617/151‑1 (MQ=255)
      aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                                   >  1:5263735/1‑151 (MQ=255)
       ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                  >  1:46140/1‑151 (MQ=255)
         tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATACGTAAGGTTGGTTTTCTCCCACTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCTgcgcgc                                                                                                                                                >  2:3943020/1‑151 (MQ=255)
            aaaaaGAGATGAAAGGGCGTTATCCCAATCTTGTCTGGGCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTGTTCGGAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGCGAGGGCCGGGGTGAGGGCCTCTGCGCGCTCg                                                                                                                                             >  2:2979658/1‑151 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                            >  2:2340266/1‑151 (MQ=255)
             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                            >  1:4273624/1‑151 (MQ=255)
              aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                                           >  2:420853/1‑151 (MQ=255)
              aaaGAGATGAAAGGGCGTTATCCCAAACATGACTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGAACGtg                                                                                                                                           >  1:30952/1‑150 (MQ=255)
                agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGtgga                                                                                                                                         >  1:779686/1‑148 (MQ=255)
                agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggca                                                                                                                                         >  1:606239/1‑147 (MQ=255)
                 gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGCCGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggaa                                                                                                                                        >  2:611247/1‑147 (MQ=255)
                  agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGGGGGCATCAGCGCGCACGtggacc                                                                                                                                       >  1:1476872/1‑146 (MQ=255)
                       aacgggCGTTATCCAAAACATGTCTGGCGGTACGCCCCGGAAATTACGGCACCGACGCGAGGGACTAAAATGTTTTCGTATGTTTGGTTTTCTCCCTCTCC‑CTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                  <  2:1054917/148‑1 (MQ=25)
                       aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacac                                                                                                                                  >  1:2232178/1‑151 (MQ=255)
                           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacacccgt                                                                                                                              >  2:4286471/1‑149 (MQ=255)
                             ctttAACCCAACCTTGGCTGGCGGCACGCCCCGCAAAATCCTGCACCCCCGCCCCGGACCAAAATTTATTGTTCTGGTTGTTTCTCTCCCTCCCC‑GTGTGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                            <  1:1204146/149‑1 (MQ=255)
                             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCAAAc                                                                                                                            >  2:5297160/1‑151 (MQ=255)
                                  tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAAcccaca                                                                                                                       >  2:2224500/1‑148 (MQ=255)
                                   cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCAAAaccacac                                                                                                                      >  2:168311/1‑147 (MQ=255)
                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGACCCGCACCCAAAACCAC‑Acccc                                                                                                                  >  1:3385325/1‑150 (MQ=37)
                                         cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCTC‑ACCCTCaa                                                                                                               >  2:1297069/1‑151 (MQ=255)
                                            gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACACTAACCCAC‑TCCCAGAAggg                                                                                                            >  2:1251531/1‑151 (MQ=255)
                                              cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGCCCCGCCCCCGAACCCAC‑TCCCGCATggggg                                                                                                          >  1:829127/1‑149 (MQ=25)
                                                ggCGGGACGACCCGGCAAATACGGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGTTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGagag                                                                                                        <  2:2981830/151‑1 (MQ=255)
                                                 gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCTTATGGTTGGTTTTCTCCCTCACC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAgg                                                                                                       <  2:1153133/151‑1 (MQ=255)
                                                   cGGAGGCCCCGCCAAAAACTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGTTTTTCCCCCCCCCC‑CTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAgggg                                                                                                     <  2:2565359/151‑1 (MQ=255)
                                                    ggccgacCCGGCAAATACTGCACCGCCGCGACGGACGAAAATTTATTCGTACGGTTGGTTTTCTCCCTCTCC‑CTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGa                                                                                                    <  2:4187990/148‑1 (MQ=38)
                                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCGAAACCTC‑TCCCGCAAGGGAGAGGGGa                                                                                                    >  2:349675/1‑151 (MQ=255)
                                                        gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGACCACCCCCAAAACCCC‑ACCCGCAAGGGAGCGGGGGCCGa                                                                                                >  1:2093058/1‑151 (MQ=25)
                                                         ccccttaaaaTCCGGCCCCGAGCCGACGGACTACATAGTATTCGTAGTGTTGGTTGTCCCCCTCTCC‑GGGTGGGAGCGGGCGGGGGGGCGGGCATCAGCGCGCACGTTCCCCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGAt                                                                                               <  1:5002913/144‑1 (MQ=255)
                                                         acCCGGCAAATACTGCCCCGCCGCGACGGACGAAAATGTATTCTTACGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGAt                                                                                               <  1:782009/151‑1 (MQ=38)
                                                         acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAGCCCCCAACCCCCAAGGGGGGGGGGGCCgt                                                                                                >  1:932651/1‑150 (MQ=25)
                                                          cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCGAACCCCC‑TCCCCCCAGGGCGGGGGGGCCGaag                                                                                              >  2:1492015/1‑149 (MQ=25)
                                                          cccGGAAAATCCTCCCCCGCCGCGACGCACCAAAATTTTTCCTTATGGTTGGTTTTCCCCCCCCCC‑CTGGGGGCGTGGGCGGGGGGTTGGCCATCAGCGCGCACGGTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATc                                                                                              <  1:4175831/151‑1 (MQ=255)
                                                           ccGGCAAATCCTGCCCCGCCGCGCCGGACGAAAAAGTATTCGTACGGTTGGTGTTCGCCCACTCC‑CGGTGGTTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCg                                                                                             <  2:2743349/151‑1 (MQ=255)
                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCAAACCCAA‑ACCCACAAGGGTGGGGGGGCCCAtac                                                                                             >  1:310845/1‑149 (MQ=25)
                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCGACCCCAAACCCC‑GACCACACTGGAGGGGGGAGCGGTCg                                                                                             >  1:5140122/1‑151 (MQ=21)
                                                            cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACCACCGAACCCGC‑GCCCACCAGGGCGAGGGGACCGAtaca                                                                                            >  2:3843877/1‑148 (MQ=25)
                                                                aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAACCAC‑GCCCCACAGGGGGGGGGGGCCCATCGAggcc                                                                                        >  1:404233/1‑148 (MQ=25)
                                                                 aaTACTCCACCGACGCTACGGCCTAAAAATTTTTCGTAAGGTTGGTTTTCTCCCTCCCC‑CTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                       <  1:1327006/151‑1 (MQ=38)
                                                                  aTACGGCACCGCCGCGACGGACTAAAATTTATTGTTATGGTTGGTTTTCTCCCTCTCC‑CGGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACaa                                                                                      <  2:460429/151‑1 (MQ=255)
                                                                    aCTGCACCGCGGCGCCGCACGAAAATTAATTCGTACGGTTGTTTTTCTCCCACTCC‑CGGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                    <  1:3142034/151‑1 (MQ=255)
                                                                    aCTGCACCGACGCCACGTACTAAAAAGTATTCGACCGGTGGGTTTACTCCCTCTCC‑CTGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                    <  2:1982414/151‑1 (MQ=255)
                                                                     cggCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                   <  1:5398527/149‑1 (MQ=255)
                                                                     cTGCACCGACGCGACGGACGAACAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCAACCCCCCCCCAACCCAC‑ACCCGCCAGGGGGGGGGGGGCGAGAGAGCACcaaat                                                                                   >  1:2641544/1‑149 (MQ=16)
                                                                     cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                   <  2:4892740/151‑1 (MQ=255)
                                                                        cACCGACGCTACGGACTAAAATTTTTTTTTTTGTTTGTTTTTCCCCCTCCCC‑GTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                <  1:2356649/151‑1 (MQ=255)
                                                                         aCCGCCGCGAGGGACTAAAAATTTTTCTAACGGTTTTTTTTCCCCCTCTCC‑CTGTGGGATTGGGCCGGGGTTCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  1:308411/151‑1 (MQ=255)
                                                                         aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCAC‑CCCCTCAAGTGCGTTGGGCCCCATCCAGCTCAAAatttta                                                                               >  2:2113229/1‑149 (MQ=21)
                                                                             acgcgacgGCCGAAAAAGTTTTCGTAAGGTTGTTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                           <  1:4504531/151‑1 (MQ=255)
                                                                                cgacgGACTAAAATTTTTTCTTATGTTTGGTTTTCTCCCTCTCC‑CGGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                        <  1:3178018/151‑1 (MQ=255)
                                                                                 gacggtctAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCC‑CTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                       <  2:2700790/151‑1 (MQ=38)
                                                                                  acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGACCTCACCCTAACCCTC‑TCCCACAAGGGAGAGGGGACCCATCGAGCCCAAATTTTGTGAGAAAtct                                                                      >  2:5146393/1‑151 (MQ=37)
                                                                                     gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTAC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCAC‑TCCCACAAGGGGGCGGGGACCGATCGAGCCCAAATTTTGAGAGAAAtcttat                                                                   >  2:624377/1‑149 (MQ=37)
                                                                                      actaaaatGTTTTCTTTAGGTTGGTTTT‑TGCCTCCCCGTTTTGGTTGTGGGCGGGGGGGAGGCCCTCAGCGCGCTCGTTCACCCTCACCATAACCATC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATAtcttctt                                                                  <  1:5259648/148‑1 (MQ=255)
                                                                                      actaaaaaTTATTCTTATGGTTGTTTTACTCCAACTCC‑CTGTGGGAGATGTCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGTGAGAGGGGACAGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                  <  2:2573624/148‑1 (MQ=255)
                                                                                      acgAAAAAGAATTCGTAAGGTTGGTTTTAGCCATCTCC‑CTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                  <  2:2614082/151‑1 (MQ=255)
                                                                                         aaaaaGTATTCGTCAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTg                                                               <  2:4467650/151‑1 (MQ=255)
                                                                                            caGTATTCGTACGGTTGGATTTCTCCCACACC‑CTGTCCGACAGGGCCGGGCTGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGACt                                                            <  2:382451/150‑1 (MQ=255)
                                                                                            aattATTCCAAAGGTGGGTTTACTCCCTCCCC‑CTGTGGTAGTGGGCCGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                            <  1:1741807/148‑1 (MQ=255)
                                                                                            aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCAAAACCAC‑GCCCACCAGGGAGGGGGGGCCGCGAGAGCCCAAAAATTGTGAGATAtcttattagttgt                                                            >  1:3526007/1‑145 (MQ=16)
                                                                                                    ttatggttgtctTACTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCTAGGGAGAGGGGACCGATCGACCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTACCag                                                    <  2:2746875/147‑1 (MQ=38)
                                                                                                    gtaggtttgtttTTCTCCCTCACC‑TGGTGGTAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                    <  1:2535106/145‑1 (MQ=255)
                                                                                                    gTATGGTTGTTTTTCTCCCTCCCC‑CTGGGGGTGTGGGCGGGGTTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                    <  1:2535123/151‑1 (MQ=255)
                                                                                                    gTAAGTTTGGTTTTCTCCCTCTCC‑CTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                    <  1:5187448/151‑1 (MQ=255)
                                                                                                    ctacgtttgtttCTCACCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                    <  1:5079871/145‑1 (MQ=255)
                                                                                                      atggttgtttTTCTCCCTCTCC‑CTGTGGGAGTGGCCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                  <  1:4142067/149‑1 (MQ=39)
                                                                                                           tggttTTCTCCCTCTCC‑CTGGGGTAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                             <  2:1613867/151‑1 (MQ=255)
                                                                                                            ggttTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCACACCCAAAACCAC‑TCCCACAAGGGGGTGGGGGCCGCTCGAGCCCAAAATTTTTGTTATAAATTAATCTGTGTTGTAACAGAATAAAAg                                            >  2:1094917/1‑151 (MQ=12)
                                                                                                             gttTTCTCCCACACC‑CTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGAACGAGCACACATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                           <  2:4967351/151‑1 (MQ=255)
                                                                                                                   tccctctcc‑ctGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                     <  2:879136/151‑1 (MQ=255)
                                                                                                                    ccctctcc‑ctGTGGTAGCGGGCCGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGATGAACAGAAAATCgg                                    <  2:2904532/151‑1 (MQ=255)
                                                                                                                    ccctctcc‑ctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                    <  2:3500732/151‑1 (MQ=255)
                                                                                                                    ccctctcc‑caGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                    <  2:416890/151‑1 (MQ=255)
                                                                                                                    ccctcacc‑ccGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑ACCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                    <  2:3500600/151‑1 (MQ=255)
                                                                                                                        cgcc‑ctGGGGGAGTGGGCCGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                                <  1:2760087/149‑1 (MQ=255)
                                                                                                                         tcc‑ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCGC‑GCCCACCAGGGGGGGGGGACCGATAGGGCCCCAATTTTGTGTGATATATTATTCTTGGTTGGTACACAAGAACAGAAAATCGGGCgtt                               >  2:89722/1‑151 (MQ=16)
                                                                                                                         tcc‑ctGTGGCTTAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                               <  2:1618745/151‑1 (MQ=255)
                                                                                                                          cc‑ctGTGGGAGAGGGCCGGGGAGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                              <  2:214110/151‑1 (MQ=255)
                                                                                                                             ctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                            <  2:2489045/151‑1 (MQ=255)
                                                                                                                             ctGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCAC‑TCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGACAATCGGGCTTTTgc                            <  1:3312377/151‑1 (MQ=255)
                                                                                                                                    tgagGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGac                     <  1:3780842/150‑2 (MQ=255)
                                                                                                                                      agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                   <  2:638533/151‑1 (MQ=255)
                                                                                                                                          ccGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc               <  2:2919404/151‑1 (MQ=255)
                                                                                                                                          ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc               <  2:4437612/151‑1 (MQ=255)
                                                                                                                                           cGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCGAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg              <  2:3782918/151‑1 (MQ=255)
                                                                                                                                            ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg             <  1:4108618/151‑1 (MQ=255)
                                                                                                                                               gTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga          <  2:4154261/151‑1 (MQ=255)
                                                                                                                                                tGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCAC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCAGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa         <  1:1778581/151‑1 (MQ=255)
                                                                                                                                                   gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa      <  2:3205380/151‑1 (MQ=255)
                                                                                                                                                   gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa      <  1:116575/151‑1 (MQ=255)
                                                                                                                                                    ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg     <  2:4094882/151‑1 (MQ=255)
                                                                                                                                                     gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc    <  2:4277931/151‑1 (MQ=255)
                                                                                                                                                       aTCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt  <  1:501588/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
ACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCC‑CTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTC‑TCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAT  >  CP009273/160918‑161218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑