Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,579,657	A→C	18.6%	V62G (GTG→GGG)	symE ←	toxic peptide regulated by antisense sRNA symR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,579,657	0	A	C	18.6%	76.9 / 2.1	53	V62G (GTG→GGG)	symE	toxic peptide regulated by antisense sRNA symR
Reads supporting (aligned to +/- strand):  ref base A (36/4);  new base C (3/8);  total (39/14)
Fisher's exact test for biased strand distribution p-value = 9.85e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGC  >  NC_000913/4579434‑4579882                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATAAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TAACtt                                                                                                                                                                                                                                  >  2:169532/1‑227 (MQ=255)      gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCTCCGGCGGCGGGCAGTTGGGCGGTTAGGGCAA‑TTCA‑G‑CCTTCTA‑TCACTTTGAcc                                                                                                                                                                                                                             >  2:438892/1‑226 (MQ=255)      gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACa                                                                                                                                                                                                                             >  2:87521/1‑227 (MQ=255)      gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTCCACACCTGGCGTAGCGAC‑TGCATTCTTTCGCTCTCCTCGGCGGCCGGGTGTGGGGCGGGGAGGGCCA‑TACA‑G‑CCTTCCA‑TCCCCTTTACAt                                                                                                                                                                                                                            >  1:183090/1‑228 (MQ=255)       ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGCCTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCCCCTGCTTTTGTTTACG‑CGCCGACGGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTg                                                                                                                                                                                                                                <  2:183090/223‑1 (MQ=255)       ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACAt                                                                                                                                                                                                                            >  2:29575/1‑227 (MQ=255)                 tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTg                                                                                                                                                                                                                 >  2:81484/1‑228 (MQ=255)                    aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGCCGGGGGGTTGGGCGTTGAGGACCA‑TACA‑G‑CCTTCCA‑TCACTTTGCCATCTACCGCCGGGCCa                                                                                                                                                                                                              >  2:142814/1‑228 (MQ=255)                         gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGCCATCTACCGCCGTGCCAGTgg                                                                                                                                                                                                          >  2:646924/1‑227 (MQ=255)                          tCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGAGGGGGGGACAA‑TTCA‑G‑CCTTCCC‑TCCCTTTGACATCTATCGCCGTGCCacgggc                                                                                                                                                                                                         >  2:163831/1‑222 (MQ=255)                             gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGGCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTTGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGGCAA‑TAAC‑G‑CCTTCCA‑TCACTTTGCCTGTTACCGCCGGGCCCTGGGCaaaa                                                                                                                                                                                                     >  1:529162/1‑228 (MQ=255)                                    gCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGGGCCAGTGGCAAAACcggc                                                                                                                                                                                                >  2:547671/1‑226 (MQ=255)                                      tCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAACTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TAAA‑G‑CCTTCCA‑TCACTTTGACATGTACCGCCGTGCCAGTGGCAAAACcggcggat                                                                                                                                                                                            >  2:231520/1‑226 (MQ=255)                                        ccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCtt                                                                                                                                                                                           >  2:214700/1‑227 (MQ=255)                                              gtTACCCGGGAGTAACCGCTGTTGCGCGGGTTGGTGTTACCCGGAATTCCCTCCAATAACCGCCGGCACCGGTTTTGGTTCCCG‑CGCCTGCGGTTGGCCCGCCGGGCG‑GGCGGCGTGGAGCAGTCCGCCCTCCCGGGGGGCGGGGGGTGGGGGGGTGAGGACAA‑TCCA‑G‑CCGTCCC‑CCCTTTGGACATCGACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGcc                                                                                                                                                                                     <  1:40931/227‑1 (MQ=255)                                                      ggAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGCGGCCAC‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCt                                                                                                                                                                             >  2:18552/1‑227 (MQ=255)                                                      gcaGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTACACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGCCGGGTGGTTGGGCGGTGAGGACCA‑TACA‑G‑CCTTCCA‑TCACATTGCCACTTTCGGCCGGGCCAGTGGGCAAACCGGACGCTTCCAGCCACttgccct                                                                                                                                                                            >  1:106020/3‑222 (MQ=255)                                                        aGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCT‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGGGGCTTCCAGCCACTGGCCTTTc                                                                                                                                                                          >  2:628120/1‑228 (MQ=255)                                                              cgcgtttcctggtttttGTGTTTCCCGGCAACCCCACCAAAAAAC‑CCCGGACCCGTTTTGTTTTTCGCCGCCGACGGTTGGTCCACTGGGCGTGGCGAC‑TGAATCAGTCCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAAAA‑TACA‑G‑CCTTCCA‑TCCCTTTGACATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTg                                                                                                                                                                    <  1:547671/214‑1 (MQ=255)                                                                  gTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGAGGTTGGGCGGTGAGGACAA‑TACC‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTAACCTTTCAGGGTGATgg                                                                                                                                                                >  1:174148/1‑228 (MQ=255)                                                                        ccttttttttttcacggacAACCATTCCAAAAACCCCTTGCCCTTGTTTTTTTTTACG‑CGCCGACGGGTTGCACCCCGGGCGGAGCGCA‑CGCATATGTTCGCCTCCCTCGTCGGCGGGTGGTTGGGGGGTGAGGAAAA‑TACA‑G‑CCTGCCA‑TACCTTGGCAATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGaa                                                                                                                                                          <  1:664177/210‑1 (MQ=255)                                                                                     aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACGGACCTTTCAGGGTGATGGCGGGGATT‑CGCCTGTGAt                                                                                                                                             >  2:352987/1‑227 (MQ=255)                                                                                            aTCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGGCCTTTCCGGGCGATGGCGGGGATA‑CGGCTGTGATCCggcgac                                                                                                                                     >  2:631680/1‑224 (MQ=255)                                                                                                           gCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCCCTTTGACAACTACCGCCGTGCCAGAGGCAAAACCGCCGGCTTCCAACCACTGACCTTTCCGGGGGAGGGCGGGAATA‑CCGCCGGCATCCCGGGAGCAACCCGCCCAACCg                                                                                                                      >  1:533736/1‑228 (MQ=255)                                                                                                           gCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGCCCTTTCAGGGGGGTGGCGGGCAAC‑AGGCCGGAATCCGGGtgcgccctcccaaaaacc                                                                                                                      >  1:104402/1‑211 (MQ=255)                                                                                                           gCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGGTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTAATGGCGGGCAGA‑CGACTGGAATACCGGGAGCGACTCCCATAACCg                                                                                                                      >  2:666858/1‑228 (MQ=255)                                                                                                              tGCACCTGGTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGGCAA‑TACA‑G‑CCTTCCA‑TCACTTTGCCCTCTACCGCCGTGCCAGTGGCAAAACCTGGGGCTTCCAGCCACTGACCTTCTCGGGTGATGGAGGGAATA‑AGTCTGTAATCCGGGGAACGACACCCCTAAcccccg                                                                                                                   >  1:385365/1‑224 (MQ=255)                                                                                                              tGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGTGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAAGA‑CGGCCTGAATCCggggaacgcatcggcgacacgacg                                                                                                                   >  1:663479/1‑207 (MQ=255)                                                                                                               gCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TCCA‑G‑CCTTCCA‑TCACTCTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCACCCACTGACCTTTCAGGGTGAGGGGGGGGATA‑CGGCGCGAATCCGGGTGGCCGCCCCCCTAACCGCCgg                                                                                                                  >  1:560220/1‑228 (MQ=255)                                                                                                                cACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGCGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCTGGTAGCGGCCAGCATAACCGACGGt                                                                                                                 >  2:526613/1‑228 (MQ=255)                                                                                                                     gCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCCACGGTGACAt                                                                                                            >  2:222775/1‑228 (MQ=255)                                                                                                                     gCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGAAAA‑CGGCTGTAATCCCGGGAGCGACTCCCATAAACGACGGGGaaaa                                                                                                            >  1:283452/1‑225 (MQ=255)                                                                                                                     gCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGGTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCAc                                                                                                                                                                                                                                    >  1:404057/1‑109 (MQ=255)                                                                                                                                                     cTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGAGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGCGGCAAAACCGGCGGCTTCCAGGCACTTCACCTTCAGGGTGACGTCGGGAATACCGGCGGCATACAGGGTA‑CGAATCGCACAAACGACGGGGGGATGTCGGTTATTTGCCgggggagcacccgcgcc                                                                             >  1:202031/1‑215 (MQ=255)                                                                                                                                                      tGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGAGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑AGGCTGTAATCCGGGTAGCGACTCGCCTAACCGACGGTGACACGACGGTGAGTTGCTGGGGGGACTTCTGCTCCg                                                                            >  2:196568/1‑228 (MQ=255)                                                                                                                                                             gtcctcatcATCCGTTCCCTCTCCCTGGGGGGGGGTGGGTGGG‑GGGGTGGACAA‑AACACGCCCGCCCA‑TCCTTTTGACTGTCGCCGCCGTGCCAGTGGCAAAGCCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACgg                                                                      <  1:646924/220‑1 (MQ=255)                                                                                                                                                             gCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAAACGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGtt                                                                    >  2:593889/1‑228 (MQ=255)                                                                                                                                                             gCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGAATA‑CGCCTGTAATCCGGGTAGCCACCCCCCTAACCGCAGGGGCCATAACGGTTATTTGCGgggggggacttggcttccaaagggt                                                                     >  1:667929/1‑207 (MQ=255)                                                                                                                                                               gAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAAGCCCGGGAGCGACTCGCATAAACGAAGGTGGCATGACGGCTATTTGCCGGGGGGACTTCTGCTTCG‑GACGGgtgt                                                                  >  2:562597/1‑228 (MQ=255)                                                                                                                                                                       tCAGTTGGCGCCCCCCGGGGGGGGGTGGTTGGGGGGGGGGGATAACAACA‑G‑CCTCCCT‑ACCCTTGGCCTTCCCCCGCGGGGCCAGGGGAAAGACCGGGGGCTTCCAGCCACTGACTTTTCAGGGTGAGGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAAt                                                             <  1:631680/227‑1 (MQ=255)                                                                                                                                                                        ctgttcccctcccggggaggggggTGGTGGGGCGGGAAGGACAA‑TACC‑G‑CCTTCCA‑TCATTTTGACTTCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGa                                                          <  2:45376/207‑1 (MQ=255)                                                                                                                                                                               ctctCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCCATAGAATGCg                                                     >  2:138534/1‑226 (MQ=255)                                                                                                                                                                               ctccctccggcggcggGTGGTTGGGGGTGGAGGACAA‑TACA‑G‑CCTTCCA‑TCCCTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCgtg                                                   <  2:427659/221‑1 (MQ=255)                                                                                                                                                                                    ctcGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGCGATGGCGGGAATA‑CGGCTTGAATCCCGGTAGCGACTCCCCTAAACGGCGGGGGCATCGCGGGTTTTTTCCGGGGGGGATTTCGCTTCGCCACGGTGGTGC‑ATAAAATGCgtgctagg                                              >  1:413227/1‑223 (MQ=255)                                                                                                                                                                                        ccggcggGTGGTGGGGGAGGGAGGAAAG‑TCCG‑T‑CCTTCCA‑TCATTTTGACATATACCGTCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCTCTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCAt                                           <  2:529162/226‑1 (MQ=255)                                                                                                                                                                                              ggtgg‑tGGGGGGGGCGGAAAA‑TAAC‑G‑CCTTCCA‑TCTCTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTa                                   <  2:172552/228‑1 (MQ=255)                                                                                                                                                                                              gggggtGGGGGGGTGAGGAAAA‑TACC‑C‑CCTTCCA‑TCCCTTTGCCATTTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCt                                    <  2:283452/225‑1 (MQ=255)                                                                                                                                                                                                  gtTGGGCGTGGAGGACAA‑CACC‑G‑‑CTTCCACTCATTTTGACTTTCACCGGCGTGCCAGGGGAAAAACCGGCGGTTTCCAGCCACTGCCCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATtc                                <  1:81484/228‑1 (MQ=255)                                                                                                                                                                                                     gggCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGGGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGGTGTGCAATAGAATGCGGGTCAGTCATAACTGCTCTTCTcc                             >  2:606760/1‑228 (MQ=255)                                                                                                                                                                                                       gcgtggAGGACAA‑TACG‑T‑CCCTCCA‑TACCTTTGCAATCGACCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCa                            <  2:426783/222‑1 (MQ=255)                                                                                                                                                                                                                acaa‑caCA‑G‑CCTTCCA‑TCTATTTGACATCTCCCGGCTGGCAAGTGTACAAACCGGGGGCATCCAGCCACTGACCTTTCAGAGTGATGGCGGGAATA‑CGGCTGTAAGCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGTGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtga                  <  2:560220/228‑1 (MQ=255)                                                                                                                                                                                                                            ccTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGGCGTAATCCGGGTAGCGACTCCCATAACCGACGGTGACATGGCGGGTCTTTGCCGGGGGGACTTCTGGTTCA‑GCCGGTTGTGCAATAGAATGCGCGTCCGTGCTAACTTCATTTCCCCCGGAATAgtgatttggataa         >  1:130458/1‑221 (MQ=255)                                                                                                                                                                                                                                    tCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGGGACATGACGGTTATTTCGCGGGGGGACCTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGGCAGTCATAACTGCTTTTCTCCAGGAAAAGGGATTGTGATTAgggcgtg  >  1:526140/1‑222 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    ACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACG‑CGCCGACAGTTTGCACACCTGGCGTAGCGAC‑TGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAA‑TACA‑G‑CCTTCCA‑TCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATA‑CGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCG‑AACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGC  >  NC_000913/4579434‑4579882

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