Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,817,197	G→A	R103C (CGC→TGC)	yqaB ←	fructose‑1‑phosphate phosphatase YqaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,817,197	0	G	A	96.2%	152.3 / ‑1.6	52	R103C (CGC→TGC)	yqaB	fructose‑1‑phosphate phosphatase YqaB
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (23/27);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

CACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCGACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTAC  >  NC_000913/2816955‑2817433                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              cACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCAt                                                                                                                                                                                                                                       >  1:356677/1‑251 (MQ=255)  ccGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCGACGACCATg                                                                                                                                                                                                                                      <  1:358392/250‑1 (MQ=255)  aCGGCGTCCATGCCTGCTGCACGGGCCGCCTGCATACCGAAATCGGCGACTTAAAAGACCACACACTGCGTAGGATGCAC‑CACACTACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTA‑ATCGGTGTTTGAC‑GTGATCGACAGCGACGACA‑GCGTCAAAATAATGTGGTAAACCCAGGGGCACCAACAATACCTCACAGATGTCGCTTTCAATCCACGTTCCT‑ACAGACTTTGGGccacgcccgtg                                                                                                                                                                                                                                      >  2:62132/1‑245 (MQ=255)      cGTCCCTGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTTCGCGAACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGCTTGAC‑GTGCTCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGCTGGCGCTTTCACTCCCCGCTCCT‑ACAGCCATACGGCAACGACCATGCCa                                                                                                                                                                                                                                   >  1:558644/1‑250 (MQ=255)            tGCCTGCTGCACGGGCCGCCTGAAAACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCAAAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATAGCGTAATCCCCGGTGCGCCAGCAATGCCTCCGCGATAGCGCTTTCACTCCCCGTTCCTTAC‑GAAATTGGGCAACGACCATGCCGACttttt                                                                                                                                                                                                                            >  2:123373/1‑250 (MQ=255)                    gcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATc                                                                                                                                                                                                                   >  1:72301/1‑251 (MQ=255)                    gcACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATc                                                                                                                                                                                                                   <  2:469516/251‑1 (MQ=255)                                tGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGACCGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGAGTTTAAC‑GTGATCGGCAGCGACG‑CGTGCGTCAAACTAATGGCCTAATCCCAGGTGCGCCAGCAATGACTCAGCAATAGAGATTTCACTCCCCGTTTCT‑CCCGCAATTGGGCAACGACCCAGCCAACTTTTAAACAC‑ATCCACAAGCGGaa                                                                                                                                                                                                        >  2:343478/1‑250 (MQ=255)                                 gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAg                                                                                                                                                                                                       <  1:559068/250‑1 (MQ=255)                                     cccGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGGTGCACGCCCA‑TACGCTGCGCGCACAACAAAAAGGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCAGTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGtt                                                                                                                                                                                                  <  2:505016/250‑1 (MQ=255)                                        gAAATAGGCATCGTCAAAGACCACACACTGAGTCGGTTGAACGCCCA‑TACGCTGCGCGCACCTCACCACAGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGAAGACG‑GCGGCAAAATAAAGGCGTAATCCCAGGTGCGCAAGCAATGCCTCCGCGCTGGCGCTTTCACTACCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCACACTTTTC‑CCACCATCCACAAGAGGAACCGGTTct                                                                                                                                                                                                >  2:361280/1‑249 (MQ=255)                                              ggCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCAAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTAt                                                                                                                                                                                         >  1:105553/1‑251 (MQ=255)                                                catgttCAAAGACCACACACTACGTCGGTTGCCCGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTAAC‑GTGATCGGCAGCGACGATG‑GCGTCAAAATAATTGCGGAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCCTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTTACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCCATcc                                                                                                                                                                                       >  2:142710/5‑251 (MQ=255)                                                       aaaGACCACACACTGGGTCGGTGGCACGCCCA‑TACGCGGCGCGAACACCAAAAATGGGTCGGGCGCGG‑TTTAA‑GGTGGTTGACCGGAATCGGAAGGGACGACG‑GCGCAAAAAAAAGGGCGTAACCCAAGGTGGGCCGGCAATGCCTAAGCGATG‑CGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑AGCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCAt                                                                                                                                                                               <  2:539279/251‑1 (MQ=255)                                                         aGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTCTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCCTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCGACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATa                                                                                                                                                                              <  1:6801/251‑1 (MQ=255)                                                           accatactcggcgtgggtaggaCGCCAA‑TATGCGGCGCGAACAACAAGAATGGGTCTGGCGCGTGTTTAG‑TGTGTTGGAC‑GAGATCGGCGACGACGACG‑TCGTCAAAATAAGTGCGTAATCCCAGGTGAGCCAGCAGTGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTGTTCACAAC‑ATCAACAAGAGGAAGCGGTTCGACGCTAGCCAGCAGCATAct                                                                                                                                                                            <  2:52608/230‑1 (MQ=255)                                                             cacctacTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAAGGGGTCTGGCGCTGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActtc                                                                                                                                                                          <  2:135899/246‑1 (MQ=255)                                                                 cacTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCCACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATcacgaaccctctttcttct                                                                                                                                                                      >  2:545908/1‑233 (MQ=255)                                                                                 cGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGaa                                                                                                                                                                                                        >  2:503419/1‑201 (MQ=255)                                                                                 cGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGaa                                                                                                                                                                                                        <  1:503419/201‑1 (MQ=255)                                                                                   cccA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcg                                                                                                                                                    <  1:424515/251‑1 (MQ=255)                                                                                   cccA‑TACGCTGCGCGCACAACAAAAATGGGTCTGGCGCGGGTTTAT‑GGGGTTCTAC‑GTGTTCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTTCGCCAGCAATGCCTCAGCGATGGCTCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACGTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCGAGCAGCATACTTCTTACTGTGTCTTTTGTTTCAcg                                                                                                                                                    <  2:72301/251‑1 (MQ=255)                                                                                      a‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACgcgc                                                                                                                                                 >  1:303779/1‑251 (MQ=255)                                                                                      a‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACgcgc                                                                                                                                                 >  1:370910/1‑251 (MQ=255)                                                                                           gCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAAc                                                                                                                                             <  1:164077/251‑1 (MQ=255)                                                                                                         aaaaaTGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCgag                                                                                                                               <  1:524843/251‑1 (MQ=255)                                                                                                                     ggCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGAtt                                                                                                                   <  1:108527/251‑1 (MQ=255)                                                                                                                                         tGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTg                                                                                                 <  1:96322/250‑1 (MQ=255)                                                                                                                                          gAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTg                                                                                                 >  1:235931/1‑249 (MQ=255)                                                                                                                                                     gaaGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTTCGCCAGCAATTCCCCAGCGATGGCGCTTTCCCTCCCCGTTCCT‑ACATCCATTGGGCAACGACCATGCCAACTTGTCTCCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATTCTTCTTACTGCTTCTGTGTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCa                                                                                     <  2:205899/249‑1 (MQ=255)                                                                                                                                                      cAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGcc                                                                                      >  1:67139/1‑249 (MQ=255)                                                                                                                                                        gCGACGACG‑GCGTCAAAATAATAGCGTAATCCCAGGTGCGCCAGCAATGCCCCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGATGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGt                                                                                  <  2:565891/251‑1 (MQ=255)                                                                                                                                                         cgacgaCG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTg                                                                                 <  1:109677/251‑1 (MQ=255)                                                                                                                                                          gacgcag‑gCGTCAAAATTATGGCGTAATCCCGGGGGCGCCAGCACTGCCTCAGCGATGGCGCTTTCACTCCCCGGTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTGTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAAGAATTGCCTGAGCAATTCGCCAGGTgg                                                                                <  2:493739/245‑1 (MQ=255)                                                                                                                                                                g‑gCGTCAAAATAATGGCGTAATCCCATGTGCTCCAGCACGGCCTCATCGTGGGCGCTTTCACTCCCCGTGCCA‑ACAGCCATTGGGCAACGTCCATGCCAACTGTTCTCTAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACTCCAGGTGGGCgttc                                                                          <  2:333066/251‑4 (MQ=255)                                                                                                                                                                g‑gCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGcc                                                                                                   >  1:451677/1‑226 (MQ=255)                                                                                                                                                                g‑gCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGcc                                                                                                   <  2:451677/226‑1 (MQ=255)                                                                                                                                                                       aaaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTa                                                                     >  1:554696/1‑250 (MQ=255)                                                                                                                                                                          ataatGGCGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACCCCCCGTCCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAgcg                                                                 <  2:67139/251‑1 (MQ=255)                                                                                                                                                                                 cGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAatcat                                                          <  1:250718/251‑1 (MQ=255)                                                                                                                                                                                       cccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGGGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCCCCTg                                                    >  2:428003/1‑251 (MQ=255)                                                                                                                                                                                               cgcCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTg                                                                                                                      >  1:400163/1‑177 (MQ=255)                                                                                                                                                                                                gcCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGt                                            >  1:557265/1‑250 (MQ=255)                                                                                                                                                                                                       aTGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGa                                       >  2:109484/1‑248 (MQ=255)                                                                                                                                                                                                           cTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGt                                >  1:193730/1‑251 (MQ=255)                                                                                                                                                                                                             cAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATTCTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATGCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGt                                <  1:361280/249‑1 (MQ=255)                                                                                                                                                                                                             cAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAg                              <  1:7556/251‑1 (MQ=255)                                                                                                                                                                                                                          ttcctcccccGTTCCT‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATAc                   <  2:487264/243‑1 (MQ=255)                                                                                                                                                                                                                                   cGTTCCT‑ACGGCCATTGGGCAACGACCCTGCCTACCTTTCACCAC‑ACCAACAAGAGGAAGCGTTGGGACGCTATCCAGCAGCATATTTCTTACTGCTTCTGTTTGTTCCCGCGTTAACGCATTCGGGGCGAGATCGGCCTGATTCAGCTTAATAATTTCCTGAGCCATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGACTATCGTACTGAAGACCGGAGTGCCCTAATGCATCGCGCCACg        <  2:14035/251‑1 (MQ=255)                                                                                                                                                                                                                                         t‑ACAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTAc  >  1:144095/1‑251 (MQ=255)                                                                                                                                                                                                                                           aCAGCCATTGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTa   <  1:375749/249‑1 (MQ=255)                                                                                                                                                                                                                                                   tGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGCGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGAc                                   >  2:510184/1‑209 (MQ=255)                                                                                                                                                                                                                                                   tGGGCAACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGAc                                   <  1:510184/209‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              CACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCA‑TACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTAT‑GGTGTTTGAC‑GTGATCGGCAGCGACGACG‑GCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCT‑ACAGCCATTGGGCGACGACCATGCCAACTTTTCACCAC‑ATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTAC  >  NC_000913/2816955‑2817433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑