Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,316,912	G→C	6.2%	intergenic (‑115/+58)	ribC ← / ← treSA	riboflavin synthase/trehalose synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,316,912	0	G	C	6.2%	183.4 / 2.0	60	intergenic (‑115/+58)	ribC/treSA	riboflavin synthase/trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base C (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

GGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCT  >  NC_002947/3316770‑3317049                                                                                                                                                |                                                                                                                                          ggACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAg                                                                                                                                     <  8:150039/149‑1 (MQ=255)   aCGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAggc                                                                                                                                   <  2:306643/149‑1 (MQ=255)    cGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAggcg                                                                                                                                  <  7:203702/149‑1 (MQ=255)       ccTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCCCCGCGCAggcgcca                                                                                                                               <  2:298382/149‑4 (MQ=255)       ccTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCCCCGCGCAggcgcca                                                                                                                               <  2:293287/149‑4 (MQ=255)         tGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTg                                                                                                                             >  5:266295/1‑149 (MQ=255)          ggACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGc                                                                                                                            <  7:32030/149‑1 (MQ=255)             cGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCtt                                                                                                                          <  1:139815/148‑1 (MQ=255)               aTGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTACGAACGCTATGCCCGAAGCGCAGTAA‑ATGCAAGACCAGGCATCCTGGATGGCAACCCGGCTTGCGAGTCTGGAGCCTGCGCCGGGCAGGCGGCGTGCTtctg                                                                                                                       >  7:145350/1‑149 (MQ=255)                  ccGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGc                                                                                                                    <  1:95137/149‑1 (MQ=255)                          cATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCgg                                                                                                            >  2:303625/1‑149 (MQ=255)                          cATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGCCGGCGTGCTTGAGGGAGAAACCgg                                                                                                            >  5:390833/1‑149 (MQ=255)                           aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGt                                                                                                           <  6:266295/149‑1 (MQ=255)                               aGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGc                                                                                                       <  1:303625/149‑1 (MQ=255)                                gTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGc                                                                                                       >  2:138116/1‑148 (MQ=255)                                   gCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACg                                                                                                   >  5:117780/1‑149 (MQ=255)                                    cTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACgg                                                                                                  >  5:108026/1‑149 (MQ=255)                                       cTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCa                                                                                                <  6:108026/148‑1 (MQ=255)                                             aCGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTtgt                                                                                                                        <  2:198750/118‑1 (MQ=255)                                             aCGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTtgt                                                                                                                        >  1:198750/1‑118 (MQ=255)                                             aCGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCAc                                                                                         >  7:154248/1‑149 (MQ=255)                                               ggAAAAGAAACCTTCCGAACGCTATGGCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCATGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAGCCGGTTTGCGACGGCATGGCCACGc                                                                                       >  8:211347/1‑149 (MQ=255)                                                  aaagaaaaCTTACGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGGTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAg                                                                                    <  1:124621/149‑1 (MQ=255)                                                    agaaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGAt                                                                                  <  7:25788/149‑1 (MQ=255)                                                     gaaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGAtt                                                                                 >  6:51920/1‑149 (MQ=255)                                                       aaaCTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACg                                                                                                   <  8:196581/129‑1 (MQ=255)                                                       aaaCTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACg                                                                                                   >  7:196581/1‑129 (MQ=255)                                                          cTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAAc                                                                            <  8:366227/149‑1 (MQ=255)                                                           ttCCGAACGCTATGCCCGTAGCGCAGTAAGATCCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGCCCCCTGCGCGGGGCAGGCGCCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                             >  1:29455/1‑148 (MQ=255)                                                               gAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGcc                                                                       >  1:191590/1‑149 (MQ=255)                                                                     tATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCgctgc                                                                 >  6:154835/1‑149 (MQ=255)                                                                         cccGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGtt                                                             <  1:138116/149‑1 (MQ=255)                                                                          ccGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGtt                                                             >  1:379322/1‑148 (MQ=255)                                                                                 gcAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATcc                                                     <  6:267923/149‑1 (MQ=255)                                                                                 gcAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATcc                                                     <  8:140545/149‑1 (MQ=255)                                                                                     tAA‑ATGAAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAg                                                 >  6:221993/1‑149 (MQ=255)                                                                                           gCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGca                                              >  3:99490/1‑147 (MQ=255)                                                                                            cAAGACCAGGCTTCCTGGCTTGCATGCGGGCTTGCGAGTCTGGCGCCTGCGCGGGGCATTCGGCGTGATTGTGGTCGAAACCGGTTTGCGACGGCAttg                                                                                             <  7:70020/99‑3 (MQ=255)                                                                                            cAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAg                                           >  1:11497/1‑149 (MQ=255)                                                                                             aaGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAgc                                          >  3:11499/1‑149 (MQ=255)                                                                                               gACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAgc                                          <  8:340221/147‑1 (MQ=255)                                                                                               gACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAgc                                          >  7:340221/1‑147 (MQ=255)                                                                                                 ccAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCt                                      <  2:379322/149‑1 (MQ=255)                                                                                                     gCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTa                                  <  6:117780/149‑1 (MQ=255)                                                                                                     gCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTa                                  >  5:181274/1‑149 (MQ=255)                                                                                                     gCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTa                                  <  5:221993/149‑1 (MQ=255)                                                                                                       ttCCTGGCGGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                           <  8:252975/122‑1 (MQ=255)                                                                                                          ctggctggCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGt                             <  2:11497/149‑1 (MQ=255)                                                                                                             gctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCa                          <  6:181274/149‑1 (MQ=255)                                                                                                               tggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTc                                     >  4:337673/1‑136 (MQ=255)                                                                                                               tggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTc                                     <  3:337673/136‑1 (MQ=255)                                                                                                                ggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGt                       >  3:145660/1‑149 (MQ=255)                                                                                                                    aGCCGGCTTGCGAGTCTGGCGCCTGGGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGGATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                   <  8:145350/149‑1 (MQ=255)                                                                                                                      ccGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGAGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAcc                  >  5:83904/1‑148 (MQ=255)                                                                                                                        ggCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCg                                    <  1:192255/128‑1 (MQ=255)                                                                                                                        ggCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCg                                    >  2:192255/1‑128 (MQ=255)                                                                                                                           ttGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAgg            >  5:502/1‑149 (MQ=255)                                                                                                                                  cacgccgcCTGCCCCGCGCAGGCGCCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGt     >  7:196921/6‑149 (MQ=255)                                                                                                                                     ggCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGTGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCt  >  5:203181/1‑149 (MQ=255)                                                                                                                                      gcgcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTcc   <  5:51920/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                          GGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAA‑ATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCT  >  NC_002947/3316770‑3317049

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