Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 910,491 G→C 8.2% K333N (AAG→AAC) ‡ fruB → PTS fructose transporter subunit EI/HPr/IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947910,4910GC 8.2% 193.9 / 6.7 61K333N (AAG→AAC) ‡fruBPTS fructose transporter subunit EI/HPr/IIA
Reads supporting (aligned to +/- strand):  ref base G (21/35);  new base C (3/2);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCACGGGTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTACCGCAGAGCG  >  NC_002947/910345‑910630
                                                                                                                                                  |                                                                                                                                           
gcACGGGTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAgct                                                                                                                                           <  7:238524/149‑1 (MQ=255)
   cGGGTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGTGAAATGCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAgctgct                                                                                                                                        <  3:257988/149‑1 (MQ=255)
    gggTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTc                                                                                                                                       <  6:96873/149‑1 (MQ=255)
     ggTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCa                                                                                                                                      <  4:83213/149‑1 (MQ=255)
      gTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAg                                                                                                                                     <  5:232128/149‑1 (MQ=255)
       tCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGc                                                                                                                                    <  7:60399/149‑1 (MQ=255)
       gttctcTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAGCTTGCTCAgg                                                                                                                                    <  3:319742/144‑2 (MQ=255)
          gctgggCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTc                                                                                                                                 <  1:144378/145‑1 (MQ=255)
            cTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCgg                                                                                                                               >  7:45307/1‑149 (MQ=255)
             tGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGACTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAGCTTGCTCAgg                                                                                                                                    >  4:319742/1‑142 (MQ=255)
                 cAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCAc                                                                                                                          >  4:160263/1‑149 (MQ=255)
                         aTGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGt                                                                                                                  <  3:304488/149‑1 (MQ=255)
                            gCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCa                                                                                                                <  5:121219/148‑1 (MQ=255)
                              ttGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAgg                                                                                                              >  3:17599/1‑148 (MQ=255)
                                 cATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCAc                                                                                                           <  5:247851/148‑1 (MQ=255)
                                   tGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTgg                                                                                                        <  3:160263/149‑1 (MQ=255)
                                      gcgACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAAc                                                                                                     <  6:39028/149‑1 (MQ=255)
                                       cgACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACt                                                                                                    <  5:236598/149‑1 (MQ=255)
                                        gACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTg                                                                                                   <  5:348966/149‑1 (MQ=255)
                                           cGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTc                                                                                                <  6:350084/149‑1 (MQ=255)
                                           cGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTc                                                                                                <  5:100901/149‑1 (MQ=255)
                                                aaGGTGCTGGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCg                                                                                            >  8:207958/1‑148 (MQ=255)
                                                      cTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGa                                                                                      <  8:217608/148‑1 (MQ=255)
                                                         gCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCAt                                                                                  >  7:316912/1‑149 (MQ=255)
                                                           acaACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCg                                                                                >  7:167992/1‑149 (MQ=255)
                                                           acaACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCg                                                                                <  8:316912/149‑1 (MQ=255)
                                                                  gCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGcc                                                                           <  5:190096/147‑1 (MQ=255)
                                                                    cAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGa                                                                                      >  7:5814/1‑134 (MQ=255)
                                                                    cAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGa                                                                                      <  8:5814/134‑1 (MQ=255)
                                                                    cAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATg                                                                       >  2:162296/1‑149 (MQ=255)
                                                                     aaaGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGaa                                                                                          <  5:73608/129‑1 (MQ=255)
                                                                     aaaGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGaa                                                                                          >  6:73608/1‑129 (MQ=255)
                                                                     aaaGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGc                                                                      >  4:280638/1‑149 (MQ=255)
                                                                            ttGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCg                                                               <  5:56587/149‑1 (MQ=255)
                                                                            ttGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCg                                                               >  3:345356/1‑149 (MQ=255)
                                                                                 ggggAAATCCGTGTGCGCCTGGTCGTTAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGtgct                                                          <  8:167992/149‑1 (MQ=255)
                                                                                   ggAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGtgctgc                                                        <  1:162296/149‑1 (MQ=255)
                                                                                    gAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTgctgct                                                       <  3:280638/149‑1 (MQ=255)
                                                                                    gAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTgctgct                                                       <  4:2246/149‑1 (MQ=255)
                                                                                      aaTCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTgctgg                                                     >  5:279386/1‑149 (MQ=255)
                                                                                      aaTCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTgctgg                                                     <  3:369043/149‑1 (MQ=255)
                                                                                              tgCGCCTGGTCGTTAGCGAGCAGCCAGCCGTATCGGTTAAGAGCCTGAGGAAGCTGCTCAGCCTCGGCGGACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTAt                                              <  7:270545/148‑1 (MQ=255)
                                                                                                cgcCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCg                                            >  3:152055/1‑148 (MQ=255)
                                                                                                   cTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAAc                                         >  2:357890/1‑148 (MQ=255)
                                                                                                   cTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAAc                                         >  3:27475/1‑148 (MQ=255)
                                                                                                       tCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGc                                    <  6:279386/149‑1 (MQ=255)
                                                                                                         gATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCt                                   <  4:345356/148‑1 (MQ=255)
                                                                                                          aTAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTg                                 <  4:152055/149‑1 (MQ=255)
                                                                                                           tAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTgg                                <  1:357890/149‑1 (MQ=255)
                                                                                                           tAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTgg                                <  4:17599/149‑1 (MQ=255)
                                                                                                            aGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGa                                                                                      >  6:367410/1‑94 (MQ=255)
                                                                                                              cgcgcAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTg                             <  4:345901/149‑1 (MQ=255)
                                                                                                                         ccGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGaa                  >  4:9166/1‑149 (MQ=255)
                                                                                                                           gTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAAcc                <  2:370782/149‑1 (MQ=255)
                                                                                                                              tCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTCGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGc              >  8:61381/1‑148 (MQ=255)
                                                                                                                                    ccgagGCTGAGCAGCTTGCTCAGGCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTACCg         >  2:101421/3‑147 (MQ=255)
                                                                                                                                    ccgagGCTGAGCAGCTTGCTCAGGCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTACCg         >  4:211668/3‑147 (MQ=255)
                                                                                                                                    atgagGCTGAGCAGCTTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTACCGc        <  3:19527/146‑1 (MQ=255)
                                                                                                                                         ccTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTACCGCAGAGCg  >  5:131335/1‑149 (MQ=255)
                                                                                                                                           tGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGt                      <  7:65284/127‑1 (MQ=255)
                                                                                                                                           tGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGt                      >  8:65284/1‑127 (MQ=255)
                                                                                                                                           tGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTAc           >  3:363918/1‑138 (MQ=255)
                                                                                                                                           tGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTAc           <  4:363918/138‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GCACGGGTCGTGCTGGCCAATCCACATGGCTTGCATGCGCGACCGGCCAAGGTGCTTGCACAACTGGCCAAAGGCTTTGAAGGGGAAATCCGTGTGCGCCTGGTCGATAGCGCGCAGCCAGCCGTATCGGTGAAGAGCCTGAGCAAGCTGCTCAGCCTCGGCGCACGCCGTGGTCAGGCACTGGAACTGGTCGCCGAACCGAGCATCGCTGCCGATGCCTTGCCGGTGCTGCTGGCTGCTATCGAACAAGGCTTGGGTGAAGAGGTGGAACCGCTACCGCAGAGCG  >  NC_002947/910345‑910630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑