Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,282,186 T→C 7.0% E348E (GAA→GAG) ‡ PP_3753 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,282,1860TC 7.0% 206.7 / 5.0 71E348E (GAA→GAG) ‡PP_3753AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (2/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCG  >  NC_002947/4282039‑4282330
                                                                                                                                                   |                                                                                                                                                
gCTTTGCCATGGTGCCGGCCTGATCCTTTCTGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACAGGCAGAGCCGGTGCCATCCACCGCGGTGGttt                                                                                                                                                 >  6:245742/1‑149 (MQ=255)
gCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGttt                                                                                                                                                 >  3:87706/1‑149 (MQ=255)
gCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGttt                                                                                                                                                 >  3:83313/1‑149 (MQ=255)
 cTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTttc                                                                                                                                                >  2:230549/1‑149 (MQ=255)
    tGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCt                                                                                                                                               >  5:10048‑M2/1‑146 (MQ=255)
       cATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcg                                                                                                                                          >  4:191821‑M2/1‑143 (MQ=255)
          tgTGCCGGCCTGATCCTTTGGGGTGTCCTGTGCTGGCCTATTCGCGTGCAAGCCCTATCCTACAATTCTGTGCGTGCCCTGTATGAGCGGGTTTAGCCGCGAAAACCGGCAGAGCCGGTGCCATCCAGGGAGGTGGTTTGTTCGCggg                                                                                                                                        <  8:176181/147‑1 (MQ=25)
            tGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcgggcat                                                                                                                                     >  3:87298‑M2/1‑138 (MQ=255)
              ccGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcgggcatgc                                                                                                                                   <  6:305916‑M2/149‑14 (MQ=255)
                   ctcaTCCTTTCGGGTGGCCTGTGCTGGCCTATTAGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGCAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCCGAGGTTTCTTCGCGGGCATGCCCGCt                                                                                                                              >  1:122113/4‑149 (MQ=37)
                    tGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcgggcatgcccgctc                                                                                                                             <  7:126672‑M2/149‑20 (MQ=255)
                    tGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcgggcatgcccgct                                                                                                                              <  2:132730‑M2/148‑19 (MQ=255)
                      aTCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcgggcatgcccgctc                                                                                                                             >  1:87431‑M2/1‑128 (MQ=255)
                      aTCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgcgggcatgcccgct                                                                                                                              >  1:132730‑M2/1‑128 (MQ=255)
                             cGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGt                                                                                                                     <  4:87706/148‑1 (MQ=255)
                             cGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGt                                                                                                                     <  4:83313/148‑1 (MQ=255)
                              gggTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTa                                                                                                                    <  8:126672/148‑1 (MQ=255)
                                gTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGaa                                                                                                                 <  3:135665/149‑1 (MQ=255)
                                  ggCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACg                                                                                                               >  3:232996/1‑149 (MQ=255)
                                    cgtgtgCTGGCCTATTCGCGGGCAAGCGCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCa                                                                                                             <  8:274374/147‑1 (MQ=255)
                                           tGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACa                                                                                                                        <  5:89393/131‑1 (MQ=255)
                                           tGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACa                                                                                                                        >  6:89393/1‑131 (MQ=255)
                                               cTATTCGTGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGc                                                                                                  >  4:262564/1‑149 (MQ=255)
                                                 aggagcgGGCTTGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCcg                                                                                                >  3:170785/5‑149 (MQ=255)
                                                 aTTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCcg                                                                                                >  7:61647/1‑149 (MQ=255)
                                                  ttCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGcc                                                                                                 <  1:230549/147‑1 (MQ=255)
                                                     gcgGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCTGTGCCATCCACCGCGGTTGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCcgc                                                                                               <  5:289072/146‑1 (MQ=255)
                                                     gcgGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCgcgc                                                                                             <  4:59628/148‑1 (MQ=255)
                                                     gcgGGCAAGCCCGCTCCTACAATTCTGGGAGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGAAGAGCCGGTGCCATCCACCGCCGTGGTTTCTTCGCGGGCATGCACGCTCATACAGGTAGAACGCAGTCCCACCAGCcgc                                                                                               >  6:289072/1‑146 (MQ=25)
                                                      cgGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGccc                                                                                           >  8:197171/1‑149 (MQ=255)
                                                              cccGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttc                                                                                                                                             >  4:274943‑M2/1‑88 (MQ=255)
                                                              cccGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttc                                                                                                                                             <  3:274943‑M2/91‑4 (MQ=255)
                                                              cccGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTc                                                                                   >  6:284950/1‑149 (MQ=255)
                                                               ccGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTcc                                                                                  >  3:56868/1‑149 (MQ=255)
                                                                   tCCTACAATTCTGGGCGTACCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTgg                                                                              >  4:258628/1‑149 (MQ=255)
                                                                     cTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCt                                                                                >  7:296744/1‑145 (MQ=255)
                                                                     cTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCt                                                                                <  8:296744/145‑1 (MQ=255)
                                                                     cTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGc                                                                            >  6:310604/1‑149 (MQ=255)
                                                                     cTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGc                                                                            >  2:169348/1‑149 (MQ=255)
                                                                         aaTTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgc                                                                                                                                           >  6:134997‑M2/1‑77 (MQ=255)
                                                                         aaTTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCttcgc                                                                                                                                           <  5:134997‑M2/82‑6 (MQ=255)
                                                                         aaTTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTg                                                                        >  7:159002/1‑149 (MQ=255)
                                                                           ttCTGGGCGTACCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGcc                                                                      <  3:258628/149‑1 (MQ=255)
                                                                              tGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCgg                                                                   >  1:190725/1‑149 (MQ=255)
                                                                                       ccTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCActtct                                                          <  4:170785/149‑1 (MQ=255)
                                                                                        cTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCActtctt                                                         <  4:56868/149‑1 (MQ=255)
                                                                                                  gggTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTGGAACGCCCgg                                               <  5:284950/149‑1 (MQ=255)
                                                                                                       tACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGa                                           <  7:171850/148‑1 (MQ=255)
                                                                                                        gcccgCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCttcgcgggcatgcccgctcccacacggacctcacacaagccagatcagtggtggtgcttgccactcagggccgcatggaaattggcgctctggcgcaggtact                                         <  2:61720‑M2/148‑104 (MQ=255)
                                                                                                        gcccgCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCttcgcgggcatgcccgctcccacacggacctcacacaagccagatcagtggtggtgcttgccactcagggccgcatggaaatt                                                             <  8:55969‑M2/128‑84 (MQ=255)
                                                                                                        gcccgCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCttcgcgggcatgcccgctcccacacggacctcacacaagccagatcagtggtggtgcttgccactcagggccgcatggaaatt                                                             >  7:55969‑M2/2‑46 (MQ=255)
                                                                                                          ccGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCttcgcgggcatgcccgctcccacaaggacctcacacaagccagatcagtggtggtgcttgccactcagggccgcatggaaattggcgctctggcgcaggtactgc                                       <  2:46414‑M2/149‑106 (MQ=255)
                                                                                                             cgAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCttcgcgggcatgcccgctcccacaaggacctcacacaagccagatcagtggtggtgcttgccactcagggccgcatggaaattggcgctctggcgcaggtactgctcg                                    >  1:46414‑M2/1‑41 (MQ=255)
                                                                                                                    gggCAGAGCCGGTGCCATCCACCGCGGGGGTTTCTTCGCCGGCAGGCCCTCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTcc                             <  6:100019/148‑1 (MQ=255)
                                                                                                                    cGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCgg                                               <  8:78557/131‑1 (MQ=255)
                                                                                                                    cGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCgg                                               >  7:78557/1‑131 (MQ=255)
                                                                                                                     ggCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCt                            <  2:190725/149‑1 (MQ=255)
                                                                                                                      gCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCCCTCGGTTCCt                            >  1:139319/1‑148 (MQ=255)
                                                                                                                       cAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGa                          <  8:61647/149‑1 (MQ=255)
                                                                                                                         gagCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGaaa                        <  5:310604/149‑1 (MQ=255)
                                                                                                                         gagCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGaaa                        <  8:159002/149‑1 (MQ=255)
                                                                                                                             cGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAAcccc                    <  7:295204/149‑1 (MQ=255)
                                                                                                                               gTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGc                  >  4:17875/1‑149 (MQ=255)
                                                                                                                                tGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGcc                 <  4:148355/149‑1 (MQ=255)
                                                                                                                                   cATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGc              >  6:169871/1‑149 (MQ=255)
                                                                                                                                    aTCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTg                                                        <  8:127603/106‑1 (MQ=255)
                                                                                                                                    aTCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTg                                                        >  7:127603/1‑106 (MQ=255)
                                                                                                                                     tCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGATGGAAGGCACTCGGTTCCTGAAACCCCg                   >  2:34040/1‑142 (MQ=255)
                                                                                                                                     tCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGATGGAAGGCACTCGGTTCCTGAAACCCCg                   <  1:34040/142‑1 (MQ=255)
                                                                                                                                       cACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCgg          <  4:6419/149‑1 (MQ=255)
                                                                                                                                       cACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCAGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCgg          >  4:112238/1‑149 (MQ=255)
                                                                                                                                        aCCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCgg          <  7:80941/148‑1 (MQ=255)
                                                                                                                                        aCCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCgg          >  4:299461/1‑148 (MQ=255)
                                                                                                                                        aCCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCgg          >  8:80941/1‑148 (MQ=255)
                                                                                                                                               tggtTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCg  >  8:238978/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCG  >  NC_002947/4282039‑4282330

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: