Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,998,303 A→C 8.8% E330A (GAG→GCG)  PP_3526 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,998,3030AC 8.8% 240.9 / 13.9 80E330A (GAG→GCG) PP_3526AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base C (4/3);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

ACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAGGTC  >  NC_002947/3998156‑3998450
                                                                                                                                                    |                                                                                                                                                   
aCATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCTCCCTGTCGGTGATGGAAGTGGCCGGGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGa                                                                                                                                                     <  1:116820/148‑1 (MQ=255)
   tGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCaa                                                                                                                                                 >  2:60530/1‑149 (MQ=255)
     gAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAAc                                                                                                                                                <  5:337562/148‑1 (MQ=255)
     gAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGGGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAAc                                                                                                                                                <  1:209556/148‑1 (MQ=255)
      aGCTGCGGATAACCAAGAGCCGCCGGCTGCTGGAGTATTCCGTCCTGTCGGTGATGGAACTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTGGGCTACCCGCCGTCGCGCGAGCAACgg                                                                                                                                              <  4:70975/149‑1 (MQ=255)
       gCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGc                                                                                                                                             <  8:296625/149‑1 (MQ=255)
           cgcATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTg                                                                                                                                           <  1:268404/147‑1 (MQ=255)
            gcATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGc                                                                                                                                        <  3:82860/149‑1 (MQ=255)
                  ccAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGa                                                                                                                                  >  7:30565/1‑149 (MQ=255)
                      gagCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGa                                                                                                                              <  5:250636/149‑1 (MQ=255)
                            ccgGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCt                                                                                                                        >  2:363223/1‑149 (MQ=255)
                              ggCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCttt                                                                                                                      <  8:30565/149‑1 (MQ=255)
                                ctgctgCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGa                                                                                                                    <  1:60530/149‑1 (MQ=255)
                                 tgctgcAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGag                                                                                                                   <  1:132091/149‑1 (MQ=255)
                                  gctgcAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagc                                                                                                                  >  5:366592/1‑149 (MQ=255)
                                     gcAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccag                                                                                                               >  8:130983/1‑149 (MQ=255)
                                      cAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagtcagc                                                                                                              >  7:139005/1‑149 (MQ=255)
                                      cAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccagc                                                                                                              >  5:40201/1‑149 (MQ=255)
                                      cAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagcca                                                                                                                <  2:168907/147‑1 (MQ=255)
                                       aGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccagcc                                                                                                             >  8:162652/1‑149 (MQ=255)
                                       aGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccagcc                                                                                                             <  5:294017/149‑1 (MQ=255)
                                           ttCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagcca                                                                                                                >  5:63129/1‑142 (MQ=255)
                                           ttCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagcca                                                                                                                <  6:63129/142‑1 (MQ=255)
                                             ccGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGgcg                                                                                                       >  3:188377/1‑149 (MQ=255)
                                                aCCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTc                                                                                                    <  1:363223/149‑1 (MQ=255)
                                                  cTGTCGGTTATGGAAGTGGCGGTGGCCTGGGGGTTTGGCTCGGTTGCGCACTTCATCAATTGCTATTCGGCGTATTTCGGCTACCTCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTg                                                                                                                     <  1:256275/131‑1 (MQ=255)
                                                  cTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTg                                                                                                                     >  2:256275/1‑130 (MQ=255)
                                                  cTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGa                                                                                                  <  8:63334/149‑1 (MQ=255)
                                                    gTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCAGTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccag                                                                                                               <  8:70618/134‑1 (MQ=255)
                                                    gTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGa                                                                                                                                                     >  7:234039/1‑96 (MQ=255)
                                                    gTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccag                                                                                                               >  7:70618/1‑134 (MQ=255)
                                                    gTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGGCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGa                                                                                                                                                     <  8:234039/96‑1 (MQ=255)
                                                         tGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAAtct                                                                                           >  2:318429/1‑149 (MQ=255)
                                                          gATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAAtct                                                                                           >  8:280389/1‑148 (MQ=255)
                                                             ggAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGAATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTCGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAgg                                                                                       >  4:295058/1‑149 (MQ=255)
                                                                 gTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACg                                                                                   <  3:358931/149‑1 (MQ=255)
                                                                  tGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGa                                                                                  <  7:130983/149‑1 (MQ=255)
                                                                          gCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGa                                                                           <  7:162652/148‑1 (MQ=255)
                                                                            cTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCa                                                                                                                                                  <  4:158890/75‑1 (MQ=255)
                                                                            cTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCa                                                                                                                                                  >  3:158890/1‑75 (MQ=255)
                                                                               cGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTgg                                                                      <  6:40201/148‑1 (MQ=255)
                                                                               cGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGa                                                                     >  1:180320/1‑149 (MQ=255)
                                                                                gggTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGAt                                                                    <  2:180320/149‑1 (MQ=255)
                                                                                 ggTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCt                                                                        >  4:259779/1‑144 (MQ=255)
                                                                                 ggTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCt                                                                        <  3:259779/144‑1 (MQ=255)
                                                                                       tctcGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCAGCCGTTGCTCGCGCGACGGCgggt                                                                                                                                         <  2:199072/73‑2 (MQ=25)
                                                                                       tctcGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCAGCCGTTGCTCGCGCGACGGCgggt                                                                                                                                         >  1:199072/1‑72 (MQ=25)
                                                                                                 gCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCAt                                                   >  2:2274/1‑149 (MQ=255)
                                                                                                 gCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATATCAGGTACGAGGCTAGAGCTGGATGTACCGCAACAACGAAt                                                   >  2:117414/1‑149 (MQ=255)
                                                                                                   aCTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATca                                                 <  6:366592/149‑1 (MQ=255)
                                                                                                   aCTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGAGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATca                                                 >  3:124454/1‑149 (MQ=255)
                                                                                                    cTTCAGCAAGTGCTATGCGGCGTATTTCGGCTAGCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTATGAGGCTCGAGCTGGATGTACCGCAACACCGCATcac                                                <  1:241861/149‑1 (MQ=255)
                                                                                                       cAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATcaca                                               <  4:188377/147‑1 (MQ=255)
                                                                                                       cAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATca                                                 >  8:307644/1‑145 (MQ=255)
                                                                                                           aaGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAg                                                                                                          <  2:143202/84‑1 (MQ=255)
                                                                                                           aaGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAg                                                                                                          >  1:143202/1‑84 (MQ=255)
                                                                                                           aaGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACAcc                                                       >  8:129272/1‑135 (MQ=255)
                                                                                                           aaGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACAcc                                                       <  7:129272/135‑1 (MQ=255)
                                                                                                                cTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGAGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCt                                     <  4:124454/148‑1 (MQ=255)
                                                                                                                 tATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCttt                                   >  1:291713/1‑149 (MQ=255)
                                                                                                                       gcgTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCAGGAAGTACCGCAACACCGCAACACAGATCATGGCATtgatgcg                             >  8:141451/1‑144 (MQ=255)
                                                                                                                         gTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGgcg                             >  7:67342/1‑147 (MQ=255)
                                                                                                                            tttCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCt                        >  5:116052/1‑149 (MQ=255)
                                                                                                                             ttCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCtt                       >  5:91939/1‑149 (MQ=255)
                                                                                                                              tCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTgg                                                                      >  8:330122/1‑101 (MQ=255)
                                                                                                                              tCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTgg                                                                      <  7:330122/101‑1 (MQ=255)
                                                                                                                                ggCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAt                    >  6:311262/1‑149 (MQ=255)
                                                                                                                                 gCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAtt                   <  1:318429/149‑1 (MQ=255)
                                                                                                                                  cTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAttt                  <  6:283421/149‑1 (MQ=255)
                                                                                                                                     cccGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGTCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGc                <  8:139005/148‑1 (MQ=255)
                                                                                                                                     cccGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTg                 <  5:265409/147‑1 (MQ=255)
                                                                                                                                      ccagccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAttt                  <  2:129092/143‑1 (MQ=255)
                                                                                                                                      ccagccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAtt                   <  6:315695/142‑1 (MQ=255)
                                                                                                                                      ccagccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAtt                   >  5:315695/4‑145 (MQ=255)
                                                                                                                                      ccagccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTgaa               >  1:129092/4‑147 (MQ=255)
                                                                                                                                        cgccgTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtca             <  8:67342/149‑1 (MQ=255)
                                                                                                                                         gccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtca             >  7:196659/1‑148 (MQ=255)
                                                                                                                                          ccgTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtcatc           <  7:280389/149‑1 (MQ=255)
                                                                                                                                          ccgTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtcatc           <  1:2274/149‑1 (MQ=255)
                                                                                                                                             tCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGa        <  3:295058/149‑1 (MQ=255)
                                                                                                                                                 gcgAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAgg    >  6:268210/1‑149 (MQ=255)
                                                                                                                                                 gcgAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAgg    >  1:149162/1‑149 (MQ=255)
                                                                                                                                                   gAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAGGTc  >  7:185022/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAGGTC  >  NC_002947/3998156‑3998450

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: