Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,151,689	G→T	32.8%	intergenic (+43/+22)	yegJ → / ← yegK	uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,689	0	G	T	32.8%	17.8 / 5.6	36	intergenic (+43/+22)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base T (7/5);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

GAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTT‑ACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCT  >  NC_000913/2151491‑2151927                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aattatggcaggaaggaaaaaaACAAAAAATAAGTTATTTATGCATGGTTTTATCCGCGGGGGAAAAAAATGT‑CCGGCGTTAATACTGTTTATCTTTTTTTGGGCCCCTAACAAAAAGAAGAACCCGACGAACTAAGAGAAAAATTAGTTCGCTAA‑CCACCCTTTATTTTTATTTTTTAAC‑CGGGGGG‑TTGAT‑TTGTTTTC‑TC‑A‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGcc                                                                                                                                                                                                          <  2:1011943/221‑1 (MQ=255)  aaaggggaagggAGGGAAAAAACAAAAAAAAAGTAATTGAACACGGGTTTTATACGCG‑GGGGAAAAAATTTAACGCCGGTTTTCCATTTGACCTTTTATTGTCCCCCAATCAAAAAAAAAAAGCGGACGAGTAAAAAGAAAATTTAGTTCGCTAA‑TCAGCCGTTATATTTTTTTTTAAAA‑CTGGCGG‑TTAAA‑TTGTTGTT‑TC‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                                                                                                                                                                                                         <  1:1037307/229‑1 (MQ=255)              agtaaaaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATCCGCGGGGGAAAAAAATTT‑GCGGCGTTTATACCATTGATCCTTTATTGGTCCCCTACCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTATTTCGTTAA‑TCCGCCGTTATTTTTTTTGTTTAAA‑CGGGGGG‑TTAAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                              <  1:1060838/233‑1 (MQ=255)              aGTGAAAAACAAAGAAAGAAGATATTTCAGACTGGATGTATCCGCGGGGGGATAAAATTT‑ACGGCGTTTATCCCTTTGATCTTTTATTGGTCCCCTACCCAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTTTTGTTTAAA‑CTGGGGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                              <  1:1001456/238‑1 (MQ=255)                    aaaaaaaaaaaaaaGATATTGCAGACTGGTTTTTATCGCGGGGTGAAAAAATTG‑ACGCCGTTTATCCCATTTCTCTTTTATTGGTCCCCTACCAAAAAGAAAAAGCCGACGAGCTAAGAGCAAATTTATTTCTCTAA‑CCAGCCTTTATATTTATTTTTTAAC‑CGGGCGT‑TGGAA‑TTTTTTTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                                                                                                                       <  2:446524/231‑1 (MQ=255)                                                                        tt‑ACGGCGTTTATACCATTGATCTTTTATGGGTCCCCTACCAAAAAAAAGAGCCCAACGAGTAAAAAGAAAAATTAGTCCGTTAG‑TCACCAGTTATATTTATTTTTTAAC‑CGGCGGG‑TTAAT‑TTGTTTTC‑TT‑A‑AACCGCCATGAGAAACTCACTCTGCCCAAAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                                   <  1:260894/239‑1 (MQ=255)                                                                                    acaCCTTTGATCTTTTATTGGTCACCTACACAAAAAAAGAAGCCGAGGAGAAAAGAGAAAATTTTTTTCGATAA‑TCAGCTGTTATTTTTATTGTTTAAA‑GGGGCGT‑TTGAA‑TTGTTGTC‑TA‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                                                                                                                        <  1:1032405/237‑1 (MQ=255)                                                                                           tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttg                                                                                                                 >  1:233142/1‑239 (MQ=255)                                                                                           tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑ATGTTTTG‑TT‑C‑GAACCCCACGAGAAACTCACCATTCCCACCACGCCAGCGCCCACCTTTTATCATCATCCGTACACTCCTTTCCCGCCGGGCGGGACAAAAACAGCTTCAACAAATAAGGAGATAAATCGatgtg                                                                                                                 >  2:805371/1‑235 (MQ=255)                                                                                               ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt                                                                                                               >  1:116902/1‑237 (MQ=255)                                                                                                   tATGGGCCACCTACCAAAAAGAAGAACCCGACGAGCTAAGAGAAAAATTAGTTCGCTAA‑CCAGCCTTTATATTTATTGTTTAAA‑CGGGCGT‑TTGAT‑TTGTTTTC‑TC‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCg                                                                                                         <  1:559861/239‑1 (MQ=255)                                                                                                       ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTTTTGGC‑TT‑T‑CACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGCTTTTTCATTCTACGTACGCTCGCTTTCCGCCCGGCTTGACAAAAACTGCTTTAACAATTCAGGGAATAAATCCAGGTGTTTCTTCggttc                                                                                                     >  2:1020073/1‑235 (MQ=255)                                                                                                        gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TTGTTGTC‑TT‑C‑TACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAAAGTTTTATCATTATCCGTACGCTCGTTTACCGACGGGCTGGGCACAAAACGCTTTCACAATTCAGGTAATAAATCGAGTTGTTCCTTCGtt                                                                                                       >  2:765178/1‑236 (MQ=255)                                                                                                        gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGGA‑TTGTTTTT‑TT‑C‑AACCGCCATGAGGAACTCACTCTGCCCACAGCGACAGCGCCAACGTTTTTTCATCATCCGTACGCCCGTTTACCGCCGGGCTGGGCCACAACTGCTTTTACAATTTCGGGAATAAAACGAGTTGGTCCttgggtgct                                                                                                    >  2:446165/1‑231 (MQ=255)                                                                                                        gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGGA‑ATTTTTTT‑TT‑T‑AACCGCCATGAGAAACTCACTCTTCCCCCAGCGCCCGTGCCCACCTTTTATCATTATCCGGAAGCTCCTTTACCACCCGGCCGGGCAAAAACTTCTTTTACAATTCAGGGAATAAATCGCGGTTGTCCTGCGGtgct                                                                                                    >  2:21152/1‑239 (MQ=255)                                                                                                         tCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTTTTGTT‑TT‑C‑AACCGCCCTGAGAAACTTACTCTGCCCACCACGCCAGTGCCAACCTTTTTTCATCATCCGTACGCCCGTTTTCCGCCCGGCCGGACAAAAACTGCTTTTACACATTAGGTAATAAATCCAGTTGTTTCtttgtttct                                                                                                    >  2:866416/1‑230 (MQ=255)                                                                                                             ttaaaaaaaaGAAAAAGCCGAGAAGACAAAAGAAAAATT‑GTTCCGTAACTCGGCCTTTAT‑TTTTTTTTTTAAA‑CCGGGGT‑GTGTA‑TTTTTTTC‑TTCC‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccgc                                                                                               <  2:116902/233‑1 (MQ=255)                                                                                                                  aaaaaaaAGAAGCCGCCGGGCAAAGAAAAAAATTAGTTCGCAAA‑ACACCCTTTAT‑TTTTTTTTTTAAA‑CCGGGGGTTTGAA‑TTTTTGTC‑TC‑A‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTc                                                                                          <  1:805371/238‑1 (MQ=255)                                                                                                                          gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTGA‑TTTTTGTC‑TT‑T‑AACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACCCCCCTTTACCCCCGGGCCGGGCAAAAACTCGTTTAACCATTCCGGGAATTAATCGAGTTGGTCCTTCGTTTCTGCCCCCCGTCCA‑TTGgaa                                                                                 >  2:860063/1‑239 (MQ=255)                                                                                                                          gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTTTC‑TT‑T‑AAACGCCATTAGAAACTCACTCTGCCCACAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCAGGCTGGGCAAAAACTGCTTTAACAATTCAGGTAATAAATAGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGaaa                                                                                 >  1:746787/1‑239 (MQ=255)                                                                                                                           aagCGGACGAGTAAAGAGAAAATTTATTTCGTTAA‑CCAGCTGTTAT‑TTTTTTTTTTAAA‑AGGGGGG‑TTTAATTTGTTGTC‑TC‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGaaaa                                                                                <  1:217792/239‑1 (MQ=255)                                                                                                                              cgaaggaGCAAAGAGAAAATTTGGTCCCCTAA‑TCAGCGGTTATATTTTTTGTTAAAA‑C‑GGGGG‑GTGAA‑TTTTTGTCTTT‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGaaaaa                                                                               <  1:765178/232‑1 (MQ=255)                                                                                                                                              aaaaTAAGCTCGCAAA‑CCAGCGTTTATATGTATTTTTTAAA‑CGGGGGG‑TTGAT‑TTGTGGTC‑TT‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGcac                                                             <  1:1126355/239‑1 (MQ=255)                                                                                                                                                   tcgtccgttAA‑TCATCTGTTATATTTCTTTTTAAAA‑CGGGGGG‑TTAAT‑TTGTT‑TC‑TT‑CTAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATgag                                                        <  1:426931/231‑1 (MQ=255)                                                                                                                                                    tattccgtaaaccacccGTTATTTTTATTGTTAAAA‑CGGGGGG‑TGAAA‑TTGTTGCC‑TC‑C‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATgag                                                        <  1:1020073/224‑1 (MQ=255)                                                                                                                                                     gTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑GTTTTTTC‑TT‑C‑AAACGCCCTTAGACAGTCAACCTGCCCACACCGGCCGTGGCCACGGTTTATCATCATCCCTAAGCCGCTGTACCCCCGGGCTGGACCAACACATCTTTTAAAATTCCAGTAATAAATCCAGTTTGTCCTATGTGTCTGGCGCCCGGCGC‑TTTAAAAAAAGGGTAAAGAAAGGCACAttcgga                                                      >  2:1142029/1‑234 (MQ=255)                                                                                                                                                     gTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTAACGCCGGGGCGGACAAAAACTGCTTTAACAATTTCGGTAATAAATCGCGTTGTTTCTGCGTTGCTTCCGCCAGTCCC‑TTGAAAAACGGGGGAAAAAAAAGCACATgagga                                                      >  1:218977/1‑239 (MQ=255)                                                                                                                                                                    cGTTATA‑TTATTTTTTAAACCGGGCGG‑TTGAT‑TTGTTGTC‑TA‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGc                                         <  1:453539/237‑1 (MQ=255)                                                                                                                                                                      ttATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TTTTTTTC‑TT‑C‑AAACGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCCACGTTTTTTTCTTATCCGTACACTCGTTTACCGCCGGGCCGGGCAAAAACTGCTTTAACAATTCAGGGAATAAATAGCGTTGTTCCTGCGTTGCTGCCGCCAGACCA‑TTGGAAAACGGGGTAAAGAAAGGCACATGAGGGGGA‑TTATCCAGCctat                                     >  2:191291/1‑239 (MQ=255)                                                                                                                                                                             tattGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTTCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAAGCACATGAGGAGAA‑TTATCCAGCATATTTAACg                               >  1:1051334/1‑238 (MQ=255)                                                                                                                                                                               tgtttaaaaaggggggg‑TGGAT‑TTGTTGTT‑TA‑C‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGCATATTTAACGcc                             <  2:218977/224‑1 (MQ=255)                                                                                                                                                                                tGTTTAAA‑CTGGCGG‑TTTAA‑TTTTTTTT‑TT‑C‑TACCGCCTTGAAAAACTCACTCTGCCCCCAGCGCCAGGGCCCACGTTTTATAATCATCCGTACGCCCCTTTCCCGCCGGGCCGGGCAAAAACTGCTTTAACAATTCAGGTAATAAAACGAGTTGGTCCTGCGTTGCTGCCGCCACACCC‑TTGAAAAACGGGGGAAAGAAAAGCACATGAGGCGTA‑TTTTccctttttcttattccaa                            >  2:1147820/1‑222 (MQ=255)                                                                                                                                                                                tGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCCCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTTCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGGGGAGAA‑TTATCCAGCATATTTAACGCCAg                           >  1:152446/1‑239 (MQ=255)                                                                                                                                                                                      taaacctgggggtaaatttgtttTC‑TT‑A‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGCATATTTAACGCCAGTCg                        <  1:197892/222‑1 (MQ=255)                                                                                                                                                                                          tGGCGG‑TTAAT‑TTTTTGTC‑TT‑C‑ACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGCATATTTAACGCCAGTCGCTGGa                   <  1:583138/238‑1 (MQ=255)                                                                                                                                                                                             cGG‑TTGAA‑TTGTTGTC‑TT‑T‑AACCGCCATGAGAAACTCACTCTGCCCACAGCCCCCGTGCCAACGTTTTATCCTCATCCATACGCTGGTTTTACGCCCGGCTGGACAAAAAACACTTTAACAAATCAGGTAATAAATTCAGGTGTTCCTTCGTTTCTTCCGCCCGTTCATTTCACAAACGGGGAAAACAAAGCCA‑ATGGGGAGAATTTATCCA‑CATAttacacgaccaccgtgggataac               >  2:366050/1‑218 (MQ=255)                                                                                                                                                                                             cGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTaaag                                                                     >  1:120927/1‑186 (MQ=255)                                                                                                                                                                                                            gtC‑TT‑C‑ATCCGCCATGAGAAACTCCCTCTGCCCACAGCTCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAtgct  >  2:460266/1‑239 (MQ=255)                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    GAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTT‑ACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAA‑TCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TTGTTGTC‑TT‑C‑AACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCA‑TTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAA‑TTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCT  >  NC_000913/2151491‑2151927

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