Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,445,044 T→G 6.0% E156D (GAA→GAC) ‡ lexA‑I ← transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,445,0440TG 6.0% 208.8 / 4.4 67E156D (GAA→GAC) ‡lexA‑Itranscriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base G (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 = 9.11e-01

TCCTGGATCAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTGGGGGTGG  >  NC_002947/2444896‑2445185
                                                                                                                                                    |                                                                                                                                             
tcctGGATCAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACt                                                                                                                                               >  4:189555/1‑149 (MQ=255)
   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)
                                                                                                                                                    |                                                                                                                                             
TCCTGGATCAGCGGCGAATGACGCCGACGCTCAAGCCCTCGATCACCAGCTCCTGTTCTTTCAGGTCGACTTCGATGGGGGCGAATTCGGGGTTTTCGGCAAGCAGCCAGACTTTGCTGCCTTCGCGCTTGAAGCGCTTGACGGTGACTTCGTCACCGATGCGCGCCACAACGATCTGGCCGTTGCGGGCTTCACGGCAGGTGTGCACCGCCAGCAGGTCGCCGTCGAAGATGCCGACGTCCTTCATGCTCATGCCGTGTACGCGCAGCAAATAGTCGGCTTGGGGGTGG  >  NC_002947/2444896‑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: ‑