Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,657,745	A→T	6.0%	intergenic (‑51/‑201)	tktA ← / → sahR	transketolase A/methionine metabolism transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,657,745	0	A	T	6.0%	204.2 / 2.8	67	intergenic (‑51/‑201)	tktA/sahR	transketolase A/methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base T (2/2);  total (37/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.54e-01

ATGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAA  >  NC_002947/5657603‑5657891                                                                                                                                               |                                                                                                                                                   aTGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCg                                                                                                                                              >  3:185895/1‑149 (MQ=255)         aCCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCCTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTgg                                                                                                                                      <  5:118475/149‑1 (MQ=255)           cTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTgggg                                                                                                                                    >  6:247026/1‑149 (MQ=255)              ggTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGgcg                                                                                                                                 >  8:198462/1‑149 (MQ=255)                  gCCGCTGTTGGCCTGTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGGGGTGTCACGAAAAAGGCGGCCATTTTCGCCCCCAGTGGgcgaaaa                                                                                                                             <  7:275306/149‑5 (MQ=255)                    cGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTggg                                                                                                                           >  6:183643/1‑149 (MQ=255)                    cGCTGTTGGCCTTTTGCACGGCATCAATGCTGGGGGCACGAATGGCGTTGGCCGGTTGACGGCGGCTGGGCATGGCTGAGCTCCTAGGGGGGTTGAATAGGTGGTGTCACGAAAATGGCGGCCATGTTCGCCCACTGGGGGCGCTTggg                                                                                                                           >  6:338867/1‑149 (MQ=255)                     gCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCCCAGTGGgcgaaaatgg                                                                                                                          <  8:186553/149‑8 (MQ=255)                         ttGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCaaa                                                                                                                      <  6:200477/149‑1 (MQ=255)                           ggCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACa                                                                                                                    >  3:274318/1‑149 (MQ=255)                               ttttGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACAttt                                                                                                                 >  6:49242/1‑148 (MQ=255)                                 ttGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaa                                                                                                               <  5:324801/148‑1 (MQ=255)                                    cACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgc                                                                                                           >  3:125924/1‑149 (MQ=255)                                       ggCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgccaa                                                                                                        <  7:198462/149‑1 (MQ=255)                                        gCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgccaag                                                                                                       >  3:178702/1‑149 (MQ=255)                                          aTCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTaagcgccaagcg                                                                                                     >  5:133800/1‑149 (MQ=255)                                           tCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGCGCTTAAATGTttgccccaagcgc                                                                                                    <  4:185895/149‑1 (MQ=255)                                             cATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCg                                                                                                  <  4:274318/149‑1 (MQ=255)                                              aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                           <  3:25307/107‑4 (MQ=255)                                              aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                           <  3:29922/107‑4 (MQ=255)                                              aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                           >  4:29922/1‑104 (MQ=255)                                              aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAATGGCCGCCTTTTTCgtga                                                                                                                                           >  4:25307/1‑104 (MQ=255)                                              aTGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCgg                                                                                                 >  6:170021/1‑149 (MQ=255)                                                 cTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCtt                                                                                              <  5:170021/149‑1 (MQ=255)                                                    aGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGGGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCttt                                                                                             <  6:319378/147‑1 (MQ=255)                                                    aGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCttt                                                                                             <  5:247026/147‑1 (MQ=255)                                                    aGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTaa                                                                                           <  5:49242/149‑1 (MQ=255)                                                      ggCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGc                                                                                         <  4:125924/149‑1 (MQ=255)                                                      ggCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGc                                                                                         >  4:499/1‑149 (MQ=255)                                                        cACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTc                                                                                                                                               <  1:354813/93‑1 (MQ=255)                                                        cACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTc                                                                                                                                               >  2:354813/1‑93 (MQ=255)                                                                gCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTga                                                                                >  3:5670/1‑148 (MQ=255)                                                                gCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTga                                                                                <  4:5670/148‑1 (MQ=255)                                                                 cGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTg                                                                                 >  7:91604/1‑146 (MQ=255)                                                                  gTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGCCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTgaaga                                                                             <  1:112981/149‑1 (MQ=255)                                                                        aCGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCa                                                                         >  4:198716/1‑147 (MQ=255)                                                                         cGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAACCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCagag                                                                      >  7:98304/1‑149 (MQ=255)                                                                                               cTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAg                                                >  5:311841/1‑149 (MQ=255)                                                                                                 gATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCaaa                                                 >  3:260509/1‑146 (MQ=255)                                                                                                 gATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCaaa                                                 <  4:260509/146‑1 (MQ=255)                                                                                                  aTCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGc                                               <  3:499/147‑1 (MQ=255)                                                                                                   tctcCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCgtg                                            >  3:180307/1‑149 (MQ=255)                                                                                                       atAGGGTGCTTGAATAGGTGGTGTCAAGAAAAAGGCGGCCATTTTATCCCACTGGGGGCGATTGTGGCAAACATTTAATCGAAAAGCGCCGGCTTTAAGACGAAACTGAAGAAGCAGAGCGAGACTCGGGGAGACCCAAAGCgtgaga                                         >  5:74055/2‑145 (MQ=255)                                                                                                        tAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACt                                       >  3:20315/1‑149 (MQ=255)                                                                                                          gggggCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAACCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACt                                       <  8:98304/147‑1 (MQ=255)                                                                                                          gggggCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACtt                                      >  8:201126/1‑148 (MQ=255)                                                                                                                ttGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTaa                                                                                           <  4:96686/89‑1 (MQ=255)                                                                                                                ttGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTaa                                                                                           >  3:96686/1‑89 (MQ=255)                                                                                                                    aTAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCgg                           >  3:322603/1‑149 (MQ=255)                                                                                                                    aTAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCg                            <  5:183643/148‑1 (MQ=255)                                                                                                                     tAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGc                          <  4:322603/149‑1 (MQ=255)                                                                                                                        gtggtgTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTg                       >  4:346989/1‑149 (MQ=255)                                                                                                                        gtggtgTCACGAAAAAGGAGGCCATTTTCGCACACTGGGGGCGATTGGGGCAAACATTTAAGCGCCAAACGACGGATTTAAGCAGAAACTGAAGAAGCAGAGCGCGACACGGGGCGACCCAAACCGAGAAACATGAGGCATTCGGCtt                        >  2:139282/1‑148 (MQ=255)                                                                                                                         tggtgTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGc                      <  4:20315/149‑1 (MQ=255)                                                                                                                           gtgtCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCATAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCa                     <  1:139282/148‑1 (MQ=255)                                                                                                                            tgtCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGc                   <  3:346989/149‑1 (MQ=255)                                                                                                                            tgtCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGc                   <  6:133800/149‑1 (MQ=255)                                                                                                                              tCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCTGCTGTAAGCCGCAACTGAAGAAGCAGAGCGCGCCTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCt                  >  3:80695/1‑148 (MQ=255)                                                                                                                                aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAg                    <  6:42111/144‑1 (MQ=255)                                                                                                                                aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAg                    >  5:42111/1‑144 (MQ=255)                                                                                                                                aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCAt               >  7:163163/1‑149 (MQ=255)                                                                                                                                aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCAt               >  8:169436/1‑149 (MQ=255)                                                                                                                                aCGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCAt               >  1:163183/1‑149 (MQ=255)                                                                                                                                 acgaaaaGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATc              >  5:231978/4‑149 (MQ=255)                                                                                                                                       aGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCACATCGCCAat        <  4:180307/149‑1 (MQ=255)                                                                                                                                         gcggcCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAata       >  8:109185/1‑148 (MQ=255)                                                                                                                                              cATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCaaa  >  7:293677/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                   ATGGGGGCACCTGGGTGGCCGCTGTTGGCCTTTTGCACGGCATCCATGCTGAGGGCACGAATGGCGTTGGCACGTTCACGACGGCTGGGCATCGCTGATCTCCTAGGGGGCTTGAATAGGTGGTGTCACGAAAAAGGCGGCCATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAA  >  NC_002947/5657603‑5657891

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