Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,984,700 G→T 6.1% A140A (GCC→GCA ycaC‑II ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,984,7000GT 6.1% 262.5 / 4.5 82A140A (GCC→GCAycaC‑IIputative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (46/31);  new base T (3/2);  total (49/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAG  >  NC_002947/5984556‑5984847
                                                                                                                                                |                                                                                                                                                   
cagcagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGACCGATGcct                                                                                                                                                 <  5:379507/149‑3 (MQ=255)
   cagGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGTCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCgg                                                                                                                                              >  1:97775/1‑149 (MQ=255)
       ttACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcac                                                                                                                                          <  1:202743/149‑1 (MQ=255)
         aCCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcacca                                                                                                                                        >  5:168577/1‑149 (MQ=255)
          cccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcacca                                                                                                                                        <  1:12908/148‑1 (MQ=255)
           ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       >  7:65365/1‑148 (MQ=255)
           ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       >  5:122515/1‑148 (MQ=255)
           ccAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTcaccac                                                                                                                                       >  8:77567/1‑148 (MQ=255)
             aGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGCACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGaa                                                                                                                                    >  5:269132/1‑149 (MQ=255)
             aGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGaa                                                                                                                                    >  4:269372/1‑149 (MQ=255)
             aGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGACCGATGCCTCGGGcaccttcaa                                                                                                                                    <  8:73422/149‑6 (MQ=37)
                       tGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCaaaa                                                                                                                          <  6:168577/149‑1 (MQ=255)
                         tCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAAcc                                                                                                                        >  3:129457/1‑149 (MQ=255)
                           gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCg                                                                                                                                               <  8:276282/124‑1 (MQ=255)
                           gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCg                                                                                                                                               >  7:276282/1‑124 (MQ=255)
                           gTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCt                                                                                                                       <  5:322356/148‑1 (MQ=255)
                                gcCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCGTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCAGCGGTCACCACGAACACTTCAACACCTTCCGc                                                                                                                  >  2:98363/1‑148 (MQ=255)
                                   aGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAgc                                                                                                              >  3:109548/1‑149 (MQ=255)
                                     tCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAgcgc                                                                                                            >  2:206086/1‑149 (MQ=255)
                                     tCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAgcgc                                                                                                            >  2:249484/1‑149 (MQ=255)
                                        cgGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGa                                                                                                         <  2:97775/149‑1 (MQ=255)
                                             cAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCACAACCTTCCGCCAGCGCCGACAggg                                                                                                    >  7:205094/1‑149 (MQ=255)
                                             cAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAgg                                                                                                     <  3:269372/148‑1 (MQ=255)
                                               gTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTc                                                                                                  >  5:64801/1‑149 (MQ=255)
                                                  cGCAGGTCACCGAGAACCATTTCATCATCTGCGCACCGGCGTGGGTCATGCGTACCCAAGCCGGTTGCTGCACGGTTTCGTTGAAGGTGCCCGTGGCATCGGTAAGCAGGAACTCTTCAAAACGTTCCGCCAGCGCCGACAGGGTCgg                                                                                                <  1:69412/148‑1 (MQ=255)
                                                   gCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCggg                                                                                               <  8:158965/148‑1 (MQ=255)
                                                     aGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAg                                                                                            <  6:64801/149‑1 (MQ=255)
                                                       gCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGCGAAGGc                                                                                          >  6:244292/1‑149 (MQ=255)
                                                        ccACCGAGAACCAGTTCATCATCTGCGCACCCGGCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGcc                                                                                         <  5:199470/149‑1 (MQ=255)
                                                         cACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCa                                                                                        >  5:347161/1‑149 (MQ=255)
                                                          aCCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCggg                                                                                               >  4:193724/1‑141 (MQ=255)
                                                          aCCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCggg                                                                                               <  3:193724/141‑1 (MQ=255)
                                                               gaACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGATTGCTGAACGGTTTCGTTGAAGTTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGcaca                                                                                  >  5:381028/1‑149 (MQ=255)
                                                                  ccAGTTCATCATCTACGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCACACCCTTCCGCCAGCGCCGACATGGTCGGGAAGGCCACGccaccat                                                                               >  4:15069/1‑143 (MQ=38)
                                                                   cAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATc                                                                              >  6:273767/1‑149 (MQ=255)
                                                                   cAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATc                                                                              >  2:357050/1‑149 (MQ=255)
                                                                         atcatcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAACACCTTCCGCCAGCGCCGACAGGGTCGGGAACGCCACGCACACATCGGTcac                                                                        >  5:365542/1‑149 (MQ=255)
                                                                         atcatcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGGTTCGTTGAAGGTGCGCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTcac                                                                        <  6:240557/149‑1 (MQ=255)
                                                                            atcTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTcacca                                                                      >  5:271227/1‑148 (MQ=255)
                                                                                gcgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAgg                                                                                           <  1:163939/123‑1 (MQ=255)
                                                                                gcgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAgg                                                                                           >  2:163939/1‑123 (MQ=255)
                                                                                gcgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCg                                                                 >  5:204949/1‑149 (MQ=255)
                                                                                 cgcACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACAcc                                                                  <  6:347161/147‑1 (MQ=255)
                                                                                      ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCgat                                                            >  3:137850/1‑148 (MQ=255)
                                                                                      ccGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCgat                                                            >  1:318159/1‑148 (MQ=255)
                                                                                         ccTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATc                                                        <  4:109548/149‑1 (MQ=255)
                                                                                             ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCCGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGtt                                                    <  6:269132/149‑1 (MQ=255)
                                                                                             ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAg                                                      <  7:259639/147‑1 (MQ=255)
                                                                                             ggTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGtt                                                    <  6:122515/149‑1 (MQ=255)
                                                                                                    cGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCg                                             <  7:77567/149‑1 (MQ=255)
                                                                                                     gTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCgg                                            >  1:192531/1‑149 (MQ=255)
                                                                                                        cccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGTGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGccc                                         <  5:244292/149‑1 (MQ=255)
                                                                                                         ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCa                                                       >  7:103128/1‑134 (MQ=255)
                                                                                                         ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCa                                                       <  8:103128/134‑1 (MQ=255)
                                                                                                         ccAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCg                                        >  3:14412/1‑149 (MQ=255)
                                                                                                               cGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCt                                  >  5:121174/1‑149 (MQ=255)
                                                                                                                  ttGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGa                               >  7:111333/1‑149 (MQ=255)
                                                                                                                    gctgcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGa                               <  6:271227/147‑1 (MQ=255)
                                                                                                                     ctgcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGAtt                             <  2:119596/148‑1 (MQ=255)
                                                                                                                     ctgcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGAtt                             <  6:204949/148‑1 (MQ=255)
                                                                                                                       gcACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGcc                          >  8:96079/1‑149 (MQ=255)
                                                                                                                        cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCCCGCACACATCGGTCACCACCCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                         >  6:42923/1‑149 (MQ=255)
                                                                                                                        cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                         <  5:273767/149‑1 (MQ=255)
                                                                                                                        cACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCt                         >  7:32838/1‑149 (MQ=255)
                                                                                                                               tCGTTGAAGGTGCCAGAGGAATCTGTCACCACGAAAACTTCAAAACCTTCCGCCAGCGCCTACAGGGTCGGGAAGGCCACTCAGACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGa                   <  6:381028/148‑1 (MQ=255)
                                                                                                                                cGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAg                 >  8:101734/1‑149 (MQ=255)
                                                                                                                                 gTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGt                >  3:310275/1‑149 (MQ=255)
                                                                                                                                  ttGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTc               >  6:278779/1‑149 (MQ=255)
                                                                                                                                    gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTg                       <  1:79351/139‑1 (MQ=255)
                                                                                                                                    gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTg                       >  2:79351/1‑139 (MQ=255)
                                                                                                                                      tggTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCCCCACACGGGCGATGCTCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg           >  2:296962/2‑149 (MQ=255)
                                                                                                                                       ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCGTTGACGAAGTCTTCGt          >  6:86419/1‑149 (MQ=255)
                                                                                                                                       ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGt          >  7:355137/1‑149 (MQ=255)
                                                                                                                                       ggTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCAACACACCGGCGATGATCAGTTGCTTCCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGt          >  7:200541/1‑149 (MQ=25)
                                                                                                                                        gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt         >  1:148453/1‑149 (MQ=255)
                                                                                                                                        gTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt         >  6:137461/1‑149 (MQ=255)
                                                                                                                                          gCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGt       >  7:213292/1‑149 (MQ=255)
                                                                                                                                           cccGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTc      <  8:355137/149‑1 (MQ=255)
                                                                                                                                            ccGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTcc     >  3:151635/1‑149 (MQ=255)
                                                                                                                                              gAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCa   <  8:32838/149‑1 (MQ=255)
                                                                                                                                              gAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCa   <  2:192531/149‑1 (MQ=255)
                                                                                                                                               aGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAg  >  3:66278/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGTTCGCAGGCCACCGAGAACCAGTTCATCATCTGCGCACCCGCCTGGGTCATGCGTACCCAAGCCGCTTGCTGCACGGTTTCGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAG  >  NC_002947/5984556‑5984847

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: