Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,445,045 T→A 6.0% E156V (GAA→GTA) ‡ lexA‑I ← transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,445,0450TA 6.0% 207.5 / 3.1 67E156V (GAA→GTA) ‡lexA‑Itranscriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

TGGATCAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTGGGGGTGG  >  NC_002947/2444899‑2445185
                                                                                                                                                  |                                                                                                                                            
tGGATCAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCATCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCg                                                                                                                                            >  8:318688/1‑149 (MQ=255)
     cAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCAcc                                                                                                                                       >  5:66865/1‑149 (MQ=255)
     cAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACGAAGTCAcc                                                                                                                                       >  6:221614/1‑149 (MQ=255)
      aGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCg                                                                                                                                      <  3:51625/149‑1 (MQ=255)
      aGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCg                                                                                                                                      <  6:66865/149‑1 (MQ=255)
      aGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACATTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCAcc                                                                                                                                       <  8:35865/148‑1 (MQ=255)
       gcggcgAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGa                                                                                                                                     >  2:249868/1‑149 (MQ=255)
          gcgAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACGAAGTCACCgtcaa                                                                                                                                  <  5:221614/149‑5 (MQ=255)
          gcgAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACGAAGTCACCgtcaa                                                                                                                                  <  7:126702/149‑5 (MQ=255)
          gcccacgACGCCGCCGCTCACCCCCGCCACCCCCCGCACCTGTTCTTTCCCGCCGCCTTCGCAGGGGGCGAATTCGGGGGTTGCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATgc                                                                                                                                  <  8:168943/143‑1 (MQ=255)
                 aCGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCca                                                                                                                            >  1:128479/1‑148 (MQ=255)
                   gCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACaa                                                                                                                         <  3:250006/149‑1 (MQ=255)
                    ccGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAAc                                                                                                                        <  6:91790/149‑1 (MQ=255)
                    ccGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAAc                                                                                                                        >  1:333615/1‑149 (MQ=255)
                          cTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTg                                                                                                                  <  7:163931/149‑1 (MQ=255)
                           tCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTgg                                                                                                                 >  7:151673/1‑149 (MQ=255)
                             aaGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGcc                                                                                                               <  3:189555/149‑1 (MQ=255)
                                  cTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGc                                                                                                          >  8:208061/1‑149 (MQ=255)
                                         aCCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCtt                                                                                                    <  1:249868/148‑1 (MQ=255)
                                          ccAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCa                                                                                                  >  8:566/1‑149 (MQ=255)
                                           cAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGT‑‑CCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCa                                                                                                  <  3:162346/146‑1 (MQ=255)
                                           cAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGT‑‑CCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCa                                                                                                  >  4:162346/1‑146 (MQ=255)
                                                 cTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCa                                                                                             <  7:208061/147‑1 (MQ=255)
                                                 cTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAgg                                                                                           >  1:149945/1‑149 (MQ=255)
                                                 cTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAg                                                                                            <  1:78910/148‑1 (MQ=255)
                                                   gTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGgtg                                                                                         >  8:172317/1‑149 (MQ=255)
                                                    ttCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGgt                                                                                          <  8:83035/147‑1 (MQ=255)
                                                    ttCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGgt                                                                                          <  1:222200/147‑1 (MQ=255)
                                                       tttCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCa                                                                                     >  4:121552/1‑149 (MQ=255)
                                                        ttCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCa                                                                                     <  7:318688/148‑1 (MQ=255)
                                                         tCAGGTCGCCTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTATGCTGCCTTCGCGCTTGAAGCCCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGc                                                                                      <  6:121157/146‑1 (MQ=255)
                                                         tCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGc                                                                                      <  7:172317/146‑1 (MQ=255)
                                                         tCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCAcc                                                                                   >  4:25517/1‑149 (MQ=255)
                                                                gACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCcagca                                                                            <  2:333615/149‑1 (MQ=255)
                                                                  cTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAgg                                                                          <  2:128479/149‑1 (MQ=255)
                                                                       aTGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGcc                                                                     <  4:338900/149‑1 (MQ=255)
                                                                          ggggCGAATGCGGGGTTTTCGGCTAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTGCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGt                                                                   <  7:187047/148‑1 (MQ=255)
                                                                           gggCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTc                                                                  >  4:14801/1‑148 (MQ=255)
                                                                            ggCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCg                                                                 >  5:340715/1‑148 (MQ=255)
                                                                             gCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCg                                                                 >  8:28351/1‑147 (MQ=255)
                                                                             gCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCg                                                                 <  7:28351/147‑1 (MQ=255)
                                                                                 aTTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATg                                                           <  8:31603/149‑1 (MQ=255)
                                                                                     ggggTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGa                                                       >  1:126126/1‑149 (MQ=255)
                                                                                        gTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGt                                                    >  5:284598/1‑149 (MQ=255)
                                                                                           ttCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTcc                                                  <  8:151673/148‑1 (MQ=255)
                                                                                                 aaGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTtcatgc                                           <  3:14801/149‑1 (MQ=255)
                                                                                                 aaGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTtcatgc                                           <  3:25517/149‑1 (MQ=255)
                                                                                                    cagcCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTtcatgctcn                                        <  2:149945/149‑2 (MQ=255)
                                                                                                      gcCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTtcatgctcatg                                      >  4:223194/1‑149 (MQ=255)
                                                                                                           aCTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCgtg                                 >  2:178611/1‑149 (MQ=255)
                                                                                                                gctgcCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTAcgc                            >  5:337094/1‑149 (MQ=255)
                                                                                                                   gcCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCgt                                  <  6:94414/140‑1 (MQ=255)
                                                                                                                   gcCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCgt                                  >  5:94414/1‑140 (MQ=255)
                                                                                                                    ccTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCgca                         <  6:135862/148‑1 (MQ=255)
                                                                                                                     cTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCgcag                        >  8:135260/1‑148 (MQ=255)
                                                                                                                      ttCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCgcagc                       >  8:45122/1‑148 (MQ=255)
                                                                                                                         gcgcTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCaaa                    <  1:178611/148‑1 (MQ=255)
                                                                                                                          cgcTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATa                  <  7:135260/149‑1 (MQ=255)
                                                                                                                            cTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGt                >  8:116149/1‑149 (MQ=255)
                                                                                                                             ttGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCCCATGCCGTGTACGCGCAGCAAATAGTc               >  6:217616/1‑149 (MQ=255)
                                                                                                                               gAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTc               <  3:223194/147‑1 (MQ=255)
                                                                                                                                aaGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGc            >  4:5360/1‑149 (MQ=255)
                                                                                                                                 aGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCt           <  3:5360/149‑1 (MQ=255)
                                                                                                                                     cTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTggg       <  2:126126/149‑1 (MQ=255)
                                                                                                                                        gACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTgggg      <  3:121552/147‑1 (MQ=255)
                                                                                                                                        atcGGTGACGAAGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTGGGGGt    >  6:120484/3‑149 (MQ=255)
                                                                                                                                           ggTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTGGGGGTgg  >  2:106715/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGGATCAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTGGGGGTGG  >  NC_002947/2444899‑2445185

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: ‑