Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,555,989	C→A	5.3%	E170* (GAG→TAG) ‡	zwfB ←	glucose 6‑phosphate 1‑dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,555,989	0	C	A	5.3%	302.8 / 5.1	94	E170* (GAG→TAG) ‡	zwfB	glucose 6‑phosphate 1‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (51/38);  new base A (3/2);  total (54/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

TGTAGTGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAATGGCAT  >  NC_002947/4555842‑4556135                                                                                                                                                     |                                                                                                                                                    tGTAGTGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCg                                                                                                                                                    >  1:174763/1‑149 (MQ=255) tGTAGTGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCg                                                                                                                                                    >  7:207743/1‑149 (MQ=255)  gTAGTGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCgg                                                                                                                                                   >  3:108540/1‑149 (MQ=255)    aGTGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGgcg                                                                                                                                                 >  1:342554/1‑149 (MQ=255)      tGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCt                                                                                                                                               >  3:398497/1‑149 (MQ=255)      tGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCt                                                                                                                                               >  3:15866/1‑149 (MQ=255)        aTTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTgg                                                                                                                                             >  5:274892/1‑149 (MQ=255)         ttGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCGAGGCGCtcaa                                                                                                                                            <  6:254011/149‑4 (MQ=255)           gTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTgg                                                                                                                                             >  5:126120/1‑146 (MQ=255)           gTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTgg                                                                                                                                             <  6:126120/146‑1 (MQ=255)               cAGCACGATTCGAACAAGCCGTTGGAAAAGCGGCTGCCCAGAATATTCTGAACCGTCTCTTTGCCCAGCTAATGGTCGATGCGGAAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTc                                                                                                                                      >  5:67979/1‑149 (MQ=255)                  aaCGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGt                                                                                                                                   >  8:91071/1‑149 (MQ=255)                  aaCGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGt                                                                                                                                   >  4:168214/1‑149 (MQ=255)                  aaCGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGt                                                                                                                                   >  4:254220/1‑149 (MQ=255)                  aaCGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACAGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGAGCTGGCCAGGTCGGt                                                                                                                                   >  4:386293/1‑149 (MQ=255)                   aCGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTg                                                                                                                                  <  1:392069/149‑1 (MQ=255)                       ttCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCg                                                                                                                               >  3:297820/1‑148 (MQ=255)                          gAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAAc                                                                                                                            <  2:107056/148‑1 (MQ=255)                            aCAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGaggcgc                                                                                                                         <  2:342554/149‑6 (MQ=255)                             cAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGtt                                                                                                                        >  1:95915/1‑149 (MQ=255)                              aaGCCGTTGGAAAAGCGGCTGACCAGAATACTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGttt                                                                                                                       >  3:47330/1‑149 (MQ=255)                                gCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGttttt                                                                                                                     >  7:104821/1‑149 (MQ=255)                                   gTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCaa                                                                                                                  >  1:275900/1‑149 (MQ=255)                                   gTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCGTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCaa                                                                                                                  >  6:358664/1‑149 (MQ=255)                                    ttGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAac                                                                                                                 <  3:254220/149‑1 (MQ=255)                                          aaGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAaccaccac                                                                                                           <  4:398497/149‑1 (MQ=255)                                            gCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCAccc                                                                                                         >  1:86610/1‑149 (MQ=255)                                             cGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCTGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCg                                                                                                        >  1:272261/1‑149 (MQ=255)                                              ggCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCgg                                                                                                       <  4:15866/149‑1 (MQ=255)                                              ggCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCgg                                                                                                       <  8:104821/149‑1 (MQ=255)                                               gCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCAATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑c                                                                                                     <  8:207743/149‑1 (MQ=255)                                                cTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑c                                                                                                     >  2:15505/1‑148 (MQ=255)                                                   aCCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGaaa                                                                                                 <  7:91071/149‑1 (MQ=255)                                                   aCCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGaaa                                                                                                 <  1:10131/149‑1 (MQ=255)                                                   aCCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGa                                                                                                   <  5:298673/147‑1 (MQ=255)                                                     cAGATTATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGaaa                                                                                                 <  6:371725/147‑1 (MQ=255)                                                     cAGAGTATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGaaa                                                                                                 >  5:371725/1‑147 (MQ=255)                                                       gAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAgccg                                                                                             >  5:399302/1‑149 (MQ=255)                                                            ttCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCgg                                                                                         <  2:86610/148‑1 (MQ=255)                                                             tCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCa                                                                                       >  4:250846/1‑149 (MQ=255)                                                              cTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCa                                                                                       <  2:275900/148‑1 (MQ=255)                                                                 aaCCGTCTCTTTGCCCAGTTAATGGTCGATGCCGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGa                                                                                     <  2:90696/147‑1 (MQ=255)                                                                     gTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTc                                                                               <  2:95915/149‑1 (MQ=255)                                                                        tctTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTc                                                                               >  5:315965/1‑146 (MQ=255)                                                                          tttGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCa                                                                             <  6:315965/146‑1 (MQ=255)                                                                          tttGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGc                                                                           >  5:38848/1‑148 (MQ=255)                                                                             gCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAg                                                                        <  1:15505/148‑1 (MQ=255)                                                                              cccAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGcc                                                                      >  2:63868/1‑149 (MQ=255)                                                                              cccAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGcc                                                                      >  4:123873/1‑149 (MQ=255)                                                                               ccAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCg                                                                     <  3:168214/149‑1 (MQ=255)                                                                                 aGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTCCAACCACCACCCGG‑CGAAAGCCGCCgg                                                                                         >  1:265122/1‑127 (MQ=255)                                                                                 aGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTCCAACCACCACCCGG‑CGAAAGCCGCCgg                                                                                         <  2:265122/127‑1 (MQ=255)                                                                                     aaTGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGc                                                               >  2:62517/1‑149 (MQ=255)                                                                                      aTGGTCGATGCGTTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCTGGCCTGc                                                               <  3:386293/148‑1 (MQ=255)                                                                                      aTGGGCGTTGCGTTAGATCTGGGGCTCGCCCATCACCTTCAGCAGGCAGG‑GGTTGAGCGCCGCGGCGCTGGCCAGGTCGGTGCCTAACGGTTTTTCAACCACCACCCGG‑CTAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGc                                                               <  6:274892/148‑1 (MQ=255)                                                                                        ggTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGcccc                                                            >  4:380592/1‑149 (MQ=255)                                                                                         gTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGt                                                                                                                                   <  2:324171/78‑1 (MQ=255)                                                                                         gTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGt                                                                                                                                   >  1:324171/1‑78 (MQ=255)                                                                                         gTCGAGGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGcccc                                                            <  3:372878/148‑1 (MQ=255)                                                                                                    aGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGcc                                                >  4:155451/1‑149 (MQ=255)                                                                                                       tCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCAcc                                             >  3:84495/1‑149 (MQ=255)                                                                                                        cTGCCGCTCGCCCATCACCTTCAGCAGGCAGGCCGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCAcc                                             >  1:255066/1‑149 (MQ=255)                                                                                                           ccGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCgg                                         >  7:199951/1‑149 (MQ=255)                                                                                                           ccGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCg                                          <  3:250846/148‑1 (MQ=255)                                                                                                            cGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCggg                                        >  7:145524/1‑149 (MQ=255)                                                                                                            cGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCggg                                        <  3:173035/149‑1 (MQ=255)                                                                                                            cGCTCGCACATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCCCCCGG‑CGAAAGCCGCCGGCAGACTCGGCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCGCggg                                        >  7:173872/1‑149 (MQ=255)                                                                                                             gCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCgggg                                       <  7:165809/149‑1 (MQ=255)                                                                                                                 gCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGgaag                                    >  7:160014/1‑148 (MQ=255)                                                                                                                 gCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGgaag                                    <  2:272261/148‑1 (MQ=255)                                                                                                                 gCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGgaag                                    <  2:174763/148‑1 (MQ=255)                                                                                                                       cacaTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGACGAACGGTTTTTCAACCACCACACGG‑CGAAAGCCGCCGGCAGACTCGTCCAGAAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGcgcg                             >  5:329820/1‑149 (MQ=255)                                                                                                                        aCCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGTGAAGAAGCgcgg                            <  1:62517/149‑1 (MQ=255)                                                                                                                        aCCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCgcgg                            <  8:160014/149‑1 (MQ=255)                                                                                                                         ccTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCgcggc                           >  5:327381/1‑149 (MQ=255)                                                                                                                          cTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCgcggcg                          >  3:256114/1‑149 (MQ=255)                                                                                                                          cTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCgcggc                           >  5:335090/1‑148 (MQ=255)                                                                                                                          cTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCgcggc                           >  8:149631/1‑148 (MQ=255)                                                                                                                            tCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGAc                        >  6:159674/1‑149 (MQ=255)                                                                                                                             cagcagGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGAc                        <  1:268839/148‑1 (MQ=255)                                                                                                                              agcagGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGt                      >  5:50492/1‑149 (MQ=255)                                                                                                                                 aggcagg‑cGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGc                   <  8:199951/149‑1 (MQ=255)                                                                                                                                   gcagg‑cGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGc                   >  7:170316/1‑147 (MQ=255)                                                                                                                                   gcagg‑cGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCa                 <  6:38848/149‑1 (MQ=255)                                                                                                                                     agg‑cGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGc                   >  5:175106/1‑145 (MQ=255)                                                                                                                                          ggccaGCGCCGAGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAg            >  7:221890/5‑148 (MQ=255)                                                                                                                                          ggccaGCGCCGAGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAg            >  7:337421/5‑148 (MQ=255)                                                                                                                                            ccaGCGCCGAGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGaa          >  5:366998/3‑148 (MQ=255)                                                                                                                                              aGCGCCTCGGCGCTGGCCAGGTCTGTGCCGTACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGGCGGCAGACTCGTCCAGCAAGCCTGCCTGCCGGAGTCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAAt       <  6:329820/149‑1 (MQ=255)                                                                                                                                              aGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAAt       <  6:399302/149‑1 (MQ=255)                                                                                                                                              aGCGCCGAGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAAt       <  8:221890/149‑1 (MQ=255)                                                                                                                                               gcgcCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAATg      <  4:297820/149‑1 (MQ=255)                                                                                                                                                 gcCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGGTCGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAATgg     >  1:100665/1‑149 (MQ=255)                                                                                                                                                   cTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAATGGCAt  >  7:164760/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGTAGTGATTGTTCCAGAACGATTCGAACAAGCCGTTGGAAAAGCGGCTGACCAGAATATTCTGAACCGTCTCTTTGCCCAGGTAATGGTCGATGCGGTAGATCTGCCGCTCGCCCATCACCTTCAGCAGGCAGG‑CGTTGAGCGCCTCGGCGCTGGCCAGGTCGGTGCCGAACGGTTTTTCAACCACCACCCGG‑CGAAAGCCGCCGGCAGACTCGTCCAGCAAGCCGGCCTGCCCCAGGCGCTGCGCCACCTCGGGGAAGAAGCGCGGCGACGTGGCCAGATAGAAAATGGCAT  >  NC_002947/4555842‑4556135

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