Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,992,567	T→C	20.3%	E84G (GAG→GGG)	pgsA ←	phosphatidylglycerophosphate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,992,567	0	T	C	20.3%	116.6 / 2.2	56	E84G (GAG→GGG)	pgsA	phosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base T (28/11);  new base C (7/5);  total (37/19)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACCGCAGCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGG  >  NC_000913/1992331‑1992802                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 gTACCGCAGCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGACCATCGTTGCCGCCGGAAATGTCACCCACC‑AGCTTGGATAATGCTCg                                                                                                                                                                                                                                                >  2:847588/1‑239 (MQ=255)   aCCGCAGCCACAAAGAAAAGGGCAATCCCGGCGTACCTAACCCAAATGGTCGGCCGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCg                                                                                                                                                                                                                                                <  2:439172/237‑1 (MQ=255)            acaAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGAGATAATTCCCGGT‑T‑ACCAGcg                                                                                                                                                                                                                                    >  2:723192/1‑237 (MQ=255)            acaAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGACGATCATCGTGGCCACCGGTCATGTCACCCACA‑ACCTGAGATAAAACACGGT‑A‑ACCAacac                                                                                                                                                                                                                                   >  1:46979/1‑239 (MQ=255)                   aaaGTGCAATCCCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCg                                                                                                                                                                                                                            <  2:337154/239‑1 (MQ=255)                            tACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGGGATAATGCCCGGT‑T‑ACCAGCCCAATGGCAATAGCCACg                                                                                                                                                                                                                   >  2:17327/1‑239 (MQ=255)                            atccGGCGTACTCAACCAAAAGTTCCGGGACCCAAAGTCGGTATGCCAAGGCCACCATCGGGCCAGTGGTTTCACCTTTCACAATCCAGGAGAC‑GGCGACACTTCTGCGTTTACCCAACT‑CCTCCATCCATTCGCGTTGCGCAGAAATAATAATGT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGG‑T‑ACCAGCCCCATGGCGATAGCCACg                                                                                                                                                                                                                   <  1:562044/237‑1 (MQ=255)                                 gCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCTGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCCCC‑AGCTGCGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCAcaaaa                                                                                                                                                                                                                >  2:184867/1‑233 (MQ=255)                                 gCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCACCCATTCGCGTAGCGCAGAAATAATAATTT‑CCCGGGCCATCATCGTTGCCGCCGGTACTTCAACCCACC‑AGCTGTGATAATGCACGGT‑T‑CCAAGCACCACGGCGAAAACCCAGAGaaa                                                                                                                                                                                                              >  1:145279/1‑238 (MQ=255)                                 gCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGACACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAAAACTAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGGGATAATGCTCGGT‑T‑ACCAGCACCAGGGCGATAGCCACGaaaaa                                                                                                                                                                                                              >  1:233187/1‑238 (MQ=255)                                 gCGTACTCAACCCAAATGTTCGGACGACACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAGATAATAATTT‑CACGGGCAACCATCGTTGCAGCCGGTATTGTCACCCCCC‑AACTGTGCAAATGCCTGTG‑T‑ACCAGCCCCACGGGGGAAGCCAccaaa                                                                                                                                                                                                                >  1:655797/1‑233 (MQ=255)                                      cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATACTAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCCCGGT‑T‑ACCAGCACCATGCCGATACCCACCAGAACTTTa                                                                                                                                                                                                          >  1:132626/1‑238 (MQ=255)                                            ccAAAGTTTGGGCCGCCACAGCAGCAATGCCACACCAACCATGTGGGAAGGGGTTTTCCCTTTCCCAACCCA‑GTGACAGGCCACGCTACTGCGTTTACCCAACC‑CCGCCATCCACTCGCGTAGCGCAGAAATAATAATTG‑CGCGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCGCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAGCTTTATCTGCca                                                                                                                                                                                                   <  1:859127/239‑1 (MQ=255)                                              aaaTGTTCGGACGCCACAGCAGACATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGAGTTTACCCAACT‑CCGCCATCCATTTCCGTACCGCAGAAATAATAATTT‑CACGGGCGCTCACCGTGTCCGCCGGTAATGTAACCCACC‑ACCTGAGATAATGCCCGGT‑T‑CCACCCACCATAGCGAGAGCCACCAAAACTATATCTGCcaac                                                                                                                                                                                                 >  1:16841/1‑237 (MQ=255)                                                aTGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGTTAGCCACGAGAACTTTATCTGCCACAgg                                                                                                                                                                                               >  2:74811/1‑239 (MQ=255)                                                     cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCCCGAGAACTTTATCTGCCACAGGGTcca                                                                                                                                                                                          >  2:329920/1‑237 (MQ=255)                                                     cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCCCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCAAAGGGGCac                                                                                                                                                                                          >  2:130624/1‑238 (MQ=255)                                                        cccccACAGCAGCTAGGCCATAGCGACCATCTGGGCGGTGTTTTTCACTTCCCCAATCGAGGAGAC‑GGCCACGAGACTGCGTTTGCCCAACT‑CCGCCAGCCATGCCCGTGGTGCAGAAATAAGAATTC‑AAGGGGCGATCAGCGTTGCCGCCGGTAATGTCACCCACC‑AGGTGTTATAATGCCCGGT‑T‑ACCAGCACCATGGCGATAGCAGCGAGAACTTTATCTGCCACAGGGTCAAGGa                                                                                                                                                                                       <  1:330185/236‑1 (MQ=255)                                                                 cagcCATGCCAACGCCTCCTTCTGGGCAGGGTTTTTCACTTCCCCAATCCAGGGGAC‑CGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCaa                                                                                                                                                                              <  2:233187/239‑1 (MQ=255)                                                                             cGCCCCAATCTGGGCAGGGGTTTTCCCTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTActct                                                                                                                                                                 <  2:895384/239‑1 (MQ=255)                                                                                  ccaTCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCCCGAAAACTTTATCTGCCAAAGGGCCAGGGAAAGAACCAAACC‑GGGACCTCTGGtt                                                                                                                                                             >  1:895384/1‑238 (MQ=255)                                                                                  ccaTCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCCCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAAAACTTTATCTGCCACAGGGTCCAGGAAAGAACCAAACC‑GGGTACCCTGGTTc                                                                                                                                                            >  1:817148/1‑239 (MQ=255)                                                                                           gtgTGGTTTCCAGTCTCCCAATCCAGGGGAC‑GGCCACGCTCGTGCGTTGACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTG‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTGATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAgcggcg                                                                                                                                                    <  2:710065/236‑1 (MQ=255)                                                                                           ccgttgTTTCCACTTTCCCAAACCAGGAAAG‑GGCCCAGTTACGTGGTTTACCCAAGC‑CCGCCATCCATTCGGGTAGCGAAGAAATAATAATTT‑CAAGGGCGATCAGCGTTGCCGCCGGTAATGTCACCAACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAgcggcg                                                                                                                                                    <  1:373282/233‑1 (MQ=255)                                                                                                  ttccccTTCCCCAATCCAGGAGCC‑GGCCCCGCTTCTGCGTTTACCCAACC‑CCGCCATCCCTTCGCGTAGCGCAAAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGaa                                                                                                                                            <  1:48606/234‑1 (MQ=255)                                                                                                        tttCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCa                                                                                                                                       >  2:48606/1‑238 (MQ=255)                                                                                                          tCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAACTAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACTAGCAACATGGCGATAGCCACGAGACCTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAACGGCGTCACAGAAACACATCg                                                                                                                                    >  2:866968/1‑239 (MQ=255)                                                                                                                tCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTTcccacaacaccaaacacaac                                                                                                                              >  2:262905/1‑221 (MQ=255)                                                                                                                           ggCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CAcc                                                                                                                   >  2:651609/1‑238 (MQ=255)                                                                                                                           ggCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACAAGCACCATGGCGATAGCCACGAGACCTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGCACTCTGGTTCCAGCGGCGTGCCAGTAAACCATCGAACcagacagccccc                                                                                                                     >  2:237702/1‑233 (MQ=255)                                                                                                                                cGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGCAATGTCACCCCCC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCCCCATGGCGCTAGCCACGAGCACTTTATCTGCCACAGGGTCAAGGAAAGCACCAACAC‑GGGCACCCTGGCTCCAGCGGCGTGCCAGAAAACCATCGAACCAGCCAAC‑‑CAccgccggg                                                                                                             >  1:438566/1‑237 (MQ=255)                                                                                                                                  gcattccGTTTACCCAACTCCCCCCCACCATTC‑CGAGACGCCAAAAAAAAAATTTCCCCGGCCGAGCACCGTGCCCGCCGTACATGTCCCCCTCC‑GGCTGTGAAAATGCCCGGT‑T‑CCCAGCACCATGGCGATAGCCACGAGACCTTTATCTGCCACAGGGTCAAGAAAAGCACCAACCC‑GGGGACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCATTCAGT‑‑CACCGCCGCGa                                                                                                            <  1:469361/233‑1 (MQ=255)                                                                                                                                          tttACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTCCCAAAAAACAATCGAACCAGTCAGC‑‑CACC‑CCCCGACACCGaaaaa                                                                                                  >  1:71139/1‑238 (MQ=255)                                                                                                                                           ttACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CAACGCCGCGCCGCAAAAAAt                                                                                                  >  1:372292/1‑239 (MQ=255)                                                                                                                                             aCCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGACCTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCCGAAAACCATCGAACCCGTAAGT‑‑ACCCCCCCCGACGCAGAAAAtaa                                                                                                >  1:280897/1‑237 (MQ=255)                                                                                                                                             aCCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAAc                                                                                                                                 >  2:750156/1‑208 (MQ=255)                                                                                                                                                         ccatccatTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCCGAAAATGAGCGCggcggcc                                                                                     >  2:12481/1‑238 (MQ=255)                                                                                                                                                         ccatccatTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAACTGAGCGCGGCGGCa                                                                                     >  2:787242/1‑239 (MQ=255)                                                                                                                                                                  cctttggcccAGAAATAATATGTC‑CACGGGGTCTCCTCGGTGCCGCCGGTAAGGT‑ACCGACCGAGTTGTATAAACGCCGGGT‑T‑ACCGGACCCCTGGTGATGCCCCGGAAAGCTTTATCGGCCACAGGGTAAAGGAAAGCACCAAACC‑GGGTACCCTGGTTCCAGCGGCGTGCCAGAAAACCAGCGAACCAGCCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGAc                                                                            <  1:202719/230‑1 (MQ=255)                                                                                                                                                                      tAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGCCGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAgg                                                                        >  2:578440/1‑239 (MQ=255)                                                                                                                                                                                   ataatTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGCCCAGGGGCCAGGCAGATa                                                           >  2:431263/1‑239 (MQ=255)                                                                                                                                                                                   ataatTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGAACAGGTGACAGGCAGATa                                                           >  2:629804/1‑239 (MQ=255)                                                                                                                                                                                   ataatTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGCTAGT‑‑CACCGCCGCGACGCAGAAAATTAACCCGGGGGCAAACGGAGACCCGGGGACAGGCAGATa                                                           >  2:178386/1‑239 (MQ=255)                                                                                                                                                                                     aatTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGGGACAGGCAGATaaa                                                         >  2:523772/1‑239 (MQ=255)                                                                                                                                                                                                      ttcGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCAAAGGGTCAAGGAAAGCACCAACCC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGaa                                         <  2:280897/238‑1 (MQ=255)                                                                                                                                                                                                      tgctgcgcacccgGTAATTCCCCCCTCC‑ATCGGTGAAATTGC‑CGGTCT‑ACCAACCCCAGGGCGATACCCACGAGATCTTTACCGGCCCCGGGGTCAAGGAAAGC‑CCAAACCGGGGTACGCGGGTTCCAGCGGAGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGaa                                         <  1:787242/229‑1 (MQ=255)                                                                                                                                                                                                                ccgGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAgg                               >  2:586384/1‑239 (MQ=255)                                                                                                                                                                                                                 cggtatctccccccccc‑AGGTTTAATAATGCTCGGT‑T‑ACCCGCCCCATGGCGATCGCCAGGAAACTTTTACCGGCCACGGGGTCAAGAAAAGCACCAACCC‑GGGTATCCGGGTTCCAGCGGCGTGCCAGAAACCCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGa                              <  1:586384/228‑1 (MQ=255)                                                                                                                                                                                                                    tgaTGTCCCCCCCC‑AGCTGTGATAATGCTCGGT‑T‑CCCAGCCCGATGGCGATAGCCACGAGAACTTTATCGGCCACAGGGTCAAGAAAAGCCCCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGTCACACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGa                              <  2:597576/236‑1 (MQ=255)                                                                                                                                                                                                                          cccacaca‑aGCGGTGATCATGCGCGGT‑T‑ACCAGCACCATGGCGATAGCCGCGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGa                     <  2:145279/235‑1 (MQ=255)                                                                                                                                                                                                                           ccccccccatgtgtgATAATGCCCGTT‑TCACCA‑CACCATGGCGATAGCCCCGAAAACTTTAGCGGCCAAAGGGCCAAGGAAAGCCCCAAACC‑GGGTACTTTGGTTCCAGCGGCGCGCCAGAAACCCAGCGACCCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGa                     <  1:578440/228‑1 (MQ=255)                                                                                                                                                                                                                            ccccaccacggttgAAAACGCCCGGT‑T‑CCCAACACCTTGGCAATAGCCACGGAAACTTTATCTGCCACAGGGTCAAGGAAAGCACAAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGaa                    <  2:246660/227‑1 (MQ=255)                                                                                                                                                                                                                                 c‑aGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGCTaa                                                          >  2:666647/1‑195 (MQ=255)                                                                                                                                                                                                                                 c‑aGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGCGACAGGCAGATAAAAGACCAATACAAAGAATGGGTTAAGGATAACACGGAACAgtgg              >  2:533403/1‑238 (MQ=255)                                                                                                                                                                                                                                    gCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGGCAGGCAGATAAAAGACCAACAAAAAAAAAGGGGAAAAGATGacaagagaaaggtaaa            >  1:397586/1‑226 (MQ=255)                                                                                                                                                                                                                                              tGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGAGGACACGGAACAGTGTAAGccagtccct  >  2:158079/1‑231 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GTACCGCAGCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAC‑GGCCACGCTACTGCGTTTACCCAACT‑CCGCCATCCATTCGCGTAGCGCAGAAATAATAATTT‑CACGGGCGATCATCGTTGCCGCCGGTAATGTCACCCACC‑AGCTGTGATAATGCTCGGT‑T‑ACCAGCACCATGGCGATAGCCACGAGAACTTTATCTGCCACAGGGTCAAGGAAAGCACCAAACC‑GGGTACTCTGGTTCCAGCGGCGTGCCAGAAAACCATCGAACCAGTCAGT‑‑CACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGG  >  NC_000913/1992331‑1992802

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