Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,628,644	A→T	C99S (TGT→AGT)	yfgF ←	cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,644	0	A	T	98.8%	358.5 / NA	86	C99S (TGT→AGT)	yfgF	cyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (52/33);  minor base C (1/0);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGG  >  NC_000913/2628421‑2628879                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGaatgata                                                                                                                                                                                                                                   >  2:1651864/1‑234 (MQ=255)     aGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACtgag                                                                                                                                                                                                                                          >  1:638057/1‑223 (MQ=255)                               aGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc                                                                                                                                                                                                >  1:1025603/1‑239 (MQ=255)                                  taaaaaGACATGCCCCCACCCAGTTTTTCTCTGCTCGCAAGTAACCCAAAAAATGCAGCAAACGGAAAAAGTATCAA‑AAACAGCGTTGCAGA‑AAACCAAATCTGCCAGAATATACGCTGGAAAATT‑AACCGGGTATCC‑CCATGTAAAACG‑TTTTTCCCCGGGGGAGAAAAGCCCCTT‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                               <  2:1525073/219‑1 (MQ=255)                                            tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                                                                                                                                                     >  1:1525073/1‑221 (MQ=255)                                                    aCCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCata                                                                                                                                                                           >  2:1655239/1‑239 (MQ=255)                                                         aTTTTCCCTGCCCGCCACTAATCC‑AAAAAAGCAACGAAACAGAAAAATAAAAGAAAACAG‑GTTGGAGGAAAACACATTCTCCCAGAATATCCGC‑GGAAAATTAAAACGGGAATCC‑CTATGAAAACCG‑TTGTTGC‑CGGGGGGAAAAAGGCCCTGTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                                                                                                                            <  2:742720/217‑1 (MQ=255)                                                                   cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                                  >  1:742720/1‑217 (MQ=255)                                                                   cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaaa                                                                                                                                                             >  2:1898066/1‑238 (MQ=255)                                                                                  aCAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                                                                                                                                               >  1:2125204/1‑237 (MQ=255)                                                                                         cagcaAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGaa                                                                                                                                                                                                                                        >  1:1696892/1‑141 (MQ=255)                                                                                         cagcaAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAt                                                                                                                                      >  2:2371561/1‑239 (MQ=255)                                                                                            caAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGaatga                                                                                                                                                                                                                                     >  1:387752/1‑141 (MQ=255)                                                                                               aCTGAAAAATTATCAA‑AACCCGGTTTGAAGA‑AAACACATTCTGCCAAAATTTCCGTTGGGAAATA‑AACCGGGTATCC‑CAATGT‑AAACGTTTGTTCCGCGGGGGAGTAAAGCCCTTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc                                                                                                                                <  2:533426/239‑1 (MQ=255)                                                                                                           tCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAc                                                                                                                    >  1:1379799/1‑239 (MQ=255)                                                                                                                      cGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatcat                                                                                                            >  1:627251/1‑237 (MQ=255)                                                                                                                      cGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatcat                                                                                                            >  1:1262367/1‑237 (MQ=255)                                                                                                                       gTTGCAGA‑AACCAAAATCTGCCAAAATATCCGCGGGGAAATT‑AAACGTTTACCC‑CCATGAGAAGGT‑TTGTTACGCGGGGAAGAAAAGGCCCTT‑AAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                               <  2:245990/137‑1 (MQ=255)                                                                                                                       gTTGAAGA‑AAACACAATCTGCAAGATTATACGCGGGAAAATA‑AACCGGTTATCA‑CATGGTAAAACG‑TTGTTACGCCGGGGAGAAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatca                                                                                                                <  1:1360224/232‑1 (MQ=255)                                                                                                                        gttgaaaaaaaCAAATTCTGCGAAAATACGCGCAGAGAAATA‑AAACGGGTATCA‑CCATTTAAAAGT‑TTGTTACCCGGGGGAAAAAAGTCCCTG‑AACCCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcat                                                                                                                  <  2:1379799/222‑1 (MQ=255)                                                                                                                         tGCAGG‑AAACACAATCGGCCAGAATATCCGCTGGAAAATA‑AACCGGGTCTCC‑CCTTGTAAATCT‑TTTTTACGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTg                                                                                                        <  2:2267474/238‑1 (MQ=255)                                                                                                                              g‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcatc                                                                                                  >  1:2267474/1‑239 (MQ=255)                                                                                                                                      aaTCTGCCAGAATATACGCTGGAAAATT‑AAGCGGTTACCC‑CCATGTGAACCG‑TTTTTACCCGGGGGAAAAAAGCCCCTG‑TAACCGCCCCAACTGAGaatgataat                                                                                                                                                                                                                                 <  2:2378825/105‑1 (MQ=255)                                                                                                                                           gCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTTAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCATCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGc                                                                                      >  1:2057365/1‑239 (MQ=255)                                                                                                                                                aaTATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATaaaa                                                                                                                                                                                                                       >  1:2378825/1‑105 (MQ=255)                                                                                                                                                 atagACGCGGGAAAATT‑AAACGGTTATCA‑CCATGTAAAATT‑TTGTTACGCGGGGGAGTAAAGCCCCTG‑TACCCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcat                                                                                                   <  2:2224654/220‑1 (MQ=255)                                                                                                                                                      cGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAagag                                                                           >  2:2170631/1‑239 (MQ=255)                                                                                                                                                      cGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAagag                                                                           >  2:2281144/1‑239 (MQ=255)                                                                                                                                                          gggAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaaa                                                                       >  2:1921193/1‑239 (MQ=255)                                                                                                                                                          gggAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaaa                                                                       >  1:35888/1‑239 (MQ=255)                                                                                                                                                           ggAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATaaa                                                                                                                                                                                                                        >  1:833353/1‑93 (MQ=255)                                                                                                                                                           ggAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCaga                                                                                         >  1:2224654/1‑220 (MQ=255)                                                                                                                                                                 t‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCa                                                                 >  2:1063913/1‑238 (MQ=255)                                                                                                                                                                   aaagcggTATCACCCATGTAAAATG‑TTGTTGCGCGGGGAATAAAAGCCCCTG‑AAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCAc                                                                <  2:847509/234‑1 (MQ=255)                                                                                                                                                                      cGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAg                                                                                                                                                                                                                   >  1:727549/1‑88 (MQ=255)                                                                                                                                                                          tATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaagaa                                                        >  2:555295/1‑239 (MQ=255)                                                                                                                                                                                   tgtgAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGaat                                                                                                                                                                                                                                       >  1:399145/1‑56 (MQ=255)                                                                                                                                                                                   tgtgAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAAAAGgcgc                                                >  2:1131437/1‑239 (MQ=255)                                                                                                                                                                                    gtgAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGaatgata                                                                                                                                                                                                                                   >  2:1568834/1‑59 (MQ=255)                                                                                                                                                                                        aaCG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAg                                           >  1:2249787/1‑239 (MQ=255)                                                                                                                                                                                        aaCG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAg                                           >  2:1843750/1‑239 (MQ=255)                                                                                                                                                                                            ttttttACTCCGGAGAGAAAAGCCCCTG‑AAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTc                                                                  <  2:1873157/209‑1 (MQ=255)                                                                                                                                                                                               gttACGCCGGGGCAAAAAGACCTGG‑AAGACGCCCCAACTGAGAATGATAATGATAATAAAATGAGCCCTTAATGATATCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAGGCCGccc                                                                                                                               <  2:2062425/149‑1 (MQ=255)                                                                                                                                                                                                ttACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAAGTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAAt                                    >  2:2018145/1‑239 (MQ=255)                                                                                                                                                                                                      cGGGGAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaagaa                                                        >  1:1873157/1‑213 (MQ=255)                                                                                                                                                                                                          gAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa                                                                                                                     >  1:2062425/1‑148 (MQ=255)                                                                                                                                                                                                          gAGTAAAGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGaaa                          >  2:2390416/1‑239 (MQ=255)                                                                                                                                                                                                                aGACCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGc                    <  1:1898066/239‑1 (MQ=255)                                                                                                                                                                                                                  aCCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatca                                                                                                             >  2:208704/1‑148 (MQ=255)                                                                                                                                                                                                                  aCCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatca                                                                                                             <  1:208704/148‑1 (MQ=255)                                                                                                                                                                                                                  aCCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatca                                                                                                                <  2:975864/145‑1 (MQ=255)                                                                                                                                                                                                                  aCCCTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatca                                                                                                                >  1:975864/1‑145 (MQ=255)                                                                                                                                                                                                                   cccTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCagag                                                                                        <  2:1583238/168‑1 (MQ=255)                                                                                                                                                                                                                   cccTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCagag                                                                                        >  1:1583238/1‑168 (MQ=255)                                                                                                                                                                                                                   cccTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg                     <  1:1063913/235‑1 (MQ=255)                                                                                                                                                                                                                   cccTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg                     <  2:2125204/235‑1 (MQ=255)                                                                                                                                                                                                                   cccTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg                     <  1:2371561/235‑1 (MQ=255)                                                                                                                                                                                                                     cTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGc                                                                                             >  2:2208421/1‑161 (MQ=255)                                                                                                                                                                                                                     cTG‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGAt               >  1:143554/1‑239 (MQ=255)                                                                                                                                                                                                                      ctttAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGc                                                                                             <  1:2208421/158‑1 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                                  <  2:2046510/74‑1 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                                  >  1:2046510/1‑74 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                                                                                                                                               <  2:2334426/93‑1 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                                                                                                                                               >  1:2334426/1‑93 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatca                                                                                                             <  2:207791/143‑1 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatca                                                                                                             >  1:207791/1‑143 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg                     <  1:829803/231‑1 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg                     >  2:829803/1‑231 (MQ=255)                                                                                                                                                                                                                       g‑TAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTa             >  2:329395/1‑239 (MQ=255)                                                                                                                                                                                                                               ccccAATTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAaagaag       >  2:337682/1‑238 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGcc                                                                                                                                                                                                                 <  2:552229/36‑1 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGcc                                                                                                                                                                                                                 >  1:552229/1‑36 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                               >  1:135496/1‑38 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                               <  2:135496/38‑1 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                                  <  2:1332629/67‑1 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                                  >  1:1332629/1‑67 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCg                                                              >  2:2152175/1‑183 (MQ=255)                                                                                                                                                                                                                               ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCg                                                              <  1:2152175/183‑1 (MQ=255)                                                                                                                                                                                                                               ccccAACCGAGAATGTTAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACCAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGAATTAATGTGGGCGAACATGGGCAAGATTAaagaag       >  2:234265/1‑238 (MQ=255)                                                                                                                                                                                                                                cccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGc     <  1:2390416/239‑1 (MQ=255)                                                                                                                                                                                                                                 ccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                                                                                                                                    <  2:2122700/63‑1 (MQ=255)                                                                                                                                                                                                                                 ccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                                                                                                                                    >  1:2122700/1‑63 (MQ=255)                                                                                                                                                                                                                                   aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                                                                                                                                    <  2:1095473/61‑1 (MQ=255)                                                                                                                                                                                                                                   aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                                                                                                                                    >  1:1095473/1‑61 (MQ=255)                                                                                                                                                                                                                                   aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGGAGAAATATGCGCGa                                              >  2:1439488/1‑195 (MQ=255)                                                                                                                                                                                                                                   aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGGAGAAATATGCGCGa                                              <  1:1439488/195‑1 (MQ=255)                                                                                                                                                                                                                                   aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAgg  <  2:1025603/239‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            AGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAG‑AAACAGCGTTGCAGG‑AAACACAATCTGCCAGAATATACGCTGGGAAATT‑AAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAAAGAAGGCAGG  >  NC_000913/2628421‑2628879

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