Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,827,679	T→C	6.8%	W198R (TGG→CGG)	mtgA →	monofunctional biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,827,679	0	T	C	6.8%	165.8 / 2.4	59	W198R (TGG→CGG)	mtgA	monofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base C (2/2);  total (29/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 = 7.01e-01

TCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGG  >  NC_002947/5827533‑5827821                                                                                                                                                    |                                                                                                                                               tcGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTg                                                                                                                                               >  6:42445/1‑148 (MQ=255) tcGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTg                                                                                                                                               >  4:221310/1‑148 (MQ=255)   gAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCTGGGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGACGCCAGCATGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGa                                                                                                                                             <  8:17023/148‑1 (MQ=255)   gAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTg                                                                                                                                               <  4:90966/146‑1 (MQ=255)               aaCAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGccgacc                                                                                                                                <  2:240073/149‑1 (MQ=255)                aCAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGccgacc                                                                                                                                <  3:221310/148‑1 (MQ=255)                 cAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGccgacatg                                                                                                                              >  2:126281/1‑146 (MQ=255)                   gcgcCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CACTCAAGTGGAGTGCCAGCCGACCGAGt                                                                                                                            <  1:8589/149‑1 (MQ=255)                               gCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGCGGAGTGCCAGCCGACCGAGTGCCTACGTGGc                                                                                                                 >  4:139671/1‑148 (MQ=255)                               gCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGCGGAGTGCCAGCCGACCGAGTGCCTACGTGGc                                                                                                                 >  4:134762/1‑148 (MQ=255)                                  aGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgc                                                                                                             >  1:220555/1‑149 (MQ=255)                                   gggCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGTCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGTGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgcc                                                                                                            <  6:249801/149‑1 (MQ=255)                                     gCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgccgg                                                                                                          <  2:178825/149‑1 (MQ=255)                                               tGCCCATGCAGCGGCGCGTTATCAATTTGGCGTCGACGCCATCCGGCTCAGCCGCCATCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAATCCTCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTg                                                                                                  <  1:93831/148‑1 (MQ=255)                                                   cAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGCGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCg                                                                                            <  3:134762/149‑1 (MQ=255)                                                   cAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGCGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCg                                                                                            <  3:139671/149‑1 (MQ=255)                                                    aGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGc                                                                                           >  6:172482/1‑149 (MQ=255)                                                    aGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGc                                                                                           <  3:105818/149‑1 (MQ=255)                                                       cAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCgg                                                                                        <  5:42445/149‑1 (MQ=255)                                                       cAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCgg                                                                                        <  6:246370/149‑1 (MQ=255)                                                            gcgcGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGAt                                                                                   >  5:116392/1‑149 (MQ=255)                                                            gcgcGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGAt                                                                                   >  6:45436/1‑149 (MQ=255)                                                                gTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGc                                                                               <  2:220555/149‑1 (MQ=255)                                                                 ttATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGc                                                                               >  3:241989/1‑148 (MQ=255)                                                                  tATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGc                                                                               >  3:165517/1‑147 (MQ=255)                                                                  tATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGc                                                                               >  1:224443/1‑147 (MQ=255)                                                                     cACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCg                                                                                         >  2:290497/1‑134 (MQ=255)                                                                     cACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCg                                                                                         <  1:290497/134‑1 (MQ=255)                                                                        tttGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTggg                                                                       >  7:109102/1‑149 (MQ=255)                                                                        tttGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTggg                                                                       >  2:148976/1‑149 (MQ=255)                                                                        tatGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTggg                                                                       >  6:338084/3‑149 (MQ=255)                                                                          tGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggc                                                                      <  5:172482/148‑1 (MQ=255)                                                                                 gACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAg                                                              <  5:45436/149‑1 (MQ=255)                                                                                       aGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCt                                                        >  4:213457/1‑149 (MQ=255)                                                                                          cGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGAt                                                     <  2:224443/149‑1 (MQ=255)                                                                                            gCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATg                                                    >  8:130925/1‑148 (MQ=255)                                                                                             cTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATg                                                    >  2:180712/1‑147 (MQ=255)                                                                                                         caggcagCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActcctc                                      <  3:248168/149‑1 (MQ=255)                                                                                                         caggcagCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActcctc                                      >  4:248168/1‑149 (MQ=255)                                                                                                         caggcagCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActcctc                                      <  1:148976/149‑1 (MQ=255)                                                                                                          aggcagCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                                     <  2:202085/149‑1 (MQ=255)                                                                                                              agcgcagcTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGa                                  >  3:302598/1‑148 (MQ=255)                                                                                                                  cagcTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCt                             >  2:119365/1‑149 (MQ=255)                                                                                                                   agcTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCtt                            <  7:130925/149‑1 (MQ=255)                                                                                                                    gcTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                           <  3:213457/149‑1 (MQ=255)                                                                                                                            cGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAg                   <  8:109102/149‑1 (MQ=255)                                                                                                                                   gcCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                           >  7:286896/1‑134 (MQ=255)                                                                                                                                   gcCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                           <  8:286896/134‑1 (MQ=255)                                                                                                                                    ccAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccggcc           <  4:302598/149‑1 (MQ=255)                                                                                                                                         cc‑cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTc      >  4:66000/1‑149 (MQ=255)                                                                                                                                          c‑cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcg     >  5:60622/1‑149 (MQ=255)                                                                                                                                          c‑cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcg     >  5:186898/1‑149 (MQ=255)                                                                                                                                             agtCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGAAGCAGCAGGGCTCCTCCAGAACGCGTTGAATGCAGGCCCGGCGCTTTcgcg   >  4:99984/3‑149 (MQ=255)                                                                                                                                             aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcgcg   >  4:223593/1‑149 (MQ=255)                                                                                                                                             aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcgcg   <  1:180712/149‑1 (MQ=255)                                                                                                                                             aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAG‑AAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg  <  4:241989/149‑1 (MQ=255)                                                                                                                                              aTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccgg             <  4:11956/138‑1 (MQ=255)                                                                                                                                              aTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccgg             >  3:11956/1‑138 (MQ=255)                                                                                                                                                   gTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCtt                            <  2:274405/118‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               TCGAGGTCTATTTGAACAGCGCCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCC‑CAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGG  >  NC_002947/5827533‑5827821

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