Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,401	T→A	7.8%	N121Y (AAC→TAC) ‡	dctB ←	C4‑dicarboxylate transport sensor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,598,401	0	T	A	7.8%	179.5 / 2.8	60	N121Y (AAC→TAC) ‡	dctB	C4‑dicarboxylate transport sensor protein
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base A (3/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACCGCGCTGGCGAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGG  >  NC_002947/1598256‑1598545                                                                                                                                                       |                                                                                                                                                 ttCACCGCGCTGGCGAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCt                                                                                                                                               >  1:160340/1‑149 (MQ=255)    aCCGCGCTGGCGAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTggc                                                                                                                                            >  3:327865/1‑149 (MQ=255)             gCGAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGc                                                                                                                                   <  7:241746/149‑1 (MQ=255)             gCGAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGc                                                                                                                                   >  4:37850/1‑149 (MQ=255)               gAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGCCATTGCCGCCAGCAACTggcggctgcc                                                                                                                                     >  6:30497/1‑143 (MQ=255)               gAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGGTGGGGGGGAAACGGAGGTAGTTGGAGCCCACGTAGATGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATgg                                                                                                                                    >  1:112547/1‑149 (MQ=255)                          ccGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCgg                                                                                                                      >  4:277708/1‑149 (MQ=255)                           cGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGGCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGt                                                                                                                     <  7:270729/149‑1 (MQ=255)                           cGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGATGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGt                                                                                                                     <  3:68589/149‑1 (MQ=255)                                        gTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGa                                                                                                         <  4:327865/148‑1 (MQ=255)                                         tCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAgg                                                                                                       >  8:170637/1‑149 (MQ=255)                                               ccACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTggag                                                                                                 >  1:300327/1‑149 (MQ=255)                                                           aGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTCGCTGCTGGGCA‑‑GCCGCCTGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGGCGAGCAGCTCGAGGGTGGAGGAATTggcggc                                                                                      <  4:4983/148‑1 (MQ=255)                                                               gCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCAt                                                                                 >  7:11592/1‑149 (MQ=255)                                                                ccACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCAtt                                                                                <  2:160340/149‑1 (MQ=255)                                                                  aCTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTc                                                                                                                    >  1:82115/1‑111 (MQ=255)                                                                  aCTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTc                                                                                                                    <  2:82115/111‑1 (MQ=255)                                                                   cTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGAt                                                                             <  8:22638/149‑1 (MQ=255)                                                                   cTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGAt                                                                             >  6:347221/1‑149 (MQ=255)                                                                         tGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTc                                                                       >  3:104592/1‑149 (MQ=255)                                                                           gCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTc                                                                       <  5:347221/147‑1 (MQ=255)                                                                           gCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTc                                                                       <  3:277708/147‑1 (MQ=255)                                                                            cTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAg                                                                    <  7:170637/149‑1 (MQ=255)                                                                             tGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATg                                                                            <  2:121924/140‑1 (MQ=255)                                                                             tGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATg                                                                            >  1:121924/1‑140 (MQ=255)                                                                               cGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCt                                                                  <  2:300327/148‑1 (MQ=255)                                                                                  gTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAg                                                              <  2:112547/149‑1 (MQ=255)                                                                                  gTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAg                                                              <  7:316332/149‑1 (MQ=255)                                                                                       cTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTggc                                                                                                                                            <  1:10711/66‑1 (MQ=255)                                                                                       cTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTggc                                                                                                                                            >  2:10711/1‑66 (MQ=255)                                                                                       cTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTgg                                                                                                   >  3:56381/1‑107 (MQ=255)                                                                                       cTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTgg                                                                                                   <  4:56381/107‑1 (MQ=255)                                                                                        ttGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAg                                                        >  6:287526/1‑149 (MQ=255)                                                                                          gAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAggg                                                      >  8:247492/1‑149 (MQ=255)                                                                                           aaGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGc                                                     >  2:46178/1‑149 (MQ=255)                                                                                                   gCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTggc                                                                                         >  4:175182/1‑105 (MQ=255)                                                                                                   gCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGaa                                             <  7:247492/149‑1 (MQ=255)                                                                                                         aC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCtctc                                       <  5:287526/149‑4 (MQ=255)                                                                                                          c‑‑‑cGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTtc                                      <  1:46178/149‑1 (MQ=255)                                                                                                          c‑‑‑cGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGAACGCCCTGAACCtcaag                                      <  8:11592/149‑5 (MQ=255)                                                                                                                   ttGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGa                                 <  4:104592/148‑1 (MQ=255)                                                                                                                           gACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGt                                                                                                                          <  8:62522/51‑1 (MQ=255)                                                                                                                               tatctgctgGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCaa                    <  5:174405/146‑1 (MQ=255)                                                                                                                                  ctgctgGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTg                                               >  1:91406/1‑119 (MQ=255)                                                                                                                                  ctgctgGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTg                                               <  2:91406/119‑1 (MQ=255)                                                                                                                                    gctgGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAAc                                            >  4:239058/1‑120 (MQ=255)                                                                                                                                    gctgGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAAc                                            <  3:239058/120‑1 (MQ=255)                                                                                                                                       gggCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGa            >  1:304695/1‑149 (MQ=255)                                                                                                                                        ggCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCaa                    <  3:140792/140‑1 (MQ=255)                                                                                                                                        ggCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCaa                    >  4:140792/1‑140 (MQ=255)                                                                                                                                           attgccgccAGCAACTGGCGGCTATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCttt                                        >  1:210107/4‑119 (MQ=255)                                                                                                                                           attgccgccAGCAACTGGCGGCTATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCttt                                        <  2:210107/116‑1 (MQ=255)                                                                                                                                            ttgccgccAGCAACTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGAt           >  2:134457/3‑147 (MQ=255)                                                                                                                                               ccgccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCt      >  1:51040/1‑149 (MQ=255)                                                                                                                                               ccgccAGCAACTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGc       <  1:134457/148‑1 (MQ=255)                                                                                                                                                 gccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGc                                                     <  1:182542/100‑1 (MQ=255)                                                                                                                                                 gccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGc                                                     >  2:182542/1‑100 (MQ=255)                                                                                                                                                  ccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAAc                                            >  3:135561/1‑108 (MQ=255)                                                                                                                                                  ccAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAAc                                            <  4:135561/108‑1 (MQ=255)                                                                                                                                                   cAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCgg  <  2:51040/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 TTCACCGCGCTGGCGAGGAAGTAGCCCGGCACGCCACTGGTCACGCCCACGGCATAGAAGCGGCCACTGCCCTGGCTGCGGGTTTGCTTGAAGTAGGGGCGGAAAC‑‑‑CGTAGTTGGAGCCGACGTAGCTGCTGGGCA‑‑GCCGCCAGTTGCTGGCGGCAATGGCCAGCCCGGTGCGGTCGAGCAGCTCGAGGGTGGAGGAATTGGCGGCGCCATTGATGCGTTCGAGCTTGAGGTTCAGGGCGTTCTGAACCTTTTCGTCGACCGGGCCGCGCAAGGCTGCGATCAGCTCCGG  >  NC_002947/1598256‑1598545

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