Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,838,507	G→C	5.9%	R701P (CGG→CCG)	PP_3388 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,838,507	0	G	C	5.9%	212.0 / 3.7	68	R701P (CGG→CCG)	PP_3388	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base C (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CTGCTGGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGC  >  NC_002947/3838359‑3838652                                                                                                                                                     |                                                                                                                                                   ctgctgGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTcgcgcg                                                                                                                                                    <  7:246619/149‑1 (MQ=39)       gAGGGGATGCTGATGGCCTGCAAAGCGCGCGAATTGCCTGCCAGTGATGCGCGGGTGAAGACGGAACTTCTGGCGGCCGCGATGACTACGGCTGGGGCGGGCTTCGAAGTTGGGGCCAGGTAGGTGGAGGAGGTGGTCGCGCGGTACg‑g                                                                                                                                             >  4:46771/1‑149 (MQ=11)         ggggATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑gcg                                                                                                                                           <  3:46771/149‑1 (MQ=39)            gATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑Gcgcg                                                                                                                                         >  5:78838/1‑148 (MQ=39)              tGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCCGTACC‑GCGCGacc                                                                                                                                      <  5:118780/149‑3 (MQ=255)                 tGATGGCCTTCAAAGCGCGTGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAgt                                                                                                                                   >  4:162254/1‑149 (MQ=255)                  gATGGCCTTCAAAGCGCGTGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAgtg                                                                                                                                  >  4:261386/1‑149 (MQ=255)                  gATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAgtg                                                                                                                                  >  3:252652/1‑149 (MQ=39)                   aTGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAgt                                                                                                                                   <  8:101894/147‑1 (MQ=39)                    tGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGtgtg                                                                                                                                <  3:186459/149‑1 (MQ=39)                        cTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGTGTGAAGACCGATCTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCg                                                                                                                            <  8:249397/149‑1 (MQ=255)                        cTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCg                                                                                                                            >  1:48461/1‑149 (MQ=39)                        cTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACAGCTGCGGCGGGCTTCGAAGTTGGAGCTAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCg                                                                                                                            >  7:71434/1‑149 (MQ=255)                         ttCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTcgcgcg                                                                                                                                                    >  1:149986/1‑125 (MQ=40)                         ttCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTcgcgcg                                                                                                                                                    <  2:149986/125‑1 (MQ=40)                           cAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCgg                                                                                                                           <  4:252652/147‑1 (MQ=39)                           cAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGAc                                                                                                                         >  8:134278/1‑149 (MQ=39)                           cAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGAc                                                                                                                         >  8:39259/1‑149 (MQ=39)                              aGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACgag                                                                                                                      <  2:78874/149‑1 (MQ=39)                                 gcgcgAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACgaggag                                                                                                                   <  5:251869/149‑1 (MQ=39)                                   gtgAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGagca                                                                                                                 <  3:162254/147‑1 (MQ=39)                                           ccTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGc                                                                                                         >  2:294233/1‑149 (MQ=37)                                            cTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCt                                                                                                        >  4:192658/1‑149 (MQ=37)                                             tGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCtt                                                                                                       >  1:328095/1‑149 (MQ=38)                                              gACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCCGTACC‑GCGCGacc                                                                                                                                      >  4:96980/1‑115 (MQ=255)                                              gACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCCGTACC‑GCGCGacc                                                                                                                                      <  3:96980/117‑3 (MQ=255)                                                cAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAg                                                                                                    <  4:120724/149‑1 (MQ=39)                                                cAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACAGCTGCGGCGGGCTTCGAAGTTGGAGCTAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAg                                                                                                    <  8:71434/149‑1 (MQ=255)                                                   tGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGc                                                                                                 >  4:84715/1‑149 (MQ=39)                                                     aTGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑gc                                                                                                                                            <  2:46322/104‑1 (MQ=40)                                                     aTGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑gc                                                                                                                                            >  1:46322/1‑104 (MQ=40)                                                       gcgcGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGttt                                                                                             >  3:69151/1‑149 (MQ=39)                                                         gcgGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGa                                                                                           <  7:39259/149‑1 (MQ=39)                                                          cgGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGa                                                                                           >  6:7616/1‑148 (MQ=39)                                                                aaGACCGAACTTCTGGCGGCCGCGATTACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTAt                                                                                     <  1:294233/148‑1 (MQ=255)                                                                 aGACCGAACTTCTGGCGGCCGAGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATg                                                                                    <  7:291561/148‑1 (MQ=255)                                                                  gACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGc                                                                                                            >  8:126737/1‑123 (MQ=38)                                                                  gACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGc                                                                                                            <  7:126737/123‑1 (MQ=38)                                                                     cGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGgcg                                                                               >  7:31123/1‑149 (MQ=39)                                                                     cGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGACGCTTTAGGGCTTTTGAACCTATGGGgcg                                                                               >  8:206349/1‑149 (MQ=255)                                                                               gcggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGATCCAGTTACGTGGATCATGTGGTCGCTCGGTACG‑GCGCGAACAGTGTGACAGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGc                                                                     <  7:206349/149‑1 (MQ=255)                                                                               gcggcCGCGATGACTACGGCTGCGGCGGGCTTCAAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGc                                                                     >  5:26127/1‑149 (MQ=37)                                                                                cggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCt                                                                    <  6:78838/149‑1 (MQ=37)                                                                                 ggcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCtg                                                                   <  5:7616/149‑1 (MQ=37)                                                                                  gcCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATAAGGTGGTCGCGCGGTACG‑GAGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCtgt                                                                  >  1:98054/1‑149 (MQ=255)                                                                                   ccGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCtgt                                                                  >  6:248192/1‑148 (MQ=38)                                                                                    cgcgATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGa                                                                >  7:273936/1‑149 (MQ=39)                                                                                         tGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAggt                                                           >  6:121265/1‑149 (MQ=255)                                                                                            cTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAggtggt                                                        <  5:121265/149‑1 (MQ=255)                                                                                               cGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTc                                                     <  3:192658/149‑1 (MQ=255)                                                                                                ggCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCg                                                    <  2:48461/149‑1 (MQ=255)                                                                                                      ggcggGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACGTGCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTc                                               >  1:87138/1‑149 (MQ=255)                                                                                                         gggCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCt                                           <  6:166088/149‑1 (MQ=255)                                                                                                            cTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTGGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggt                                        <  3:84715/149‑1 (MQ=255)                                                                                                            cTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAggt                                                           <  1:51106/130‑1 (MQ=255)                                                                                                            cTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAggt                                                           >  2:51106/1‑130 (MQ=255)                                                                                                             ttCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggtg                                       >  7:109854/1‑149 (MQ=255)                                                                                                               cGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGa                                                                >  5:305402/1‑122 (MQ=40)                                                                                                               cGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGa                                                                <  6:305402/122‑1 (MQ=40)                                                                                                               cGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCtggtg                                       <  2:249402/147‑1 (MQ=255)                                                                                                                gAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTggg                                     <  2:328095/148‑1 (MQ=255)                                                                                                                      ggAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGc                              >  7:70628/1‑149 (MQ=255)                                                                                                                          ccAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGAt                          <  3:261386/149‑1 (MQ=255)                                                                                                                            aGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GAGCGAACAGTGTGCCCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTATGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGGTTCGTACTCGCGTGGTGGGACTTCGcctggat                        >  7:212083/1‑144 (MQ=17)                                                                                                                             gTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCg                       >  4:78890/1‑149 (MQ=255)                                                                                                                                  gTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTca                  <  8:273936/149‑1 (MQ=255)                                                                                                                                      aTCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCAATGCGGCTcagcag              <  6:26127/149‑1 (MQ=255)                                                                                                                                         aGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGTCTAGGAGCAGCGGCAGCTTTAGTGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGTTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTa           <  2:98054/149‑1 (MQ=21)                                                                                                                                           gtggtCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTaa          >  8:69976/1‑148 (MQ=255)                                                                                                                                           gtggtCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGACTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTaa          >  8:74288/1‑148 (MQ=255)                                                                                                                                            tgttCGCGCCGTACC‑GCGCGAACAGTGTGACCGGACGAGGAGAAGCGGCCGCTATAGGGCGTTTGAAACTATGGGGAGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTAGCCTGGTGGGACTTCGCCGATTCGTCTCAGCAGTTAata        >  6:178408/4‑149 (MQ=18)                                                                                                                                                  cgcgGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGc  <  8:31123/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CTGCTGGAGGGGATGCTGATGGCCTTCAAAGCGCGCGAATTGCCTGACAGTGATGCGCGGGTGAAGACCGAACTTCTGGCGGCCGCGATGACTACGGCTGCGGCGGGCTTCGAAGTTGGAGCCAGTTACGTGGATCAGGTGGTCGCGCGGTACG‑GCGCGAACAGTGTGACCGGACGAGGAGCAGCGGCCGCTTTAGGGCGTTTGAAACTATGGGGCGCCAGTTTAGCTGTGACAGGTGGTTTCGTACTCGCCTGGTGGGACTTCGCCGATGCGGCTCAGCAGTTAATATACCGC  >  NC_002947/3838359‑3838652

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