Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,816 T→G intergenic (+53/+700) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8160TG96.4% 172.2 / ‑3.0 57intergenic (+53/+700)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (24/30);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGC  >  NC_002947/4740679‑4740958
                                                                                                                                         |                                                                                                                                                
cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                                      >  2:221320/1‑149 (MQ=255)
cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                                      <  5:227523/149‑1 (MQ=255)
         cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                              <  4:429206/148‑1 (MQ=255)
            gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                                          >  5:146261/1‑149 (MQ=255)
                ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                                       <  1:253223/148‑1 (MQ=255)
                  ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                                    <  6:272306/149‑1 (MQ=255)
                       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGTCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                               <  8:16202/149‑1 (MQ=255)
                       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                               <  6:185386/149‑1 (MQ=255)
                       cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                               <  3:215429/149‑1 (MQ=255)
                             aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                         <  3:216245/149‑1 (MQ=255)
                              gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCg                                                                                                        <  5:425283/149‑1 (MQ=38)
                              gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCg                                                                                                        <  4:106205/149‑1 (MQ=38)
                                      ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTGGCGGGCGTGCCCGCGAAGAGGGCGGGacaca‑‑c                                                                                              >  4:255028/1‑149 (MQ=255)
                                       cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGCCCGCGAAGAGACCGGCACAGa‑‑ca                                                                                             >  8:202682/1‑149 (MQ=255)
                                                 aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAAcc                                                                                   >  4:309390/1‑149 (MQ=255)
                                                        aCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGttt                                                                            >  6:185909/1‑149 (MQ=255)
                                                              ttGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTTCT‑GTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGAGGCCGGCGAAGA‑‑CAACACAAAACCTCAGgttttatac                                                                      >  5:270796/1‑145 (MQ=255)
                                                               tGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTcc                                                                      <  5:185909/148‑1 (MQ=255)
                                                                 cccTCAAAGACCCCAAGTACGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTccgcc                                                                   <  6:146261/149‑1 (MQ=255)
                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTATccgccc                                                                  >  5:257186/1‑149 (MQ=255)
                                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTATccgccc                                                                  >  7:249969/1‑149 (MQ=255)
                                                                   cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCGTGGCCGCGCAGAGGCCGGCCCAGC‑‑CACCACCAAACCGAAGTTTTTGTacgcccc                                                                 >  8:223419/1‑148 (MQ=255)
                                                                    tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTccgcccgc                                                                >  6:345434/1‑149 (MQ=255)
                                                                     gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTGTTTCTCCGCCCGCt                                                               >  5:218562/1‑149 (MQ=255)
                                                                      aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGACGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGA‑‑CACCACAAAACCTCAGTTTTTCTccgcccgca                                                               >  5:163103/1‑147 (MQ=255)
                                                                       aaGACCGCAAGTAAGCTGCACTAGCCGAAAGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTAGTCCTGTGCCGCCCTCATCGCGGACTTCCCCCCGAAGAGGCCGGCACACA‑‑AACCACAAAACCTAACTTATAATCCGCCCGCTCa                                                             >  4:276714/1‑149 (MQ=255)
                                                                            cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                         <  3:130411/102‑1 (MQ=39)
                                                                            cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                         >  4:130411/1‑102 (MQ=39)
                                                                            cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCGCAGGTTTTGTACGCCCGCGCACGCAg                                                        >  7:236334/1‑149 (MQ=255)
                                                                                aGTATGCTGCTGTTGCCTAACGCAAGAAGGCTGCCTTGGGTCAG‑CTTTTTTTATTTGTCTTGTGCCGGGCTCTTCGCTGCCGTGCCCGCGATGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGACCGCTCACGCAGCCCt                                                    <  7:202682/149‑1 (MQ=255)
                                                                                aGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCt                                                    <  3:309390/149‑1 (MQ=255)
                                                                                  tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGC‑‑CACCACAAAACCTCAGGTTTTCTCCGCCCGCTCACCCAGcccat                                                   >  5:461182/1‑146 (MQ=255)
                                                                                    cgCTGCAGTAGCCCAACTCAAGACGGGGGCCGGCGGGC‑TTCGTTTTTTATTTGTCGTGGGCCGGCCTCTTCGCGGCCTGGCCGGGGATGTGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                <  6:163103/148‑1 (MQ=255)
                                                                                     gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGAGTTGGTGCAAAACAACAGATTTTCTCAGCCCGCTCACTCACCCCtat                                                  >  1:387741/1‑146 (MQ=255)
                                                                                      cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                              <  5:345434/149‑1 (MQ=255)
                                                                                      cTGCAGTACCCTAACGCAAGATGGCTTCCCTCGGGCAGCGTTTTTTT‑TTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                              <  2:387741/149‑1 (MQ=255)
                                                                                            tAGCCTACCGCAATAAGGGTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                         <  8:249969/148‑1 (MQ=255)
                                                                                                  accGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                  <  7:223419/147‑1 (MQ=255)
                                                                                                  aaCGCAAGACGGCGGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTCCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                  <  1:276732/149‑1 (MQ=255)
                                                                                                  aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                  <  1:243057/149‑1 (MQ=255)
                                                                                                       aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAACAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                             >  7:31201/1‑149 (MQ=255)
                                                                                                            ggCTGCCCTCGTGCAGCC‑TTTTTTATTTGTCCGTTCCCGGCCCCTGCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                        <  2:158104/149‑1 (MQ=255)
                                                                                                            ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                              <  1:56316/127‑1 (MQ=255)
                                                                                                            ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                              >  2:56316/1‑127 (MQ=255)
                                                                                                             gCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCACAGTTTTTATACCCCCCCTAACCCAGCCCGTTAACGGTATTGAACCGCGCAaccacccc                       >  8:114869/1‑147 (MQ=255)
                                                                                                                 cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                   <  7:114869/149‑1 (MQ=255)
                                                                                                                 cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                   >  6:130747/1‑149 (MQ=255)
                                                                                                                      gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg                <  5:130747/147‑1 (MQ=255)
                                                                                                                      ggccAGCC‑TTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt              <  2:328566/146‑1 (MQ=255)
                                                                                                                       ggccagcTTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCTGCCTTGCCCGGGAAGTGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAAGCAGCCCTTTCAGGGTATTGAACGGCTCATCCACCACGAACttg                <  5:234057/141‑1 (MQ=255)
                                                                                                                       ggCAGCTTTTTTTTTTTGCCCTGG‑GCCGCCCTCTTCGCGCCGTGGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg             <  3:255028/149‑1 (MQ=255)
                                                                                                                        gCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCAAAGTTTTTCTCCGCCCGCTCCCGCAGCCCATTCAGGGTATTGAACGGCGCATCCACCACCAACTTGTTgg            >  8:10567/1‑149 (MQ=255)
                                                                                                                          ttggttttttttttttTCCTGTGCGGCCCTTTTCGGGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc           <  7:173984/145‑1 (MQ=255)
                                                                                                                          aagcctttttaATTTTTCCTGTGCCGGCCTCTTCCCCCCCTATCCCCCGAAGCGGCCGGCACAGA‑‑CACCACAAAACCTCAGAATATCACCGCCCCCGCACGCAGCCCCTTGAAGCTATTGAAGCCCGAAACCACCACGAACATGTAGGc           >  8:196394/6‑149 (MQ=255)
                                                                                                                           agcgttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcc          <  7:10567/145‑1 (MQ=255)
                                                                                                                              cttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGA‑‑CACCCCAAAACCGCCGTTTTTCTCCGCCCGCTCCCGCAGCCCGTTCAGGGGATTGAACGGCGCCTACACCACGAACTTGTTGGccac        >  8:130629/2‑148 (MQ=255)
                                                                                                                              cttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc       <  6:218562/148‑1 (MQ=255)
                                                                                                                                   ttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGc  >  5:326188/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                
CAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGA‑‑CACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGC  >  NC_002947/4740679‑4740958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑