Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,282,185 T→C 7.0% T349A (ACC→GCC) ‡ 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,1850TC 7.0% 202.5 / 5.2 71T349A (ACC→GCC) ‡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 = 6.58e-01

GGCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCG  >  NC_002947/4282038‑4282330
                                                                                                                                                   |                                                                                                                                                 
ggCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCggtggt                                                                                                                                                   >  6:143062/1‑148 (MQ=255)
ggCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCggtggt                                                                                                                                                   <  5:143062/148‑1 (MQ=255)
 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                                                             >  7:55969‑M2/2‑46 (MQ=255)
                                                                                                         gcccgCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCttcgcgggcatgcccgctcccacacggacctcacacaagccagatcagtggtggtgcttgccactcagggccgcatggaaatt                                                             <  8:55969‑M2/128‑84 (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                                               >  7:78557/1‑131 (MQ=255)
                                                                                                                     cGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCgg                                               <  8:78557/131‑1 (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)
                                                                                                                                                   |                                                                                                                                                 
GGCTTTGCCATGGTGCCGGCCTGATCCTTTCGGGTGGCCTGTGCTGGCCTATTCGCGGGCAAGCCCGCTCCTACAATTCTGGGCGTGCCCTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGGTTTCTTCGCGGGCATGCCCGCTCCTACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCG  >  NC_002947/4282038‑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: