Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,263,552	C→A	6.6%	A115A (GCC→GCA)	trpF →	N‑(5'‑phosphoribosyl)anthranilate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,263,552	0	C	A	6.6%	184.8 / 3.4	61	A115A (GCC→GCA)	trpF	N‑(5'‑phosphoribosyl)anthranilate isomerase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base A (2/2);  total (32/29)
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.65e-01

GCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCTGTCGGCCGAT  >  NC_002947/2263405‑2263699                                                                                                                                                    |                                                                                                                                                      gcgcTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCt                                                                                                                                                      >  8:299781/1‑149 (MQ=255)       cGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑gg                                                                                                                                               <  2:37492/148‑1 (MQ=255)                   gaTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCt                                                                                                                                                      <  3:275394/131‑1 (MQ=255)                   gaTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCt                                                                                                                                                      >  4:275394/1‑131 (MQ=255)                    aTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGt                                                                                                                                   <  4:282238/147‑1 (MQ=255)                    aTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGt                                                                                                                                   >  3:282238/1‑147 (MQ=255)                    aTCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGc                                                                                                                                 >  1:119551/1‑149 (MQ=255)                      ccTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGccc                                                                                                                               >  5:189025/1‑149 (MQ=255)                      ccTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCATGTC‑‑GGCTTTACGCGGGTGccc                                                                                                                               >  3:40834/1‑149 (MQ=255)                      ccTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCATGTC‑‑GGCTTTACGCGGGTGccc                                                                                                                               >  5:91007/1‑149 (MQ=255)                       cTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCcg                                                                                                                              >  4:327927/1‑149 (MQ=255)                          gAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCcgcg                                                                                                                            <  6:189025/148‑1 (MQ=255)                                   ccGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGc                                                                                                                   <  2:119551/148‑1 (MQ=255)                                   ccGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGc                                                                                                                   <  7:299781/148‑1 (MQ=255)                                              gCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATg                                                                                                       >  1:6221/1‑149 (MQ=255)                                               cTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTAt                                                                                                        <  1:260150/147‑1 (MQ=255)                                               cTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTAt                                                                                                        <  3:327927/147‑1 (MQ=255)                                               cTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATGTTGCTCGACACCTATg                                                                                                       >  8:114077/1‑148 (MQ=255)                                                       tCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCg                                                                                              >  2:125397/1‑149 (MQ=255)                                                          cGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGc                                                                                           >  7:94114/1‑149 (MQ=255)                                                                tGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCATGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAgggg                                                                                     <  4:40834/149‑1 (MQ=255)                                                                tGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCATGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAgggg                                                                                     <  6:91007/149‑1 (MQ=255)                                                                   aaCCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGCCACCTATGTACCGGGCGTGCCAGGGGGaa                                                                                  >  8:173643/1‑149 (MQ=255)                                                                          cAGGACTGTGACGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGa                                                                           <  8:267850/149‑1 (MQ=255)                                                                                  tgAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTc                                                                    >  4:152543/1‑148 (MQ=255)                                                                                  tgAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTc                                                                    >  8:5617/1‑148 (MQ=255)                                                                                    aaGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGAc                                                                 <  1:12453/149‑1 (MQ=255)                                                                                    aaGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGAc                                                                 >  3:150509/1‑149 (MQ=255)                                                                                    aaGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGGGCCAGGGGGAACCGGCGAGGCGTTCGAc                                                                 <  1:167610/149‑1 (MQ=255)                                                                                    aaGGCTACCACCGCCCCTGGATCAAAGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGAc                                                                 >  1:53734/1‑149 (MQ=255)                                                                                     aGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACt                                                                >  2:151316/1‑149 (MQ=255)                                                                                         tACCACCGACGCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGAGCTGGAGGCGGCCTGTC‑‑GGGTTTACGAGGGTGCCAGCGGCATTTTGGTCGAGACCTATTTACCGGGCGTGCCAGGGGGAACCGGCGAGGGGGTCGACTGGTc                                                            >  4:206412/1‑149 (MQ=255)                                                                                           ccaccGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAggg                                                                                      <  1:139642/121‑1 (MQ=255)                                                                                           ccaccGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAggg                                                                                      >  2:139642/1‑121 (MQ=255)                                                                                             accGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTCAGGGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCAt                                                          >  1:271980/1‑149 (MQ=255)                                                                                             accGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTg                                                        >  8:202861/1‑149 (MQ=255)                                                                                             accGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTg                                                        <  7:202861/149‑1 (MQ=255)                                                                                              ccGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTgg                                                       >  2:159744/1‑149 (MQ=255)                                                                                               cGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCt                                                                                                                  >  5:118127/1‑89 (MQ=255)                                                                                               cGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCt                                                                                                                  <  6:118127/89‑1 (MQ=255)                                                                                                 ccccTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTAc                                                    <  8:191558/149‑1 (MQ=255)                                                                                                 ccccTGGATCAAGGCACTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTAc                                                    <  8:94114/149‑1 (MQ=255)                                                                                                  cccTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTAcc                                                   >  1:31948/1‑149 (MQ=255)                                                                                                   ccTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCg                                                  <  4:144605/149‑1 (MQ=255)                                                                                                       gATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGgcgc                                              >  4:208399/1‑149 (MQ=255)                                                                                                              gCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTgg                                       <  3:208399/149‑1 (MQ=255)                                                                                                               cccTGCGTTTGCGCCCTGGCGATGACCTGGAGGCGGCCTTTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGTGAGCGCCAGGGGGAACCGGCGAGGAGATCGACTGGTCATTGGTACCGGCGCGTCTTGGc                                      <  3:206412/149‑1 (MQ=255)                                                                                                               cccTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGc                                      <  1:125397/149‑1 (MQ=255)                                                                                                                 cTGGGTGTGCGCCCTGGCTTTTACCTGGAGGCGGCCTGTC‑‑GGCTTTACGGGGGTGCCCGGGTCATTTTGCTCTACACATATTTACCGGGCGGGCCTGGGGGAACCGGCTAGGCGTTCGACTGTTCATTGTTACCGGCGGGTCTTGGCaa                                    <  2:53734/149‑1 (MQ=255)                                                                                                                    cGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGc                                  <  1:159744/148‑1 (MQ=255)                                                                                                                       gtgCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCAt                                                          <  2:312060/121‑1 (MQ=255)                                                                                                                       gtgCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCAt                                                          >  1:312060/1‑121 (MQ=255)                                                                                                                         gcgcCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAggg                                                                                      >  2:90376/1‑91 (MQ=255)                                                                                                                         gcgcCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAggg                                                                                      <  1:90376/91‑1 (MQ=255)                                                                                                                                 gCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTgc                     <  2:31948/148‑1 (MQ=255)                                                                                                                                  cGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTgcgg                   <  3:152543/149‑1 (MQ=255)                                                                                                                                         cTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGATTAGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCt            >  1:283141/1‑149 (MQ=255)                                                                                                                                         cTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCt            >  2:177884/1‑149 (MQ=255)                                                                                                                                           ggAGGCGGACTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCTGt          >  6:46849/1‑149 (MQ=255)                                                                                                                                                cggcCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCTGTCGGcc     <  2:271980/149‑1 (MQ=255)                                                                                                                                                  gcCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCAATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACAGGTCATTGGTACCGGCGCGTCATGGCAACCCCATCATCCATCCCGGCGGGCAGGCGGccct   >  2:85298/1‑147 (MQ=255)                                                                                                                                                    cTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCTGTCGGCCGAt  >  5:133386/1‑148 (MQ=255)                                                                                                                                                    cTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTACGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCTGTCGGCCGAt  >  2:2296/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCGCTGCGAGCTGAACGAGATCCTCGAAGTGGTCCCGCTGGACCTGCTACAGTTTCACGGCGATGAAACCCCGCAGGACTGTGAAGGCTACCACCGCCCCTGGATCAAGGCCCTGCGTGTGCGCCCTGGCGATGACCTGGAGGCGGCCTGTC‑‑GGCTTTACGCGGGTGCCCGCGGCATTTTGCTCGACACCTATGTACCGGGCGTGCCAGGGGGAACCGGCGAGGCGTTCGACTGGTCATTGGTACCGGCGCGTCTTGGCAAGCCGATCATCCTTGCGGGCGGGCTGTCGGCCGAT  >  NC_002947/2263405‑2263699

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