Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,098,171	C→G	7.5%	M407I (ATG→ATC)	acsA‑I ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,098,171	0	C	G	7.5%	292.4 / 11.8	94	M407I (ATG→ATC)	acsA‑I	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base G (3/4);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

ACCAGGTGTCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGA  >  NC_002947/5098023‑5098316                                                                                                                                                      |                                                                                                                                                  accaGGTGTCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCcatc                                                                                                                                                   <  1:92519/149‑1 (MQ=255)      gtgtCGACGATCGGGCAACGCTCCATGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTCGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGgc                                                                                                                                              <  3:245680/149‑1 (MQ=255)       tgtCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGgcg                                                                                                                                             >  5:208636/1‑149 (MQ=255)         tCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCACGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGgcgc                                                                                                                                            >  6:403758/1‑148 (MQ=255)           gacgaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCATTCGCGCCATGATGgccg                                                                                                                                            >  6:266388/1‑145 (MQ=255)            acgaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAAt                                                                                                                                        >  2:9998/1‑149 (MQ=255)              gaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCTTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                      >  4:418314/1‑149 (MQ=255)              gaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                      <  7:400941/149‑1 (MQ=255)              gaTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                      <  6:279356/149‑1 (MQ=255)               aTCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                      >  2:315972/1‑148 (MQ=255)                tCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCa                                                                                                                                    <  1:9998/149‑1 (MQ=255)                tCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCATTCGCGCCATGATGGC‑CGAAGGGCa                                                                                                                                    <  3:119081/149‑1 (MQ=255)                  gggCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACTCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCATTCGCGCCATGATGGC‑CGAAGGGCAgg                                                                                                                                  <  8:271932/149‑2 (MQ=255)                  gggCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCATTCGCGCCATGATGGC‑CGAAGGGCAgg                                                                                                                                  <  5:266388/149‑2 (MQ=255)                           tCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCtgt                                                                                                                         >  7:284342/1‑149 (MQ=255)                              ttGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCtgtg                                                                                                                        >  5:345494/1‑147 (MQ=255)                                gCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTa                                                                                                                      <  5:55746/147‑1 (MQ=255)                                  cGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTaga                                                                                                                    <  7:162684/147‑1 (MQ=255)                                  cGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTaga                                                                                                                    >  8:162684/1‑147 (MQ=255)                                   gACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGa                                                                                                                 >  4:203843/1‑149 (MQ=255)                                    aCGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGAt                                                                                                                >  5:298832/1‑149 (MQ=255)                                     cGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATg                                                                                                               <  5:37148/149‑1 (MQ=255)                                            gtagtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGAcc                                                                                                         <  3:395826/148‑1 (MQ=255)                                            gtagtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCt                                                                                                        >  4:102936/1‑149 (MQ=255)                                              agtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTtg                                                                                                      >  7:371414/1‑149 (MQ=255)                                               gtaCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTtgt                                                                                                     <  6:208636/149‑1 (MQ=255)                                                taCCAGTTCCAGGCCTCAGGGTTTATCGGCTCACCCACGGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCACCAGCTGGCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTtgt                                                                                                     <  5:359110/148‑1 (MQ=255)                                                taCCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTtgtg                                                                                                    <  5:403758/149‑1 (MQ=255)                                                   cAGTTCCAGGCGTCTCGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCtt                                                                                                 <  4:321488/149‑1 (MQ=255)                                                    aGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACAGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTg                                                                                                >  5:149796/1‑149 (MQ=255)                                                     gTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                      >  5:56306/1‑110 (MQ=255)                                                     gTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATgg                                                                                                                                      <  6:56306/110‑1 (MQ=255)                                                        ccAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTtgtg                                                                                                    <  2:86603/141‑1 (MQ=255)                                                        ccAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTtgtg                                                                                                    >  1:86603/1‑141 (MQ=255)                                                         cAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTcgac                                                                                           >  8:149468/1‑149 (MQ=255)                                                           ggCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCtt                                                                                                 <  7:388325/141‑1 (MQ=255)                                                           ggCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCtt                                                                                                 >  8:388325/1‑141 (MQ=255)                                                            gCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGAt                                                                                        >  2:153465/1‑149 (MQ=255)                                                             ccTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATc                                                                                       <  5:6923/149‑1 (MQ=255)                                                              cTCAGGGTAGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTAG‑GCCATCATGGCGCGAATCGCAGTTGGCGCTGAGTAGAGGATGTTGACCTAGTGCTTGTCGACGATc                                                                                       >  6:173116/1‑148 (MQ=255)                                                                     ttGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCtt                                                                                                       <  7:124110/125‑1 (MQ=255)                                                                     ttGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCtt                                                                                                       >  8:124110/1‑125 (MQ=255)                                                                        aTCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCtt                                                                                                       >  7:158625/1‑122 (MQ=255)                                                                        aTCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCtt                                                                                                       <  8:158625/122‑1 (MQ=255)                                                                        aTCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATAcgc                                                                             <  3:135616/148‑1 (MQ=255)                                                                           ggCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGt                                                                         <  7:286346/149‑1 (MQ=255)                                                                            gCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCGTGTCGACGATCTTCGATACGCGGGt                                                                         >  6:123608/1‑148 (MQ=255)                                                                                acccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGTGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAAt                                                                      <  1:315972/147‑1 (MQ=255)                                                                                acccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCTGGTAATg                                                                     <  5:173116/148‑1 (MQ=255)                                                                                  ccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTg                                                                                                                                >  7:228224/1‑87 (MQ=255)                                                                                  ccaccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTg                                                                                                                                <  8:228224/87‑1 (MQ=255)                                                                                    accGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCgg                                                                <  2:45249/149‑1 (MQ=255)                                                                                     ccGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCggg                                                               <  2:9597/149‑1 (MQ=255)                                                                                            cagcagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTgg                                                        <  6:298832/149‑1 (MQ=255)                                                                                            cagcagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTgg                                                        <  6:149796/149‑1 (MQ=255)                                                                                             agcagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTggg                                                       >  5:331041/1‑149 (MQ=255)                                                                                               cagACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCa                                                     <  3:203843/149‑1 (MQ=255)                                                                                                agACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCAc                                                    <  3:102936/149‑1 (MQ=255)                                                                                                  aCGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCggg                                                               <  3:72951/136‑1 (MQ=255)                                                                                                  aCGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCggg                                                               >  4:72951/1‑136 (MQ=255)                                                                                                  aCGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGc                                                  >  7:61862/1‑149 (MQ=255)                                                                                                  aCGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGc                                                  >  7:64725/1‑149 (MQ=255)                                                                                                    gCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTaa                                                                       >  7:311825/1‑126 (MQ=255)                                                                                                    gCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTaa                                                                       <  8:311825/126‑1 (MQ=255)                                                                                                    gCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCt                                                >  3:421178/1‑149 (MQ=255)                                                                                                            gAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAAc                                          <  7:149468/147‑1 (MQ=255)                                                                                                            gAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAg                                        >  4:104839/1‑149 (MQ=255)                                                                                                             aaCCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCgg                                                                           >  4:321484/1‑113 (MQ=255)                                                                                                             aaCCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAg                                        >  2:303433/1‑148 (MQ=255)                                                                                                              aCCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGc                                       >  4:203581/1‑148 (MQ=255)                                                                                                                     gcgcCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAggtggt                               <  8:64725/149‑1 (MQ=255)                                                                                                                      cgcCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGgtggtg                              <  3:203581/149‑1 (MQ=255)                                                                                                                       gcCTTCAACATCTTCCTGCCCTTCG‑GCCATCATGGGGCGAATTGCAGTTGGCTCTGTGTAGATGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGgtggtg                              <  3:116074/148‑1 (MQ=255)                                                                                                                       gcCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGtggtgg                             <  6:331041/149‑1 (MQ=255)                                                                                                                         cTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGgcg                           >  3:129366/1‑149 (MQ=255)                                                                                                                           tCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAggt                                  >  6:350200/1‑140 (MQ=255)                                                                                                                           tCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAggt                                  <  5:350200/140‑1 (MQ=255)                                                                                                                           tCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGcc                         >  8:361052/1‑149 (MQ=255)                                                                                                                              aCAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGtt                      <  3:104839/149‑1 (MQ=255)                                                                                                                              aCAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGtt                      >  5:54015/1‑149 (MQ=255)                                                                                                                              aCAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGtt                      <  4:421178/149‑1 (MQ=255)                                                                                                                                aGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTgg                    >  7:36511/1‑149 (MQ=255)                                                                                                                                  ctgcctgcCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGc                   >  6:334912/1‑148 (MQ=255)                                                                                                                                     cctgccCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGc               >  3:214162/1‑149 (MQ=255)                                                                                                                                         cccTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGc           >  7:11885/1‑149 (MQ=255)                                                                                                                                          ccTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGcc          >  1:387232/1‑149 (MQ=255)                                                                                                                                          cattcgcgCCATGATGGC‑CGAAGGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGc                   <  6:225411/133‑1 (MQ=38)                                                                                                                                          cattcgcgCCATGATGGC‑CGAAGGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGc                   >  5:225411/8‑140 (MQ=38)                                                                                                                                          cattcgcgCCATGATGGC‑CGAAGGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGcc          >  5:240845/8‑149 (MQ=38)                                                                                                                                            ttCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGTCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCg         >  6:328033/1‑148 (MQ=255)                                                                                                                                            ttCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGcc                         <  2:133994/132‑1 (MQ=255)                                                                                                                                            ttCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGcc                         >  1:133994/1‑132 (MQ=255)                                                                                                                                              cG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTa       <  5:328033/148‑1 (MQ=255)                                                                                                                                                  ccATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAg                                    >  3:136412/1‑116 (MQ=255)                                                                                                                                                  ccATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAg                                    <  4:136412/116‑1 (MQ=255)                                                                                                                                                   catcatGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGa  >  2:76483/1‑149 (MQ=255)                                                                                                                                                   catcatGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGa  <  6:54015/149‑1 (MQ=255)                                                                                                                                                   catcatGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGa  >  2:363657/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ACCAGGTGTCGACGATCGGGCAACGCTCCTTGCCGACGGTCTTGTAGTACCAGTTCCAGGCCTCAGGGTTGATCGGCTCACCCACCGACCCCAGCAGACGCAGGCTGGAACCGTCGGCGCCTTCAACAGCTGCCTGCCCTTCG‑GCCATCATGGCGCGAATGGCAGTTGGCGCTGTGTAGAGGATGTTGACCTTGTGCTTGTCGACGATCTTCGATACGCGGGTAATGTCCGGGTAGTTGGGCACGCCTTCGAACAGCAAGGTGGTGGCGCCGTTGGCCAGCGGGCCGTAGACGA  >  NC_002947/5098023‑5098316

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