Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,332,245	C→G	6.8%	intergenic (+28/‑21)	dksA → / → gluQ	RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,332,245	0	C	G	6.8%	183.3 / 3.6	59	intergenic (+28/‑21)	dksA/gluQ	RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base G (2/2);  total (32/27)
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.42e-01

GTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAAC  >  NC_002947/5332097‑5332391                                                                                                                                                      |                                                                                                                                                   gTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGcatc                                                                                                                                                    <  3:20442/149‑1 (MQ=255) gTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTCAACGGGGCGcatc                                                                                                                                                    <  6:299597/149‑1 (MQ=255)    cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGcatcat                                                                                                                                                  <  7:310229/148‑1 (MQ=255)        tCGTGCGGCATCGAGATCGGCATGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGcc                                                                                                                                             <  4:192166/149‑1 (MQ=255)          gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                                           <  8:37052/149‑1 (MQ=255)          gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGccc                                                                                                                                            <  1:13855/148‑1 (MQ=255)           caaggCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                                          <  5:215292/146‑1 (MQ=255)                cATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCa                                                                                                                                     >  6:48582/1‑149 (MQ=255)                 aTCGCGATCGGCCCGCGCCGCCTGGCACCCCGCCCTACCCCCGACCTGTGCATCCCCTGCAAGGAAATTGCCCAGAAAAACGCTAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATGCGCCCCGttcag                                                                                                                                     <  1:297869/149‑4 (MQ=255)                  tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                                                   <  8:162304/149‑1 (MQ=255)                    gagaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                                                  <  1:138707/148‑1 (MQ=255)                     agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTta                                                                                                                                 >  8:300434/1‑148 (MQ=255)                       aTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTatat                                                                                                                              <  7:167008/149‑1 (MQ=255)                       aTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTatat                                                                                                                              <  2:57208/149‑1 (MQ=255)                                        ggAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGa                                                                                                             >  5:174442/1‑149 (MQ=255)                                         gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAATAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAc                                                                                                            <  7:300434/149‑1 (MQ=255)                                          aaGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                                           <  6:335077/149‑1 (MQ=255)                                           aGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATGTTTCTTCCACCTGAACGGGGCGCATCATGGGCCCCGTTTCATTtatc                                                                                                                               >  4:227460/1‑127 (MQ=255)                                             cccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGa                                                                                                         >  1:134827/1‑148 (MQ=255)                                              ccgcccTACCGCCGAGCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACt                                                                                                       >  6:40757/1‑149 (MQ=255)                                                gcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTcc                                                                                                     <  2:353404/149‑1 (MQ=255)                                                                 tgCTTCGACTGCAAGGAAATTGCCGAGAATAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCt                                                                                     <  2:151398/148‑1 (MQ=255)                                                                   cTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTc                                                                                   <  5:48582/148‑1 (MQ=255)                                                                      cGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGccc                                                                               >  2:265105/1‑149 (MQ=255)                                                                          tGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccac                                                                           >  8:182283/1‑149 (MQ=255)                                                                            cAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccc                                                                         >  2:306863/1‑149 (MQ=255)                                                                             aaGGAAATTGTCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGccccccaccc                                                                        >  3:312129/1‑144 (MQ=255)                                                                             aaGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccacccc                                                                        >  2:252506/1‑149 (MQ=255)                                                                             aaGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccacccc                                                                        >  2:270975/1‑149 (MQ=255)                                                                                       ccGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTc                                                              <  7:311888/149‑1 (MQ=255)                                                                                       ccGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTc                                                              <  3:18470/149‑1 (MQ=255)                                                                                       ccGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTc                                                              >  1:360439/1‑149 (MQ=255)                                                                                        cGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTcc                                                             >  7:180598/1‑149 (MQ=255)                                                                                         gagaAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACAACGGGCGCTTAGACCCCACCCCCAGCGGCTTCCt                                                            >  7:185262/1‑149 (MQ=255)                                                                                          agaAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGTCTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                           <  1:270975/149‑1 (MQ=255)                                                                                          agaAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                           <  5:40757/149‑1 (MQ=255)                                                                                            aaaaa‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCa                                                         >  3:186294/1‑149 (MQ=255)                                                                                              aaa‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACt                                                       >  7:305533/1‑149 (MQ=255)                                                                                                  ggAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCg                                                    <  7:302631/149‑1 (MQ=255)                                                                                                   gAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCg                                                                   >  7:201118/1‑133 (MQ=255)                                                                                                   gAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCg                                                                   <  8:201118/133‑1 (MQ=255)                                                                                                   gAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCgg                                                   >  5:220617/1‑149 (MQ=255)                                                                                                     aaaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGc                                                  >  5:362856/1‑148 (MQ=255)                                                                                                        aaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCCCCCCCAGCGGCTTCCTGCCCTTCGGCTCGc                                              >  2:138320/1‑149 (MQ=255)                                                                                                                gCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgcc                                      >  6:82151/1‑149 (MQ=255)                                                                                                                  aaaGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATc                                                                                            <  8:184099/93‑1 (MQ=255)                                                                                                                  aaaGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATc                                                                                            >  7:184099/1‑93 (MQ=255)                                                                                                                     ggCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCt                                 >  2:47870/1‑149 (MQ=255)                                                                                                                        tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                               >  5:78801/1‑68 (MQ=255)                                                                                                                        tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                                               <  6:78801/68‑1 (MQ=255)                                                                                                                                   ccACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                     <  6:174442/147‑1 (MQ=255)                                                                                                                                   ccACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCt                     <  1:306863/147‑1 (MQ=255)                                                                                                                                     aCCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                 <  1:252506/149‑1 (MQ=255)                                                                                                                                       tgaaaCGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGc               >  6:274371/4‑149 (MQ=255)                                                                                                                                            cGGGGCGCATGATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGACTCGATGCCcgcg          >  1:341725/1‑149 (MQ=255)                                                                                                                                            cGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGAt                 <  8:248926/142‑1 (MQ=255)                                                                                                                                            cGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCagcg          >  1:370783/1‑149 (MQ=255)                                                                                                                                            cGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGACTCCTGGCTCGAt                 >  7:248926/1‑142 (MQ=255)                                                                                                                                             ggggCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCgcgc         >  6:347558/1‑149 (MQ=255)                                                                                                                                                gcgcATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                           <  7:128494/96‑1 (MQ=255)                                                                                                                                                gcgcATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTg                                                           >  8:128494/1‑96 (MQ=255)                                                                                                                                                  gcATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTc     >  3:361811/1‑148 (MQ=255)                                                                                                                                                    atcatGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAAc  >  5:182810/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAA‑GGAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTCCTGGCTCGATGCCCGCGCCGTCAAC  >  NC_002947/5332097‑5332391

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