Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,253,812	G→T	6.4%	G44V (GGC→GTC)	leuB →	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,253,812	0	G	T	6.4%	199.6 / 2.7	62	G44V (GGC→GTC)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base T (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GAATGGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCT  >  NC_002947/2253664‑2253959                                                                                                                                                      |                                                                                                                                                    gAATGGATAGAGGAAAGCATGAGCTAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTgg                                                                                                                                                     >  7:257637/1‑149 (MQ=255) gAATGGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCACCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTgg                                                                                                                                                     >  2:303901/1‑149 (MQ=255)  aaTGGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GATAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGc                                                                                                                                                    >  1:252989/1‑149 (MQ=255)  aaTGGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGc                                                                                                                                                    <  4:57394/149‑1 (MQ=255)     ggATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGt                                                                                                                                                 >  8:215645/1‑149 (MQ=255)       aTAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GATAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTgc                                                                                                                                               <  2:252989/149‑1 (MQ=255)       aTAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTgc                                                                                                                                               >  6:94379/1‑149 (MQ=255)        tAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTgcc                                                                                                                                              <  7:338822/149‑1 (MQ=255)        tAGAGGAAAGCATGAGCAAGCAGAATCTGATTCTCCCAGCTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATggc                                                                                                                                              >  6:132154/1‑147 (MQ=255)             gAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCAt                                                                                                                                         <  6:125823/149‑1 (MQ=255)               aaGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCg                                                                                                                                       <  1:59017/149‑1 (MQ=255)               aaGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCg                                                                                                                                       <  5:128415/149‑1 (MQ=255)                aGCATGAGCAAGCATATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAATGTACTGGAGCTGGCCAACAGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCg                                                                                                                                       <  1:108996/149‑1 (MQ=255)                    tGAGCAAGCAGATTCTGATTCTCCCAGGGGACGGTATCGGCCCGGAACTCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAg                                                                                                                                  <  4:190824/149‑1 (MQ=255)                           gcaGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGAACGgtg                                                                                                                           >  5:182765/1‑149 (MQ=255)                            caGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGgtgt                                                                                                                          <  5:290967/149‑1 (MQ=255)                            caGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATGGCGGCACCGCCAATCACgtcgt                                                                                                                          <  7:98848/149‑5 (MQ=37)                             agcttctgattctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCCGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGCCCAGGCCGGtgtg                                                                                                                         >  4:96356/4‑149 (MQ=255)                             aGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGtgtg                                                                                                                         <  5:94379/149‑1 (MQ=255)                              gattctgattctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGc                                                                                                                        <  8:257637/149‑1 (MQ=255)                               attctgattctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATGGCGGCACCGCCAATCACGTCGTGcg                                                                                                                       >  2:226196/1‑148 (MQ=25)                                     attctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGc                                                                                                                 <  3:96211/149‑1 (MQ=255)                                     attctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGc                                                                                                                 >  8:222907/1‑149 (MQ=255)                                     attctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATGGCGGCACCGCCAATCACGTCGTGcgcca                                                                                                                    <  5:132154/146‑5 (MQ=25)                                      ttctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGcc                                                                                                                <  6:180132/149‑1 (MQ=255)                                         acccAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTACCGCCATCGACAAGCACGGTGTGCCGCTGGCcga                                                                                                              <  1:303901/147‑1 (MQ=255)                                                  aCGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCt                                                                                                     <  6:92929/148‑1 (MQ=255)                                                  aCGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTg                                                                                                    <  7:215645/149‑1 (MQ=255)                                                  aCGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTg                                                                                                    >  7:331556/1‑149 (MQ=255)                                                        tCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGt                                                                                              >  3:196291/1‑149 (MQ=255)                                                          ggCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTgc                                                                                            >  2:355378/1‑149 (MQ=255)                                                               ggAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTgcgcgc                                                                                        >  3:11940/1‑148 (MQ=255)                                                                   aTCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGc                                                                                   <  3:96356/149‑1 (MQ=255)                                                                       tGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGAt                                                                               >  5:60280/1‑149 (MQ=255)                                                                          ccGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATg                                                                              >  6:183037/1‑147 (MQ=255)                                                                             aGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGtg                                                                         <  4:196291/149‑1 (MQ=255)                                                                                   tCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGctgctg                                                                   >  4:79662/1‑149 (MQ=255)                                                                                     aaGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTggg                                                                 >  6:148187/1‑149 (MQ=255)                                                                                        gTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGGGCTGCTGGGGGc                                                              >  8:179489/1‑149 (MQ=255)                                                                                               aGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTgggc                                                       <  3:79662/149‑1 (MQ=255)                                                                                                gCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTgggcg                                                      <  5:161345/149‑1 (MQ=255)                                                                                                 cTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTgggcgg                                                     <  5:148187/149‑1 (MQ=255)                                                                                                     ccAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGcc                                                  >  4:80563/1‑148 (MQ=255)                                                                                                       aaC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCg                                                 >  6:247482/1‑147 (MQ=255)                                                                                                       aaC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGaa                                               <  1:355378/149‑1 (MQ=255)                                                                                                         c‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGt                                             <  8:331556/149‑1 (MQ=255)                                                                                                            aCAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                           <  3:80563/149‑1 (MQ=255)                                                                                                            aCAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                           <  5:247482/149‑1 (MQ=255)                                                                                                                                 gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                           >  6:187678/1‑128 (MQ=255)                                                                                                                                 gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                           <  5:187678/128‑1 (MQ=255)                                                                                                                                 gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGCGACAtt                      >  5:177924/1‑149 (MQ=255)                                                                                                                                 gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACAtt                      >  4:376494/1‑149 (MQ=255)                                                                                                                                  ccTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgag                               <  8:329521/139‑1 (MQ=255)                                                                                                                                  ccTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgag                               >  7:329521/1‑139 (MQ=255)                                                                                                                                    tGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGc                   >  6:45506/1‑149 (MQ=255)                                                                                                                                      gcgcACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGcc                  <  3:376494/148‑1 (MQ=255)                                                                                                                                       cgcACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCg                >  5:127199/1‑149 (MQ=255)                                                                                                                                          acgacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGcc                                                                                                                <  7:128714/50‑1 (MQ=255)                                                                                                                                          acgacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGcc                                                                                                                >  8:128714/1‑50 (MQ=255)                                                                                                                                            gacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAgcg           >  6:151310/1‑149 (MQ=255)                                                                                                                                            gacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAgcg           <  6:60280/149‑1 (MQ=255)                                                                                                                                             acgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgagcg                             <  6:107857/130‑1 (MQ=255)                                                                                                                                             acgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgagcg                             >  5:107857/1‑130 (MQ=255)                                                                                                                                                   ttGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCctg    <  6:182765/149‑1 (MQ=255)                                                                                                                                                     ggcggcACCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCctgc   <  6:155065/148‑1 (MQ=255)                                                                                                                                                      gCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCctgct  >  6:84011/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GAATGGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAAC‑GACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCT  >  NC_002947/2253664‑2253959

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