Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,973,872	A→T	9.4%	intergenic (+29/‑119)	gyrA → / → serC	DNA gyrase subunit A/3‑phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,973,872	0	A	T	9.4%	170.3 / 8.3	64	intergenic (+29/‑119)	gyrA/serC	DNA gyrase subunit A/3‑phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base T (3/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

GCCGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTGCG  >  NC_002947/1973726‑1974018                                                                                                                                                    |                                                                                                                                                   gCCGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAgg                                                                                                                                                  >  4:19961/1‑149 (MQ=255)       cGAGGACTAACTCGATGAAGTGAGGGAGACGTACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGcc                                                                                                                                             <  6:203946/148‑1 (MQ=255)       cGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGcc                                                                                                                                             <  4:156847/148‑1 (MQ=255)        gAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCTTGGTCGcccc                                                                                                                                           >  3:309056/1‑149 (MQ=255)        gAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGcccc                                                                                                                                           >  2:297837/1‑149 (MQ=255)              gAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCTTGGTCGCCCCgctcgg                                                                                                                                     <  4:309056/149‑6 (MQ=255)                     aTGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGAGTAAACCCGAGCGGGGC‑GACCTAGGTCGCCCCGTTTGACTGTGTa                                                                                                                              <  6:278006/149‑1 (MQ=255)                      tGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAg                                                                                                                             <  1:182454/149‑1 (MQ=255)                         aGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCTCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGa                                                                                                                           <  4:350794/148‑1 (MQ=255)                         aGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGa                                                                                                                           <  4:345280/148‑1 (MQ=255)                             aTCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCgg                                                                                                                       <  3:19961/148‑1 (MQ=255)                               cGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCt                                                                                                                    >  5:329295/1‑149 (MQ=255)                                  gacggacgAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGc                                                                                                                 >  3:36647/1‑149 (MQ=255)                                            gaggGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGc                                                                                                        <  4:210529/148‑1 (MQ=255)                                                 tgtCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGgcgc                                                                                                  >  5:101768/1‑149 (MQ=255)                                                    cGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCAc                                                                                               <  7:133819/149‑1 (MQ=255)                                                        gCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGAt                                                                                           >  3:289727/1‑149 (MQ=255)                                                         cTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATc                                                                                          >  8:144576/1‑149 (MQ=255)                                                              gcgcCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCa                                                                                     <  4:57824/149‑1 (MQ=255)                                                                 ccAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCTTGGTCGCCCCgctcg                                                                                                                                      <  5:151354/97‑5 (MQ=39)                                                                 ccAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCTTGGTCGCCCCgctcg                                                                                                                                      >  6:151354/1‑93 (MQ=39)                                                                 ccAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATcc                                                                                         <  5:168567/142‑1 (MQ=255)                                                                 ccAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATcc                                                                                         >  6:168567/1‑142 (MQ=255)                                                                    gACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAAt                                                                               <  6:329295/149‑1 (MQ=255)                                                                       aatgaaaatgaaGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAAtgtg                                                                            <  4:289727/149‑1 (MQ=255)                                                                         tgaaaatgaaGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGc                                                                          >  3:171578/1‑149 (MQ=255)                                                                                   gaagCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTa                                                                >  5:66336/1‑149 (MQ=255)                                                                                       cTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGAt                                                            >  6:103888/1‑149 (MQ=255)                                                                                        tGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGAtt                                                           >  3:93137/1‑149 (MQ=255)                                                                                           aGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGTCTTAACTAGGTCCCCCCCTATGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGAtttt                                                         <  1:297837/149‑1 (MQ=255)                                                                                             gTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGtt                                                      <  4:140005/149‑1 (MQ=255)                                                                                              tGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGtt                                                      <  6:101768/148‑1 (MQ=255)                                                                                                   acgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAAtgtg                                                                            <  6:264402/121‑1 (MQ=255)                                                                                                   acgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAAtgtg                                                                            >  5:264402/1‑121 (MQ=255)                                                                                                   acgaGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTAttt                                                <  4:36647/149‑1 (MQ=255)                                                                                                       ggCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCAACGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTgg                                              >  2:45934/1‑147 (MQ=255)                                                                                                        gCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGAAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTgg                                              >  1:289206/1‑146 (MQ=255)                                                                                                          cccGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCagag                                         >  4:245711/1‑149 (MQ=255)                                                                                                           ccGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGc                                        >  3:222498/1‑149 (MQ=255)                                                                                                               aaGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCgaga                                    <  2:51977/149‑1 (MQ=255)                                                                                                                agagTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGagag                                   <  5:103888/149‑1 (MQ=255)                                                                                                                 gagTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGt                                  <  4:93137/149‑1 (MQ=255)                                                                                                                    taagtaagCGTAAACCCGAGCGGGGC‑GCCCAAGGTCGCCCCCTTTGACTGTGTAGGAGCGGGCTTTCCCGCGCTGCAGGCGCCACCGAACCTGCATGGAATGTGGGCCCCGAAATAAGATTTTGTTCTATTTGGGAGAGCGAGAGTggc                               >  2:341569/1‑148 (MQ=255)                                                                                                                          agCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGa                          <  3:245711/148‑1 (MQ=255)                                                                                                                               aaaCCCGAGCGGGGC‑GACCAAGGTCGCCCCCTTTGACTGTGTAGGAGCGGGCTTTCCCGCGATGCAGGCGCCACCGATACTGCATGGAATGTGGCCACCGAATTAAGATTTTGTTAAATTTGGCAGAGCGAGAGTGGATGTGAGCAAAc                    >  6:14209/1‑149 (MQ=255)                                                                                                                                 aCCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGa                  <  4:171578/149‑1 (MQ=255)                                                                                                                                   ccgcccGGGGC‑GACCAAGCCCGCCCCGGTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGCATGTGAGCAAACGAGc                <  2:289206/144‑1 (MQ=255)                                                                                                                                   ccGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTTGGAGCGGGCTTGCCCGCGATGCAGGCGCCCCCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGc                >  2:183977/1‑149 (MQ=255)                                                                                                                                   ccGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGc                >  6:257563/1‑149 (MQ=255)                                                                                                                                       acGGGGC‑GACCTTGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCttt            <  6:15975/148‑1 (MQ=255)                                                                                                                                       acGGGGC‑GACCTTGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCtt             >  5:15975/2‑148 (MQ=255)                                                                                                                                        cGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCttt            <  5:257563/148‑1 (MQ=255)                                                                                                                                         ggggC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTa           >  7:82743/1‑148 (MQ=255)                                                                                                                                          gggC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATTCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCga                                      <  1:272443/120‑1 (MQ=255)                                                                                                                                          gggC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCga                                      >  2:277170/1‑120 (MQ=255)                                                                                                                                          gggC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCga                                      >  2:272443/1‑120 (MQ=255)                                                                                                                                          gggC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCga                                      <  1:277170/120‑1 (MQ=255)                                                                                                                                          gggC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAAc         >  8:83753/1‑149 (MQ=255)                                                                                                                                            gC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACtt       <  7:83753/149‑1 (MQ=255)                                                                                                                                               gaacAAGGTCGGCCCGTTTGACTGTGTAGGAGCGGGCTTTCCCGCGATGCAGGCTCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTa           <  8:171831/140‑1 (MQ=255)                                                                                                                                               gACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTa           >  7:171831/1‑143 (MQ=255)                                                                                                                                                aCCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTg    >  6:86696/1‑149 (MQ=255)                                                                                                                                                  cAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTg                                 <  7:307892/118‑1 (MQ=255)                                                                                                                                                  cAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTg                                 >  8:307892/1‑118 (MQ=255)                                                                                                                                                   aaGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCATCGATCCAGCAAGGAATGTTGCCACCGAAATAAGAATTTGTTCGATTTTGCAGAGCGAGAGTGGAAGTGAGCAAACGAGCCTTTAAcatattcg  >  1:93942/1‑141 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCGTCCGAGGACGAACTCGATGAAGTGATCGAGACGGACGAGGAGGGTGTCGAGGCTGTAGCGCCAGACAATGAAAATGAAGAAGCTGCAGGTGCCGACGAGGCCCCGCAAGAGTAAGTAAGCGTAAACCCGAGCGGGGC‑GACCAAGGTCGCCCCGTTTGACTGTGTAGGAGCGGGCTTGCCCGCGATGCAGGCGCCACCGATCCTGCATGGAATGTGGCCACCGAATTACGATTTTGTTCTATTTGGCAGAGCGAGAGTGGATGTGAGCAAACGAGCCTTTAACTTCTGCG  >  NC_002947/1973726‑1974018

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