Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	630,187	A→G	19.3%	G98G (GGA→GGG)	cstA →	carbon starvation protein involved in peptide utilization; APC peptide transporter family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	630,187	0	A	G	19.3%	145.2 / 3.6	70	G98G (GGA→GGG)	cstA	carbon starvation protein involved in peptide utilization; APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (13/40);  new base G (8/8);  total (21/49)
Fisher's exact test for biased strand distribution p-value = 6.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTG  >  NC_000913/629951‑630419                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             gCTCTGGGATACTTTGTTTTAAAGCGTGGGGAACAGAACAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATGTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGGCAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGc                                                                                                                                                                                                                                            <  2:365833/239‑1 (MQ=255)             ttttgtttAAACCGGGGGGAACAGAGCAACGCGCTGTGGGTTTGGGGGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTTTGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCtggt                                                                                                                                                                                                                               <  2:1067460/234‑1 (MQ=255)             tttGTTTTAAATCGGGGGGAACAGTGAAACGCTGTGTGGGTTGGGGTGGGGTCGGTCTGGATGTATCTGATCGCTTACCGTTTTTATGGGCTGTAGATCGCCAAAAATGTGCTGGCGGTTGGCCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCtggt                                                                                                                                                                                                                               <  2:915442/238‑1 (MQ=255)                 cTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTCGCTATCTGATCGCTTACCGTTTTTATGTGCGGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGGCTAGGTCGCCACCGGGCCGGCCAAAGGAAGGGCT‑GTAGTACGTGGCGCCGGACAAGAAAGTGCT‑CTTCCGTCCCCCTTTTGCGCGCAGTCGCGGAGGAG‑GGGCGCgggggggg                                                                                                                                                                                                                           >  2:192650/1‑239 (MQ=255)                 cTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCCACGCGTATGACGCCAGCGGTGCGCACTAACGACGGGCTGGGCCTAT‑TGTCGACGGGCCAGAAAGAGCT‑GTTGGGTCCCCCTTTTGCGGCCCTTGCCGGGGAGA‑GGCCGCTGgggggg                                                                                                                                                                                                                           >  2:207691/1‑239 (MQ=255)                  tttAACACGTGGGGAACAGATCAACGCGCTGTGGATTGTGTTGGCGTCGGTCTGTAACTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTgggg                                                                                                                                                                                                                           <  1:692339/238‑1 (MQ=255)                  tttAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTgggg                                                                                                                                                                                                                           <  1:207691/238‑1 (MQ=255)                      aaTCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTTACCCCACGCGTATGAGGCCCGCGGTGGGCCCTAACCACGGGCT‑GGGCTATGTGCCCGCGGGCACAAGAGTGCG‑TTTTGGTCGCCCTTTTTTGGCCACTGACCGGGGAG‑GGCCGCGGGTggggggggg                                                                                                                                                                                                                      >  2:586968/1‑234 (MQ=255)                      aaTCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGAACATTTGCCGACGGACAAGAAAGTGCC‑GTTCGGGCACCATTTTTCGGCCCTTGCCGCGGCAG‑GGCCCCCGGGGGGGCCGGt                                                                                                                                                                                                                      >  1:376820/1‑239 (MQ=255)                           tGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGTCCCGACTCGGATGACGCCAGCGGGGGGCCCTAACGCCGGGCT‑GGACCATGTGCCGACCGGCCAGCAAATGGT‑GTTCGGTTCCCGTTTTGTGGCCCTTGCCGGGGCAGCGACCGCTGGGGGGGcgcgcacgg                                                                                                                                                                                                                  >  2:654159/1‑231 (MQ=255)                                     aTCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGAACCGTATGACGCCCGCGGGGCGGCCTAACGACGGGCC‑GGACTCTGTGCCGCCGGCAAAGAAAGGGGC‑GTTCTGTCCCCCTTTTGCGGCCATTGCCGGGGCAG‑GACCGGTGGGGGGGCGGGGAC‑TGGG‑GGGGCaaaa                                                                                                                                                                                                     >  2:794728/1‑238 (MQ=255)                                          cgcgCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGGCCCGGAATACGCCAGCGGCGGGTCATAACGACGGGCC‑GGACTATGGGGCGACGGCACAGAAAGGGCTGGTCCGGTC‑CCCTTTTGCGGCCCTGTCCGGGGCAG‑GGCCGCAGGGGGGGGCGGGCC‑TGGG‑Ggcgcccacggggg                                                                                                                                                                                                >  2:584785/1‑230 (MQ=255)                                            cggtttgGTTTGTGGTGGCGTCGGTTTGTATCTATCTGACCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTAc                                                                                                                                                                                              <  1:377413/234‑1 (MQ=255)                                                       gtggtgGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGGATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGGC‑GGTCGGTCACCATTTTTGTGCCCTTGGCGGGGCAG‑GGCCGCCGGGTGGGCCGGGAC‑CAGCGGGCGGAAAGGGGCGACCCC‑CCggggg                                                                                                                                                                                   >  2:705195/1‑238 (MQ=255)                                                                      ggtaactgattGATCGGTTACCGTTTTTAGGGGCTGTATATCGCCAAAAAGGGGCTGGGGGTTGACCCGACGCGTATGCCGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGc                                                                                                                                                                   <  1:1087504/229‑1 (MQ=255)                                                                           ttattttaTCGCTTCCCGTTTTTTTGGGCTGGATATCGCCAAAAATGGGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTgggg                                                                                                                                                              <  2:830687/232‑1 (MQ=255)                                                                           ctatctGAGCGCTTACCGTTTTTAGGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGGATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTgggg                                                                                                                                                              <  1:459062/239‑1 (MQ=255)                                                                           ccaactGATCGGTTGCGGTTTTGATGGGCTGTATATCGCCAAGAATGTGCTGGAGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGGCTGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTgggg                                                                                                                                                              <  1:284128/235‑1 (MQ=255)                                                                              tctGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGGGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGGCCAAAAGGGGTACCTG‑CCGGGGGTGATTTGGCTGGCTGCTGgggggg                                                                                                                                                           >  1:721334/1‑239 (MQ=255)                                                                                tgtcgtgttttcGTTTTTTGGGGTTGTATATGCCCAAAAATGTGCTGGCGGTTGGCCCGACGCGTATGACGCCAGGGGTGCGCCATAACGACGGGCT‑GGACTTTGTGCCGGCGGACAAGAAAGTGCT‑GTTCGGTCACCTTTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGtt                                                                                                                                                         <  2:312934/228‑1 (MQ=255)                                                                                tgagcggtgtccGTTTTTATGGGGTTTATATCGCCAAAAATGTGCTGGCGGTTTGCCCGACGCGTATGACGCCAGGGGTGCGCCATAACGACGGTCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGtt                                                                                                                                                         <  2:197852/229‑1 (MQ=255)                                                                                tGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGGAC‑TGGC‑GGCGCCAAAGGGCTACCTT‑CCGGGGGTGGTCTTGCTTCCCGCTGgggggggt                                                                                                                                                         >  1:830687/1‑237 (MQ=255)                                                                                  aTCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGGATGACGCCAGCGGTGGGCCATAACGCCGGGCT‑GGACTTTTTGCCGCAGGGCAAGAAAGTGCC‑GTTCCGTCCCACTTTTGCGGGCATTGCCGGAGAAG‑GTGCCCCGGTGGGGCCGGGCCTTGGC‑GGC‑CACATTGGCGACCCGCCCGGGGGTGATCGGGCCGTCCGTCGGGGGGGTtc                                                                                                                                                        >  2:230663/1‑239 (MQ=255)                                                                                     gCTTACCGTTTTTTGGGGCTGTATATCGACAAAAATGGGCGGGCGGGTCGCCCGGCGGGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGc                                                                                                                                                    <  2:121715/239‑1 (MQ=255)                                                                                          ccttttttATGGGCTGTATGTCGCCAAAAATGGGCTGGGGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                                                                                                               <  2:944578/236‑1 (MQ=255)                                                                                          ccttttttAGGGGCTGTATATCGCCAAAAAGGGGTGGGGGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                                                                                                               <  2:627242/236‑1 (MQ=255)                                                                                          ccGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                                                                                                               <  1:655743/239‑1 (MQ=255)                                                                                          ccGTTTTTATGGGCTGGATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                                                                                                               <  1:789601/239‑1 (MQ=255)                                                                                          ccGGTTTTAGGGGCTTTATATCCGCAAAAATGGGCCGGCGGTTGCCCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                                                                                                               <  1:436029/239‑1 (MQ=255)                                                                                          aataagagacaggCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                                                                                                               <  1:402569/228‑1 (MQ=255)                                                                                               tttagtggCTGTATTACACAAAAAATGTGCTGGCGGTTGACCCGACGCGTAGGACGCCAGCGGTGCGTCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                                                                                                                                          <  1:445521/233‑1 (MQ=255)                                                                                                              tCGCCAAAAATGTGCTGGCGGTTGACCCGACGCTTATGACCCCGGCGGGGCGCCATAACGACGGGCT‑GGACTATGGGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTg                                                                                                                           <  2:1075290/239‑1 (MQ=255)                                                                                                              tCGCCAAAAATGGGCGTGCGGTTGGCCCGGCGCGTATGACGCCAGTGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTg                                                                                                                           <  1:574588/239‑1 (MQ=255)                                                                                                              tCGCAAAAGATTGGCTGGCGTGTGACTCGACGGGTATGGCGCCAGCGTGTCGCCATAACGGCGGGCT‑GGACTAGTGGCCGGCGGACAAGAAAGGGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTg                                                                                                                           <  2:401661/239‑1 (MQ=255)                                                                                                              accccAAAAATGGGCTGGGGGTTGCCCCGACGCGTAGGACGCCAGCGGTGCGCCATAACGGCGGGCT‑GGACTATGGGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTg                                                                                                                           <  2:721334/236‑1 (MQ=255)                                                                                                                           gCTGGCGGTTGACCCGGCGGGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCCCCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTa                                                                                                              <  1:700450/239‑1 (MQ=255)                                                                                                                                             gcgTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGGGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgc                                                                                             >  1:704669/1‑238 (MQ=255)                                                                                                                                               gTATGAGGCCAGCGGGGCGCCAGAACGGCGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgctc                                                                                           <  1:951353/238‑1 (MQ=255)                                                                                                                                               gTATGACGCCAGGGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgctc                                                                                           <  1:870794/238‑1 (MQ=255)                                                                                                                                                        cAGGGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAg                                                                                 <  1:648097/239‑1 (MQ=255)                                                                                                                                                         aGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCCAATGGGCTACCTG‑CCGGGGCTGATCTGGCGGCTCGCTGGGGGGGTTCGCGCCGGTGCGGGGCAGGATTTTCTGGGGGTGTTTGTTTTTACGCGGCGGGGCGGTCCGTCCCTgggggggg                                                                                >  2:574438/1‑234 (MQ=255)                                                                                                                                                                    ttatcggcgGGCT‑GGTCTTGGTCCCGACGGACAAGAAATGGTC‑GTTGGGTCCCCATTTTGCGGCCATGCGCGGGGCAG‑GTCCGCTGGGGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTCCCTG‑CCGGGGATGATCTTGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaagaa                                                                     <  2:672130/232‑1 (MQ=255)                                                                                                                                                                        cggtgggCT‑GGCCTATGTGCCGGCGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAgaga                                                                 <  1:911295/235‑1 (MQ=255)                                                                                                                                                                                ttggtttttgtgCCGTGGGCGAAGAAGGTGCT‑TTTGGGGCCCCTTTTGGGGGCCATGGCGGGAGGCG‑GGCCGGTGGGGGGGCGGGTAC‑TGGG‑GGCGCAAATGGGCTTCCTG‑CCGGGGAGGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGcc                                                          <  2:377641/231‑1 (MQ=255)                                                                                                                                                                                        tgtgCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGAAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCGGCTCCCTGGGGTGGGTCTCGCCGGTGCGGTGCAGGGTTTCATGGGGCTGTTGGGTTTTAGCCGCCGTGGCGGTGCCTCCCTGGGGGGGCCGGGCAAAGCACAGATGGGGCCAAacgccg                                                  >  2:116708/1‑239 (MQ=255)                                                                                                                                                                                             cGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGGGGTTCTCGCCGGTGCGGGGCAGGGTTTCATGGGGCTGTTTGTTTCTACGCGCCGGGACGGTCGCTCGCTGGGGGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCgggggg                                             >  1:853607/1‑237 (MQ=255)                                                                                                                                                                                             cGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCAGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGGTGATCTGGCTGCTCGCTGGGGGGGTTGTCGCCGGTGCGGTGCAGGGTTTCATGGGGCTGTTTGTTTCTACGCGCCGTGACGGGTGCCCGCTGGGGGGGCGGGGCAAAGAAAGGAGGGGGCCAACCGCGgggggg                                             >  2:439021/1‑237 (MQ=255)                                                                                                                                                                                                cGGACAGTAAAGTGCT‑TTTCGGTCACCATTTGGCGGCCTTTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTTCCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCCCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGAt                                           <  1:35865/238‑1 (MQ=255)                                                                                                                                                                                                cGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAAAGGGCTACCTG‑CCGGGGGTGGTCTGGCTGCTCGCTGGGGTGGTTTTCGCCGGGGCGGTGCCGGAGTTCATGGGGCCGTTTGTTTTTACGCGCCGGGCCGGGCGCCCGCGGGGGGGGCGGGGCAACGAAGGGATGGGGACAACCGCCggggggggt                                          >  2:159237/1‑234 (MQ=255)                                                                                                                                                                                                  gACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAACTGGGCTACCTG‑CCGGGGGTGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGGGCGGTGCAGGATTTCTTGGGGCTGTTTGTTTGTCCGGCCCGGGACGGTCGGCCCTTGGGTGGGGTTGTTCAAGAGAGGATGGGGACAACCCCCCGGGGGAtggg                                        >  2:465035/1‑236 (MQ=255)                                                                                                                                                                                                            tGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGGGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCCCGCTGGGGGGGCTGGTCAAAGAAGAGATGGGGCCAAACGCCGGGGGGATTGGGCGGGGGGcc                              >  2:977803/1‑239 (MQ=255)                                                                                                                                                                                                             ggt‑GTCCGGTCACCATTTGGCGGCCATTGCGGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGTGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:162362/235‑1 (MQ=255)                                                                                                                                                                                                             gCT‑GTTCGTTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:324619/237‑1 (MQ=255)                                                                                                                                                                                                             gCT‑GTTCGGTCCCCATTTTGCGGCCTTTGCCGGGGCAG‑GTCCGCTGGGGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:317162/237‑1 (MQ=255)                                                                                                                                                                                                             gCT‑GTTCGGTCCCCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGGGAAAATGGGCTCCCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:392727/237‑1 (MQ=255)                                                                                                                                                                                                             gCT‑GTTCGGTCACCATTTGGGGGCCATTGCCGGGGCAG‑GTCGGTTGGTGGGGCCGGTAC‑TGGG‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:289721/237‑1 (MQ=255)                                                                                                                                                                                                              gctggtgtgTCACTATTTTGCGGCCATTGCCGGAGCAG‑GTCCGGTGGGGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:654159/229‑1 (MQ=255)                                                                                                                                                                                                              cgttgtCGGTCCCCTTTTTGGGGCCATTGCCGGAGCAG‑GTCCGTTGGTGGGGCCGGTAC‑TGGC‑GGCGGAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                               <  1:794728/232‑1 (MQ=255)                                                                                                                                                                                                                          cctttttGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCa                 <  1:236549/236‑1 (MQ=255)                                                                                                                                                                                                                          ccttttgggggCCATTCCCGGGGGAG‑GGCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCa                 <  1:439021/230‑1 (MQ=255)                                                                                                                                                                                                                          ccATTTGGCGGCCATGGCCGGAGCAG‑GCCCGGTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCa                 <  1:234353/239‑1 (MQ=255)                                                                                                                                                                                                                           ttttttGCGGCCTTTGCCGGCGGGG‑GTCGGGGGTGGGGGGCGGTCG‑TGGC‑GGCGCAAATGGGCTACCGG‑CCGGGGATGATCGGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCGTGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCAt                <  2:561332/237‑1 (MQ=255)                                                                                                                                                                                                                           ctttttGGGGCCTTGGCCGGAGGAG‑GTCCGGTGGTGGGGCCGGTAC‑GGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCAt                <  1:584785/237‑1 (MQ=255)                                                                                                                                                                                                                                tGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTAAGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGGGATTGCGCTGGTGGCCTGCTTTATGATCATGGGCa           >  1:226959/1‑239 (MQ=255)                                                                                                                                                                                                                                 ggggCCGTGGCCGGGGGAG‑GTCCGGGGGGGGGGCCGGTAC‑GGGC‑GGGGAAAATGGGTGACCTG‑CGGGGGATGATCTGGCTGCCCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGAGGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCAt          <  2:843257/237‑1 (MQ=255)                                                                                                                                                                                                                                 gCGGCCATTGCCGGGGCAG‑GCCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAAGGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCAt          <  1:192534/239‑1 (MQ=255)                                                                                                                                                                                                                                     ccATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGGGAAAAGGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTAt       <  1:574438/238‑1 (MQ=255)                                                                                                                                                                                                                                       aTTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGGGATTGCGCTGGTGGCCTGCTTTTTGATCATGGGCATTATCCt    >  1:546742/1‑239 (MQ=255)                                                                                                                                                                                                                                         tcccGGGGCGG‑GTCCGGTGGGGGGGCGGGTGC‑TGGC‑GGCCGAAATGGGCCCCCGG‑CGGGGGATGATCTGGCTGCTCGTTGGGGTGGTTCTCGCGGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTg  <  2:546742/237‑1 (MQ=255)                                                                                                                                                                                                                                         tGCCGGGGCAG‑GGTCGTGGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTg  <  1:586968/239‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             GCTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCT‑GGACTATGTGCCGACGGACAAGAAAGTGCT‑GTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAG‑GTCCGCTGGTGGGGCCGGTAC‑TGGC‑GGCGCAAATGGGCTACCTG‑CCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTG  >  NC_000913/629951‑630419

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: ‑