Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,528,844 G→T 5.3% L381L (CTG→CTT ftsA → cell division protein FtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,528,8440GT 5.3% 781.8 / 23.5 226L381L (CTG→CTTftsAcell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base G (97/117);  new base T (9/3);  total (106/120)
Fisher's exact test for biased strand distribution p-value = 7.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 235 total aligned reads displayed.

GCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGGCGCAAC  >  NC_002947/1528695‑1528992
                                                                                                                                                     |                                                                                                                                                     
gcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgctg                                                                                                                                                       >  2:1178360/1‑150 (MQ=255)
gcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgctg                                                                                                                                                       >  2:3853684/1‑150 (MQ=255)
gcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgctg                                                                                                                                                       >  1:2334664/1‑150 (MQ=255)
 cgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGa                                                                                                                                                      >  2:4203702/1‑150 (MQ=255)
  gAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                                                                                                                     <  2:3985713/150‑1 (MQ=255)
  gAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                                                                                                                     <  2:2957900/150‑1 (MQ=255)
     gATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCa                                                                                                                                                  >  2:4015120/1‑150 (MQ=255)
     gATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCa                                                                                                                                                  <  1:2639961/150‑1 (MQ=255)
       tGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGt                                                                                                                                                >  2:1372718/1‑150 (MQ=255)
          aaGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCg                                                                                                                                             <  2:1832572/150‑1 (MQ=255)
           aGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                                                                                                                           <  2:2887331/150‑1 (MQ=255)
           aGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                                                                                                                           <  2:1685397/150‑1 (MQ=255)
           aGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                                                                                                                           <  1:3096565/150‑1 (MQ=255)
           aGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGa                                                                                                                                           <  1:1392929/150‑1 (MQ=255)
            ggCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGACTGCTTCAGCAGGCCgta                                                                                                                                           <  2:2210309/150‑3 (MQ=255)
                ccGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACt                                                                                                                                      <  1:3745689/150‑1 (MQ=255)
                ccGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACt                                                                                                                                      <  1:625185/150‑1 (MQ=255)
                   tCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTgg                                                                                                                                   >  1:1545011/1‑150 (MQ=255)
                   tCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTgg                                                                                                                                   <  2:3647884/150‑1 (MQ=255)
                    cGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGc                                                                                                                                  <  1:2123721/150‑1 (MQ=255)
                     gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                                                                                                                                 <  1:1029859/150‑1 (MQ=255)
                     gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                                                                                                                                 <  2:1545685/150‑1 (MQ=255)
                     gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGcc                                                                                                                                 >  1:3827694/1‑150 (MQ=255)
                       aCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCt                                                                                                                               >  1:4310239/1‑150 (MQ=255)
                          ggCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGAc                                                                                                                            <  2:2402913/150‑1 (MQ=255)
                           gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGAcc                                                                                                                           <  1:3347113/150‑1 (MQ=255)
                             cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                                                                                                                         >  2:789114/1‑150 (MQ=255)
                             cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                                                                                                                         >  2:3256030/1‑150 (MQ=255)
                              gAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGt                                                                                                                        >  2:3520564/1‑150 (MQ=255)
                                 gAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaac                                                                                                                     <  2:2622995/150‑1 (MQ=255)
                                  aaaTCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaaca                                                                                                                    <  1:2167993/150‑1 (MQ=255)
                                   aaTCTTCCACATGCCGGTACGCATGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAg                                                                                                                                          <  1:3001188/127‑1 (MQ=255)
                                      cTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaa                                                                                                                <  2:653145/150‑1 (MQ=255)
                                       ttCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac                                                                                                               >  2:518795/1‑150 (MQ=255)
                                       ttCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac                                                                                                               <  1:4349798/150‑1 (MQ=255)
                                          cacaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTcacagcaacagc                                                                                                            >  1:673371/1‑150 (MQ=255)
                                                ccGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACagcagcagc                                                                                                      <  2:1017264/150‑1 (MQ=255)
                                                         cTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACATCTAt                                                                                             <  1:2623916/150‑1 (MQ=255)
                                                           gggCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATgg                                                                                           <  2:4478021/150‑1 (MQ=255)
                                                             gCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCg                                                                                         >  1:4056757/1‑150 (MQ=255)
                                                                      acaGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGa                                                                                >  1:169464/1‑150 (MQ=255)
                                                                       caGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGaa                                                                               <  2:2000653/150‑1 (MQ=255)
                                                                       caGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGaa                                                                               <  1:3647699/150‑1 (MQ=255)
                                                                              cGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCa                                                                        <  1:433564/150‑1 (MQ=255)
                                                                              cGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCa                                                                        <  1:3676840/150‑1 (MQ=255)
                                                                              cGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCa                                                                        <  1:101379/150‑1 (MQ=255)
                                                                                     tGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTg                                                                 >  1:1004584/1‑150 (MQ=255)
                                                                                       aGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGaa                                                               <  2:830730/150‑1 (MQ=255)
                                                                                             gtggtgCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTc                                                         <  1:2735129/150‑1 (MQ=255)
                                                                                             gtggtgCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTc                                                         <  2:1763827/150‑1 (MQ=255)
                                                                                               ggtgCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCaa                                                       >  2:1942442/1‑150 (MQ=255)
                                                                                                gtgCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCaag                                                      <  2:480006/150‑1 (MQ=255)
                                                                                                gtgCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCaag                                                      <  1:2015966/150‑1 (MQ=255)
                                                                                                    gcAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGACCGCTTCAAGAAGt                                                  <  1:4466686/150‑1 (MQ=255)
                                                                                                    gcAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGt                                                  <  1:3847052/150‑1 (MQ=255)
                                                                                                    gcAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGt                                                  >  2:2211891/1‑150 (MQ=255)
                                                                                                      aaCCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTgg                                                <  1:256305/150‑1 (MQ=255)
                                                                                                          cGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTcc                                            >  2:1625677/1‑150 (MQ=255)
                                                                                                          cGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTcc                                            >  2:4277739/1‑150 (MQ=255)
                                                                                                          cGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTcc                                            >  2:4045668/1‑150 (MQ=255)
                                                                                                           gATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCa                                           <  2:3156695/150‑1 (MQ=255)
                                                                                                             tCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAgg                                         <  2:3714884/150‑1 (MQ=255)
                                                                                                             tCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAgg                                         <  2:2484616/150‑1 (MQ=255)
                                                                                                              cTATTCCACTGGTGTGGGCGTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAggg                                        <  2:2862199/150‑1 (MQ=255)
                                                                                                                aTTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCa                                      >  2:1093799/1‑150 (MQ=255)
                                                                                                                 ttCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCaa                                     <  2:2578914/150‑1 (MQ=255)
                                                                                                                    cACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTTAACGCTTCAAGAAGTGGGTCCAGGGCAACtt                                  >  2:147849/1‑150 (MQ=255)
                                                                                                                    cACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACtt                                  >  2:2255567/1‑150 (MQ=255)
                                                                                                                     aCTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTc                                 <  2:1980548/150‑1 (MQ=255)
                                                                                                                     aCTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAA‑‑‑CAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTaa                              >  2:1266686/1‑150 (MQ=255)
                                                                                                                       tGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTa                               <  2:3646403/150‑1 (MQ=255)
                                                                                                                        ggTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTaa                              <  2:2767825/150‑1 (MQ=255)
                                                                                                                        ggTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTaa                              <  1:1943806/150‑1 (MQ=255)
                                                                                                                        ggTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTaa                              <  1:3155700/150‑1 (MQ=255)
                                                                                                                         gtgtGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAg                             >  1:4482894/1‑150 (MQ=255)
                                                                                                                         gtgtGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAg                             >  1:3295662/1‑150 (MQ=255)
                                                                                                                          tgtgGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGt                            <  2:3819454/150‑1 (MQ=255)
                                                                                                                                cTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCaaa                      <  2:1012896/150‑1 (MQ=255)
                                                                                                                                cTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCaaa                      <  1:791579/150‑1 (MQ=255)
                                                                                                                                   gCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCa                   <  1:492202/150‑1 (MQ=255)
                                                                                                                                     tCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagt                 >  2:4417770/1‑150 (MQ=255)
                                                                                                                                     tCACCTACGGCCTGATGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagt                 >  1:1612719/1‑150 (MQ=255)
                                                                                                                                      cACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagta                <  2:2688037/150‑1 (MQ=255)
                                                                                                                                      cACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagta                <  2:2326358/150‑1 (MQ=255)
                                                                                                                                      cACCTACGACCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagta                >  1:3942840/1‑150 (MQ=255)
                                                                                                                                       aCCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagtag               >  1:2914705/1‑150 (MQ=255)
                                                                                                                                         cTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTg             >  2:3457476/1‑150 (MQ=255)
                                                                                                                                          ctcGGACTGCTTCAGCAGGCCGTAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTg             >  1:3412929/3‑150 (MQ=255)
                                                                                                                                           tcGGACTGCTTCAGCAGGCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGc                                                            >  1:4018083/2‑101 (MQ=255)
                                                                                                                                           tcGGACTGCTTCAGCAGGCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGt           >  1:2262201/2‑150 (MQ=255)
                                                                                                                                           tcGGACTGCTTCAGCAGGCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGt           >  1:877002/2‑150 (MQ=255)
                                                                                                                                           tcGGACTGCTTCAGCAGGCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGt           >  1:876931/2‑150 (MQ=255)
                                                                                                                                           tcGGACTGCTTCAGCAGGCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGt           >  2:3086995/2‑150 (MQ=255)
                                                                                                                                           aCGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGt           <  1:2725135/150‑1 (MQ=255)
                                                                                                                                              gCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCagtagt              >  2:535871/1‑144 (MQ=255)
                                                                                                                                               ccTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGgc       <  2:429616/150‑1 (MQ=255)
                                                                                                                                               ccTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGgc       <  2:2130773/150‑1 (MQ=255)
                                                                                                                                                   ctgAAGCAGTCCGAGGACTTGGCCCTGACTGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGGCGcaa   <  1:1178360/150‑1 (MQ=255)
                                                                                                                                                   ctgAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGGCGcaa   <  2:1726106/150‑1 (MQ=255)
                                                                                                                                                    tgAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGGCGcaac  >  1:1071860/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGATGAACCGAAGGCCCCAGTGCTTGAACGCTTCAAGAAGTGGGTCCAGGGCAACTTCTAAGTTTCAAAGCAGTAGTGGTAGGCGCAAC  >  NC_002947/1528695‑1528992

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: