Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,418,251	A→C	6.4%	*237S (TAA→TCA)	purC →	phosphoribosylaminoimidazole‑succinocarboxamide synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,418,251	0	A	C	6.4%	175.3 / 4.1	63	*237S (TAA→TCA)	purC	phosphoribosylaminoimidazole‑succinocarboxamide synthase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (2/2);  total (32/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 = 7.26e-01

GGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGCAA  >  NC_002947/1418109‑1418398                                                                                                                                               |                                                                                                                                                     ggCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCg                                                                                                                                                  <  2:145563/147‑1 (MQ=255)               cccGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAgc                                                                                                                                  <  2:113347/149‑1 (MQ=255)                 cGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAgcgc                                                                                                                                <  4:344973/149‑1 (MQ=255)                 cGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAgcgc                                                                                                                                >  7:254698/1‑149 (MQ=255)                  ggACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGcc                                                                                                                               >  6:139985/1‑149 (MQ=255)                        cTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATAc                                                                                                                         >  8:294873/1‑149 (MQ=255)                         tGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATAc                                                                                                                         >  7:209014/1‑148 (MQ=255)                                      aCAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGAttttttt                                                                                                           >  8:249888/1‑149 (MQ=255)                                        aaGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTg                                                                                                          >  2:259308/1‑148 (MQ=255)                                        aaGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTg                                                                                                          >  4:138170/1‑148 (MQ=255)                                          ggAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCGGTTAGCGGCACGTTAACGCAAGCGCCTGATAACACGCGATTTTTTTGaa                                                                                                        >  5:48536/1‑148 (MQ=255)                                                ccGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGa                                                                                                <  7:294873/149‑1 (MQ=255)                                                  gCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCAGAAACCACGCGATTTTTTTGAAAAAA‑AGAtt                                                                                              >  2:138199/1‑149 (MQ=255)                                                        agaTGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGtt                                                                                        <  6:227292/149‑1 (MQ=255)                                                         gaTGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAg                                                                                                                                                    <  8:11473/89‑1 (MQ=255)                                                          aTGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGttt                                                                                       >  3:347424/1‑148 (MQ=255)                                                              aCAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTaa                                                                                    >  6:71312/1‑147 (MQ=255)                                                               cAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTaa                                                                                    >  4:347746/1‑146 (MQ=255)                                                                aaGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCGGGAAGCGGCCCGTTCACCCCAGCGCCGGATACCCCGCGATTTTTTTTAAAAAA‑AGATTTGCGTTTTaaaat                                                                                 >  1:351751/1‑147 (MQ=255)                                                                     ccGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCt                                                                           <  4:347424/149‑1 (MQ=255)                                                                             gCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTaaaa                                                                                  >  3:150134/1‑134 (MQ=255)                                                                             gCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTaaaa                                                                                  >  1:154689/1‑134 (MQ=255)                                                                             gCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTaaaa                                                                                  <  4:150134/134‑1 (MQ=255)                                                                             gCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTaaaa                                                                                  <  2:154689/134‑1 (MQ=255)                                                                                      tGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCAt                                                          >  8:236950/1‑149 (MQ=255)                                                                                        ggCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTg                                                        <  7:249888/149‑1 (MQ=255)                                                                                           gACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTg                                                        >  6:2206/1‑146 (MQ=255)                                                                                            aCGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGgaga                                                    <  5:265709/149‑1 (MQ=255)                                                                                             cGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAAGAGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGgaga                                                    <  1:259308/149‑1 (MQ=255)                                                                                             cGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGgagag                                                   <  7:130459/149‑1 (MQ=255)                                                                                                    gAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTACCGCAAGCGCCTGGTACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTTTGTTGAGCGGCATTGGATCGATGCggg                                            >  4:202469/1‑149 (MQ=255)                                                                                                     aaGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGa                                           <  7:236950/149‑1 (MQ=255)                                                                                                        ccTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTg                                        <  3:347746/149‑1 (MQ=255)                                                                                                          tACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGt                                      >  7:165273/1‑149 (MQ=255)                                                                                                            cGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTc                                     <  8:209014/148‑1 (MQ=255)                                                                                                              aagaagTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGaa                                  >  3:273802/1‑149 (MQ=255)                                                                                                                gaagTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACg                                <  4:273802/149‑1 (MQ=255)                                                                                                                  agTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACggg                              <  8:254698/149‑1 (MQ=255)                                                                                                                  agTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACgg                               <  5:2206/148‑1 (MQ=255)                                                                                                                   gTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGa                             >  1:196839/1‑149 (MQ=255)                                                                                                                   gTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTAACGCAAGCGCCGGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGa                             >  6:308288/1‑149 (MQ=255)                                                                                                                      gCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATgcgc                                                               <  7:219462/112‑1 (MQ=255)                                                                                                                       ccAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTc                                     <  2:220727/137‑1 (MQ=255)                                                                                                                       ccAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTc                                     >  1:220727/1‑137 (MQ=255)                                                                                                                        cAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGAt                        >  3:308168/1‑149 (MQ=255)                                                                                                                         aaGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGAtt                       >  8:21672/1‑149 (MQ=255)                                                                                                                           gcgcCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTc                      >  8:175585/1‑148 (MQ=255)                                                                                                                               cTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTTAAAGAACGCTGGTATGATGCGCGGGATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGaaa                  >  5:228319/1‑148 (MQ=255)                                                                                                                               cTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAAt                 >  3:107840/1‑149 (MQ=255)                                                                                                                                  ggCGTGCCTCTTTAAGCGGCACGTTCACGCAAGCGCCTGATACCAGGCGATTTTTTTGAAAAAA‑ATATTTGCGTTTTAAAAGAACGCTGGTATAATGCGCGGCATTGGAGAGATGCCGGAGGGGTCGAACGGGACGGATTCGAAATCCg              <  1:138199/149‑1 (MQ=255)                                                                                                                                  ggCGTGCCGCTGTACGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATc                <  2:351751/147‑1 (MQ=255)                                                                                                                                  ggCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATcc               <  8:165273/148‑1 (MQ=255)                                                                                                                                  ggCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATc                <  5:71312/147‑1 (MQ=255)                                                                                                                                    cGTGCCGCTTACAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGt             <  2:20964/148‑1 (MQ=255)                                                                                                                                    cGTGCCGCTTACAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGt             >  1:20964/1‑148 (MQ=255)                                                                                                                                    cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  1:139534/149‑1 (MQ=255)                                                                                                                                    cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  3:138170/149‑1 (MQ=255)                                                                                                                                    cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  2:196839/149‑1 (MQ=255)                                                                                                                                    cGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  7:21672/149‑1 (MQ=255)                                                                                                                                     gTGCCGCTTACAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            <  1:310709/148‑1 (MQ=255)                                                                                                                                     gTGCCGCTTACAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGtt            >  2:310709/1‑148 (MQ=255)                                                                                                                                      tGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTg           >  3:31355/1‑148 (MQ=255)                                                                                                                                              tAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGCaa  >  1:180802/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     GGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCCGCAAGAAGATGGACAAGGACCGCTTCCGCCAGGGCCTGGGCGACGTGATCGAAGCCTACGAAGAAGTTGCCAAGCGCCTGGGCGTGCCGCTGTAAGCGGCACGTTCACGCAAGCGCCTGATACCACGCGATTTTTTTGAAAAAA‑AGATTTGCGTTTTAAAAGAACGCTGGTATGATGCGCGGCATTGGAGAGATGCCGGAGTGGTCGAACGGGACGGATTCGAAATCCGTTGTACTGGCAA  >  NC_002947/1418109‑1418398

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