Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,822,188 G→T 6.5% intergenic (+30/+28) PP_mr22 → / ← fdxA PrrB_RsmZ ncRNA/ferredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,822,1880GT 6.5% 242.6 / 5.4 77intergenic (+30/+28)PP_mr22/fdxAPrrB_RsmZ ncRNA/ferredoxin 1
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base T (2/3);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGAT  >  NC_002947/1822042‑1822335
                                                                                                                                                  |                                                                                                                                                   
aTGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGaa                                                                                                                                                   <  6:314716/149‑1 (MQ=255)
        gCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCttt                                                                                                                                           >  2:98761/1‑149 (MQ=255)
          tttCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCtttt                                                                                                                                          >  2:349268/1‑148 (MQ=255)
           ttCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                        >  4:319659/1‑149 (MQ=255)
            tCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTcg                                                                                                                                       >  6:170817/1‑149 (MQ=255)
            tCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTc                                                                                                                                        <  3:349148/148‑1 (MQ=255)
              aGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTcgc                                                                                                                                      <  3:164669/148‑1 (MQ=255)
              aGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTcgc                                                                                                                                      <  7:34179/148‑1 (MQ=255)
              aGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCACGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttttgt                                                                                                                                      <  4:71998/148‑4 (MQ=255)
                gATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGAGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCTTCGGGCCttttttgtttt                                                                                                                                   <  5:137724/149‑7 (MQ=255)
                     aGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAACCAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgc                                                                                                                              >  8:291028/1‑149 (MQ=255)
                     aGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgc                                                                                                                              >  6:19217/1‑149 (MQ=255)
                        ccAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGgcgc                                                                                                                            >  6:139237/1‑148 (MQ=255)
                         cAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGc                                                                                                                                               >  4:38926/1‑128 (MQ=255)
                         cAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGc                                                                                                                                               <  3:38926/128‑1 (MQ=255)
                          aGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTAt                                                                                                                         <  7:310528/149‑1 (MQ=255)
                            ggACATCGCAAGAAGCGATTTCATCAGTATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCa                                                                                                                       <  8:294989/149‑1 (MQ=255)
                                 tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTAt                                                                                                                         >  5:233278/1‑142 (MQ=255)
                                 tCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTAt                                                                                                                         <  6:233278/142‑1 (MQ=255)
                                     aagaagCGATTTCATCAGGATGATGTTTGGGGCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAg                                                                                                              >  5:271360/1‑149 (MQ=255)
                                     aagaagCGACTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAg                                                                                                              >  7:129595/1‑149 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCg                                                                                                          <  3:243510/146‑1 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCg                                                                                                          >  4:243510/1‑146 (MQ=255)
                                            gATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAAACGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCa                                                                                                                       >  7:281496/1‑133 (MQ=255)
                                            gATTTCATAAGGATGATGTTTGGTACAAGCAGGGACTACTGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCa                                                                                                                       <  8:281496/133‑1 (MQ=255)
                                             aTTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCaa                                                                                                      <  5:19217/149‑1 (MQ=255)
                                              tttCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAAt                                                                                                     <  1:98761/149‑1 (MQ=255)
                                                  atcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCtt                                                                                                  <  1:280477/148‑1 (MQ=255)
                                                  atcaGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGAGCTCCAGGTCCGCAATCTTg                                                                                                 >  8:371839/1‑149 (MQ=255)
                                                    caGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGc                                                                                                >  4:219145/1‑148 (MQ=255)
                                                       gatgatGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAgg                                                                                            >  8:158045/1‑149 (MQ=255)
                                                           atgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                         <  8:129595/148‑1 (MQ=255)
                                                            tgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGc                                                                                       <  1:193013/149‑1 (MQ=255)
                                                            tgTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGc                                                                                       >  8:222688/1‑149 (MQ=255)
                                                              tttGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTCCCGGAAGGAAAAAAAAAAAGGCCCGAAGGGCATTTTACCCTTATGGGCGGTATCAGCCTTCAAGGACGGAAATATTGCAAGGTTTGca                                                                                      >  3:132978/1‑147 (MQ=255)
                                                                 gggTCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATc                                                                                   <  6:271360/148‑1 (MQ=255)
                                                                 gggACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATcc                                                                                  >  4:277128/1‑149 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCa                                                                                              >  7:326184/1‑135 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCa                                                                                              <  8:326184/135‑1 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                                >  4:361294/1‑149 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                                <  3:319659/149‑1 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                                <  7:228946/149‑1 (MQ=255)
                                                                   gACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCa                                                                                <  7:158045/149‑1 (MQ=255)
                                                                    aCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAACCAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAt                                                                               <  7:291028/149‑1 (MQ=255)
                                                                    aCAAGCAGGGACTACGGAAAAAAATGTGGGCGGGGCAAACCGCCCCTTTTTTTGCCATAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAtt                                                                              <  6:305450/149‑1 (MQ=255)
                                                                      aaGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttc                                                                             <  7:44415/149‑1 (MQ=255)
                                                                         cAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCAttctt                                                                           >  5:281029/1‑148 (MQ=255)
                                                                           gggACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAg                                                                        <  3:361294/149‑1 (MQ=255)
                                                                            ggACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGc                                                                       >  8:149865/1‑149 (MQ=255)
                                                                              ccTACGGGCAAAAAAAGGGGCGGTTTGACCCGCCCACATTTTTTTCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                         >  2:254934/2‑129 (MQ=18)
                                                                              ccTACGGGCAAAAAAAGGGGCGGTTTGACCCGCCCACATTTTTTTCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCtt                                                                                         <  1:254934/128‑1 (MQ=18)
                                                                                  cGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCgg                                                                 <  6:281029/149‑1 (MQ=255)
                                                                                   ggAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCgg                                                                 >  6:11093/1‑148 (MQ=255)
                                                                                      aaaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAgg                                                             >  1:136954/1‑149 (MQ=255)
                                                                                       aaaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTAGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAggg                                                            >  2:106928/1‑149 (MQ=255)
                                                                                        aaaaTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGc                                                           <  5:139237/149‑1 (MQ=255)
                                                                                           aTGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTc                                                        >  2:376717/1‑149 (MQ=255)
                                                                                              tgGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCt                                                      <  6:375525/148‑1 (MQ=255)
                                                                                              tgGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCt                                                      <  5:11093/148‑1 (MQ=255)
                                                                                                ggCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTa                                                    <  7:96871/148‑1 (MQ=255)
                                                                                                ggCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTAc                                                   <  3:219145/149‑1 (MQ=255)
                                                                                                        aaaCCGCCCCTTTTTTTGCCAGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGt                                            >  7:58724/1‑148 (MQ=255)
                                                                                                         aaCCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTTGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGa                                          <  1:106928/149‑1 (MQ=255)
                                                                                                            cGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGatat                                       >  4:212948/1‑149 (MQ=255)
                                                                                                              ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTa                                     <  2:392751/149‑1 (MQ=255)
                                                                                                              ccccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTa                                     <  7:222688/149‑1 (MQ=255)
                                                                                                               cccTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAg                                    >  3:3110/1‑149 (MQ=255)
                                                                                                                          cccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTcc                         >  8:183301/1‑149 (MQ=255)
                                                                                                                          cccGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTcc                         >  4:350082/1‑149 (MQ=255)
                                                                                                                            cGTAGGAAAACAAAAAGGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGc                       <  8:259065/149‑1 (MQ=255)
                                                                                                                                     gcgaaaaaGGCCCTTCGGGCCTTTTTTGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCg                                              >  2:211397/4‑117 (MQ=255)
                                                                                                                                     gcgaaaaaGGCCCTTCGGGCCTTTTTTGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCg                                              <  1:211397/114‑1 (MQ=255)
                                                                                                                                     aCAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGc              >  4:136341/1‑149 (MQ=255)
                                                                                                                                       gaaaaaGGCCCTTCGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGt            >  5:265432/2‑149 (MQ=255)
                                                                                                                                        aaaaaGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGtt           <  4:132978/149‑1 (MQ=255)
                                                                                                                                               cccGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCg    <  2:136954/149‑1 (MQ=255)
                                                                                                                                                 cGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGAt  <  3:212948/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATGCACAGGCTTTCAGGATGAAGGCCAGGGACATCGCAAGAAGCGATTTCATCAGGATGATGTTTGGGACAAGCAGGGACTACGGAAAAAAATGTGGGCGGGTCAAACCGCCCCTTTTTTTGCCCGTAGGAAAACAAAAAAGGCCCGAAGGGCCTTTTTCGCGTCTGGGCGCTATCAGCGCTCCAGGTCGGCAATCTTGCCAGGCTTGCCATCCCATTCTTCAGCGTCCGGCAGGGCGTCCTTACGCTCGGTGATATTAGGCCAGATTTCCGCCAGCTCGGCGTTCAGCTCGAT  >  NC_002947/1822042‑1822335

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: