Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	903,624	A→G	6.8%	intergenic (‑126/‑118)	nadC ← / → PP_0788	quinolinate phosphoribosyltransferase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	903,624	0	A	G	6.8%	172.2 / 2.1	60	intergenic (‑126/‑118)	nadC/PP_0788	quinolinate phosphoribosyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (2/2);  total (30/30)
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.40e-01

GTAGGTTCGGCATGGTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAAAGAAGGCAAGGTGGTGCCCTTGG  >  NC_002947/903486‑903772                                                                                                                                            |                                                                                                                                                     gTAGGTTCGGCATGGTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCTGATCAGGGTCATTAGC‑CAGCAAGCGAAACGgc                                                                                                                                            >  8:255869/1‑149 (MQ=255) gTAGGTTCGGCATGGTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGgc                                                                                                                                            >  5:300583/1‑149 (MQ=255)             tGGTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTAt                                                                                                                                >  2:73620/1‑149 (MQ=255)             tGGTCGGCGCCACAAGCTTGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTAt                                                                                                                                <  4:280727/149‑1 (MQ=255)              ggTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATg                                                                                                                               >  6:173401/1‑149 (MQ=255)              ggTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATg                                                                                                                               >  6:80729/1‑149 (MQ=255)                  ggCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGAt                                                                                                                           >  4:9730/1‑149 (MQ=255)                  ggCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGAt                                                                                                                           >  8:151108/1‑149 (MQ=255)                      ccACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGAttttt                                                                                                                       >  8:307986/1‑149 (MQ=255)                       cacaAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGTGAAACGGCGCTGGTAGCTATGCGAttttt                                                                                                                       <  3:53919/148‑1 (MQ=255)                       aacaAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGAttttt                                                                                                                       >  4:323580/2‑148 (MQ=255)                        acaAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGAtttttt                                                                                                                      >  1:150298/1‑148 (MQ=255)                          aaGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTAcc                                                                                                                   >  6:72680/1‑149 (MQ=255)                              tAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCCGCGAAACGGCGCTGGAAGCTATGCGATTTTTTACCGAt                                                                                                                >  2:158170/1‑148 (MQ=255)                                ggTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCTGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATaa                                                                                                              <  7:255869/148‑1 (MQ=255)                                    cccGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGc                                                                                                          >  3:185450/1‑148 (MQ=255)                                    cccGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGc                                                                                                          <  5:80729/148‑1 (MQ=255)                                    cccGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGc                                                                                                          <  7:79140/148‑1 (MQ=255)                                     ccGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGa                                                                                                        >  8:16037/1‑149 (MQ=255)                                      cgcgATTATCGCGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGTCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGAc                                                                                                       <  2:62221/149‑1 (MQ=255)                                      cgcgATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGAc                                                                                                       <  7:348491/149‑1 (MQ=255)                                      cgcgATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGAc                                                                                                       <  7:16037/149‑1 (MQ=255)                                      cgcgATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTATGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCCAc                                                                                                       >  5:64450/1‑149 (MQ=255)                                          aTTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGt                                                                                                   <  7:307986/149‑1 (MQ=255)                                              taGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGATTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGAGAAACGGCGCTGGGAGCTATGCGATTTTTTACCGATAAGGGCGACTAGTAttt                                                                                               >  4:126512/1‑149 (MQ=255)                                              taGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTAttt                                                                                               >  5:245236/1‑149 (MQ=255)                                                gggCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTgg                                                                                             <  1:73620/149‑1 (MQ=255)                                                       gggCCGGTGTGTACAGCCGCCGTGGCGCTGACCTTTGGTCGGGCGGCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCACCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATa                                                                                      <  2:324730/149‑1 (MQ=255)                                                                 gtACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGc                                                                                                          >  6:294661/1‑119 (MQ=255)                                                                 gtACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGc                                                                                                          <  5:294661/119‑1 (MQ=255)                                                                     aGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGAtttt                                                                        >  8:10972/1‑149 (MQ=255)                                                                      gccgccTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGAtttt                                                                        >  4:331274/1‑148 (MQ=255)                                                                       ccgccTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGAttttt                                                                       >  8:6138/1‑148 (MQ=255)                                                                        cgccTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCa                                                                     <  3:9730/149‑1 (MQ=255)                                                                        cgccGTGGCGGTGACCGTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATGAGGGGCATTAGC‑CAGCCAGCGAAACGGCGCTGGTATCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCa                                                                     <  4:49218/149‑1 (MQ=255)                                                                        agccTTGGTGGTGACCTTTGGTCGGGCGTCGTTTGTATTCGGATCAGATCAGGGTCATTAGG‑CAGCCAGCGACACGGCGCTGGTAGCTATGCGATTTTTTACCGTTAATGGCGACTAGTATTTGGGGTCATAGCTTTGACGATTTTTCa                                                                     <  4:94926/148‑1 (MQ=255)                                                                         gccTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAg                                                                    <  4:185450/149‑1 (MQ=255)                                                                                cgcTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGGCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTc                                                             <  2:200508/149‑1 (MQ=255)                                                                                  cTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGCGCAGCCAGCGAAACGGCGCTGGTAGATATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCg                                                            >  1:98932/1‑149 (MQ=255)                                                                                  cTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCgg                                                           <  5:72680/149‑1 (MQ=255)                                                                                     aCCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCggg                                                          <  3:323580/147‑1 (MQ=255)                                                                                     aCCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGc                                                         >  7:233631/1‑148 (MQ=255)                                                                                     aCCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGc                                                         >  7:152583/1‑148 (MQ=255)                                                                                       cTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAg                                                      <  1:175626/149‑1 (MQ=255)                                                                                        tttGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACAGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGa                                                     >  3:257204/1‑149 (MQ=255)                                                                                          tGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGac                                                    <  3:331274/148‑1 (MQ=255)                                                                                                         aGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTcc                                    <  7:10972/149‑1 (MQ=255)                                                                                                           tAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCGGCGAACCGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAg                                  >  7:250787/1‑149 (MQ=255)                                                                                                           tAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCGGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAg                                  >  7:255670/1‑149 (MQ=255)                                                                                                           tAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGTCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAg                                  >  8:248967/1‑149 (MQ=255)                                                                                                            aaTCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAg                                  <  2:98932/148‑1 (MQ=255)                                                                                                             aTCGGCGCAGATCAGGGTCAGTAGC‑CAGCCGGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGa                                <  8:250787/149‑1 (MQ=255)                                                                                                              tCGGAGCAGATCAGGGTCATTAGC‑CAGCCGGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGa                                <  8:255670/148‑1 (MQ=255)                                                                                                                 gAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGc                             >  1:55111/1‑148 (MQ=255)                                                                                                                 gAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGc                             >  4:119409/1‑148 (MQ=255)                                                                                                                     aGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCaaaaa                        <  7:151108/149‑1 (MQ=255)                                                                                                                        tCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAaagaa                     <  6:245236/149‑1 (MQ=255)                                                                                                                        tCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAaagaa                     <  2:161593/149‑1 (MQ=255)                                                                                                                        tCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAaagaa                     <  6:300583/149‑1 (MQ=255)                                                                                                                          aGGGTCATTAGC‑CAGCCATCGAAACGGCGCTTGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAATGAGTCCAGGATGCAAAAAGAAgg                   <  3:252301/149‑1 (MQ=255)                                                                                                                                    gc‑cagccaGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAAAGAAGGCAAGgtggtg         >  7:335013/1‑149 (MQ=255)                                                                                                                                            agcGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAAAGAAGGCAAGGTGGTGCCCTTgg  >  3:13789/1‑149 (MQ=255)                                                                                                                                            agcGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAAAGAAGGCAAGGTGGTGCCCTTgg  >  3:93467/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                     GTAGGTTCGGCATGGTCGGCTCCACAAGCTAGGTGCCCGCGATTATAGGGCAGGGGGCCGGTGTGTACAGCCGCCTTGGCGCTGACCTTTGGTCGGGCGTCGTGAGTAATCGGAGCAGATCAGGGTCATTAGC‑CAGCCAGCGAAACGGCGCTGGTAGCTATGCGATTTTTTACCGATAATGGCGACTAGTATTTGGCGTCATAGCTTTGACGATTTTTCAGCCCTTTCGGGCAAGACACCCTGCAAGGAGTCCAGGATGCAAAAAGAAGGCAAGGTGGTGCCCTTGG  >  NC_002947/903486‑903772

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