Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	313,966	A→T	11.1%	V192E (GTG→GAG) ‡	yrfG ←	purine nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	313,966	0	A	T	11.1%	175.2 / 11.9	72	V192E (GTG→GAG) ‡	yrfG	purine nucleotidase
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base T (5/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

CTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATA  >  NC_002947/313820‑314111                                                                                                                                                   |                                                                                                                                                  ctgactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCAc                                                                                                                                                  <  6:312386/148‑1 (MQ=255) ctgactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGAGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCAcc                                                                                                                                                 >  3:142085/1‑149 (MQ=255)  tgactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCg                                                                                                                                                >  7:173608/1‑149 (MQ=255)   gactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGc                                                                                                                                               >  4:228374/1‑149 (MQ=255)     ctgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCa                                                                                                                                             <  8:9775/149‑1 (MQ=255)     ctgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCa                                                                                                                                             <  8:173608/149‑1 (MQ=255)         tCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGGTGCGCCAGccg                                                                                                                                         <  1:20760/149‑3 (MQ=255)               aaGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGccc                                                                                                                                   <  4:212317/149‑1 (MQ=255)                 gTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCca                                                                                                                                 >  5:165013/1‑149 (MQ=255)                 gTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCCGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGcccac                                                                                                                                 >  8:180510/1‑147 (MQ=255)                      tACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACAcc                                                                                                                             >  4:311727/1‑148 (MQ=255)                         aGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGGCAGc                                                                                                                                           >  5:324447/1‑131 (MQ=255)                         aGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGc                                                                                                                                           <  6:324447/131‑1 (MQ=255)                           gCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAATAACTCCCGGTAATCCTCGACCGCCGCAACCTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAccg                                                                                                                       >  1:262729/1‑149 (MQ=255)                             ctctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGCCCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAccgcc                                                                                                                     >  8:144018/1‑149 (MQ=255)                             ctctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAccgcc                                                                                                                     <  3:298657/149‑1 (MQ=255)                               ctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGcccc                                                                                                                                  >  3:103153/1‑134 (MQ=255)                               ctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGcccc                                                                                                                                  <  4:103153/134‑1 (MQ=255)                               ctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAccgcc                                                                                                                     >  5:194490/1‑147 (MQ=255)                               ctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAccgcc                                                                                                                     <  6:194490/147‑1 (MQ=255)                                ttCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGTCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGt                                                                                                                   <  8:336290/148‑1 (MQ=255)                                ttCGCGGGCGCGCCCCCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGCCCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTg                                                                                                                  >  5:141480/1‑149 (MQ=255)                                       gcgcgcCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTg                                                                                                                  <  5:7448/142‑1 (MQ=255)                                       gcgcgcCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTg                                                                                                                  >  6:7448/1‑142 (MQ=255)                                       gcccgaccGCGAATGGCCCCTTAAAGGCCTACCACTAAATCCCGGTAATCAACCACCGCCGCAAACTCTTCAGTATCCAGCGGACCTTCTTAGCTAAACGGCTCACGCACAGCCAGCAAAAGACCCACACCAAACCGATATGCAATgcg                                                                                                           >  6:306751/7‑149 (MQ=255)                                               gcgAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCcaca                                                                                                                               <  7:128294/121‑1 (MQ=255)                                               gcgAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCcaca                                                                                                                               >  8:128294/1‑121 (MQ=255)                                                cgAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCCCACCAAACCGCCGTGCACTGCGCAGAATCGc                                                                                                  >  3:190216/1‑149 (MQ=255)                                                  aaTGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCa                                                                                                <  3:228374/149‑1 (MQ=255)                                                       ggCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAGGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtg                                                                                           >  6:290879/1‑149 (MQ=255)                                                       ggCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtg                                                                                           <  6:165013/149‑1 (MQ=255)                                                       ggCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtg                                                                                           >  7:253788/1‑149 (MQ=255)                                                       ggCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtg                                                                                           <  3:343040/149‑1 (MQ=255)                                                       ggCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCAGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtg                                                                                           >  7:213498/1‑149 (MQ=255)                                                         ccTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtgt                                                                                          >  4:186021/1‑148 (MQ=255)                                                         ccTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGACCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgtgt                                                                                          >  5:182051/1‑148 (MQ=255)                                                             aCAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATc                                                                                      >  8:270343/1‑148 (MQ=255)                                                              cAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCTCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAGGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCg                                                                                     <  5:290879/148‑1 (MQ=255)                                                               aGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGAt                                                                                   >  5:299498/1‑149 (MQ=255)                                                               aGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGAt                                                                                   <  1:46102/149‑1 (MQ=255)                                                                     aCCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAAc                                                                               >  4:289378/1‑147 (MQ=255)                                                                     aCCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAAc                                                                               <  3:289378/147‑1 (MQ=255)                                                                         gTAACTCCCGGTAATCCTCGACCGCCTCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCTGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTg                                                                         <  5:306751/149‑1 (MQ=255)                                                                           aaCTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAACCCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGACCAGGCTgag                                                                       >  5:300149/1‑147 (MQ=255)                                                                                ccGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTc                                                                    <  7:270343/147‑1 (MQ=255)                                                                                          tCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACCCCAAACCGCCTTGACCTGCGCAGAACCGCCAGCGTGGAAGCGTGGAAAAGGCTGCGTTCCGGTTCGAAGcc                                                        >  7:155362/1‑149 (MQ=255)                                                                                            gACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCg                                                       <  7:144018/148‑1 (MQ=255)                                                                                            gACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGa                                                      >  6:345339/1‑149 (MQ=255)                                                                                                   gcAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGc                                               <  4:142085/149‑1 (MQ=255)                                                                                                      aaCTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTg                                            <  8:253788/149‑1 (MQ=255)                                                                                                             cAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAgg                                        <  6:182051/146‑1 (MQ=255)                                                                                                             cAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAgg                                        <  6:299498/146‑1 (MQ=255)                                                                                                                  tCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATccc                                 >  3:337184/1‑148 (MQ=255)                                                                                                                    ccGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaa                               >  3:89203/1‑148 (MQ=255)                                                                                                                     cgcgGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATTAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaaaa                             <  8:213498/149‑1 (MQ=255)                                                                                                                     cgcgGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaaaa                             >  1:282379/1‑149 (MQ=255)                                                                                                                           ccTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCgg                       <  5:345339/149‑1 (MQ=255)                                                                                                                             tGCTTTGCTATCCGGCTGGCGCACCGTCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaaaa                             <  6:32978/141‑1 (MQ=255)                                                                                                                             tGCTTTGCTATCCGGCTGGCGCACCGTCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaaaa                             >  5:32978/1‑141 (MQ=255)                                                                                                                              gCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCATAACTGCGGGCt                    >  1:124189/1‑149 (MQ=255)                                                                                                                                    cTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCt                                             >  1:120893/1‑118 (MQ=255)                                                                                                                                    cTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCt                                             <  2:120893/118‑1 (MQ=255)                                                                                                                                     tATCCGACTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTc             >  3:297070/1‑149 (MQ=255)                                                                                                                                       tCCCGCTGGCGCACCGCCATGAATTGCCCCACTCCAAACCGCCGTGCACTGCTCAGAATCGCCATCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcg            <  8:155362/148‑1 (MQ=255)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCAt                                                                                       >  5:119978/3‑71 (MQ=255)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCAt                                                                                       <  6:119978/69‑1 (MQ=38)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACa                                                                              >  5:101753/3‑80 (MQ=255)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACa                                                                              <  6:101753/78‑1 (MQ=255)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgcg          >  7:255605/3‑148 (MQ=255)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgcg          >  5:28689/3‑148 (MQ=255)                                                                                                                                         ttgCTGGCGGTGCGCCAGCAAGGGCCCCACACCAACCCGCCGTGCACTGCGCAGAATCGCCAGCGGGTCATCGATGAACAGGCGGCGTGACGGTTCGAAGCCGATATCCGCCGGCAGGGCATCCCAAAACTGCGGGCTGTCCTGcgcg          >  4:210170/3‑148 (MQ=255)                                                                                                                                                cgcACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATa  <  4:89203/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATA  >  NC_002947/313820‑314111

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