Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,123,491	A→C	I105M (ATT→ATG)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,123,491	0	A	C	98.5%	232.0 / ‑3.5	67	I105M (ATT→ATG)	glcB	malate synthase G
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (34/32);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

TTTCCAGCGGTAGAGATTCATCGAGGAAACGCCGAACCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCAATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCC  >  NC_000913/3123262‑3123720                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      tttCCAGCGGTAGAGATTCATCGAGGAAACGCCGAACCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGccc                                                                                                                                                                                                                                 <  1:470719/236‑1 (MQ=255)                   cATCGAGGAAACGCCGAACCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCGGCGGCCCTGCCGGGCTGGGGGTTCCGCTTGCCATGCCCGGGGGTTCCAcctgccg                                                                                                                                                                                                                >  2:299656/1‑230 (MQ=255)                     tCGAGGAAACGCCGAACCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCAcgcg                                                                                                                                                                                                             >  2:302952/1‑236 (MQ=255)                               gCCGAACCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTg                                                                                                                                                                                                    <  2:515958/235‑1 (MQ=255)                                    aCCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCa                                                                                                                                                                                               <  1:675741/235‑1 (MQ=255)                                     cccAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTg                                                                                                                                                                                                    <  1:95146/229‑1 (MQ=255)                                        aGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCGGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGGTTCCACCGTCACGCGCTCCGGGTGCGGCAcccg                                                                                                                                                                                           >  2:676585/1‑233 (MQ=255)                                         ggCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGCTTTCGCTGTCCATGCCCGTGGTTTCCACCGTCCCGCTCTCCGGTTGCGGGACCAgg                                                                                                                                                                                          >  2:51955/1‑235 (MQ=255)                                                  cTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCGGCGGGATGATGTCGCTGCCGTATAACGGATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGt                                                                                                                                                                                 <  1:380901/235‑1 (MQ=255)                                                      tCACCGCGTTGCGGACCGTGGCCGCTGCCCAGCGCCCCTTCCTGCGGGAGTATGTCGCTGCCGTATCACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACg                                                                                                                                                                            <  2:523465/236‑1 (MQ=255)                                                        aCCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGGTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGACCAGTTTACGCa                                                                                                                                                                          >  1:764223/1‑236 (MQ=255)                                                          cgcgTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAgg                                                                                                                                                                        >  1:206029/1‑236 (MQ=255)                                                           ggttggCGGATCGTAGCCGTTGCCCACCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGGGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGa                                                                                                                                                                       <  2:80517/231‑1 (MQ=255)                                                           gcgTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGa                                                                                                                                                                       <  1:302952/236‑1 (MQ=255)                                                           gcgTGGCGCATCGTACCCGTTGACCATCGCCCTTTCCTGCGGGATGATGGCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGa                                                                                                                                                                       <  2:764223/236‑1 (MQ=255)                                                            cgTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACCGGCTCCGGTTGCGGCACCAGGTAGCCCAGTTTACGCAGGaa                                                                                                                                                                      >  1:360997/1‑236 (MQ=255)                                                                 cGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGCtt                                                                                                                                                                 >  1:595057/1‑236 (MQ=255)                                                                        tAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCAATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAgg                                                                                                                                                           <  1:722698/235‑1 (MQ=255)                                                                        tAGCCGCAGACCATCGCCCCTGCCTGCGGGATGATGTCGCTGCCGTAAAACGAATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAgg                                                                                                                                                           <  2:702301/235‑1 (MQ=255)                                                                        tAGCCGATGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGGGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAGCCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAgg                                                                                                                                                           <  1:203122/235‑1 (MQ=255)                                                                        gaGCCGCTGACAATCGCCCCTTCCGGCGGGATGATTGCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAgg                                                                                                                                                           <  1:112754/234‑1 (MQ=255)                                                                            cGCTGACCATCGCCCCTTCCTGCGGGATGATGGCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCt                                                                                                                                                      <  1:469871/236‑1 (MQ=255)                                                                              cTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCTGGCGTTCCTTGCCGGAACCATCAGCTGCGGCCCCGCCTGGCTTGTGGTTCCGCTGTCCTTGCCCGGGGTTTCCACTGTCATCCGTTCCGGCTTCGGCCCCACGTGTCCCAGCTCTGCCCGGACAGAATTATATGctcat                                                                                                                                                      >  2:742139/1‑230 (MQ=255)                                                                                   cATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTCTAg                                                                                                                                                            >  2:761648/1‑223 (MQ=255)                                                                                             tCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGccc                                                                                                                                     >  2:47736/1‑236 (MQ=255)                                                                                             tCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGGTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCCCTTTATCTTTTACCGGccc                                                                                                                                     >  1:681434/1‑236 (MQ=255)                                                                                                tGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTCTAGGCCGCGTTATCTCTTACAGGCCCCgg                                                                                                                                  >  2:692814/1‑236 (MQ=255)                                                                                                  cGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTGTAGGCCGCTTTATCTTTTACCGGCCCCgg                                                                                                                                  >  1:522806/1‑234 (MQ=255)                                                                                                   gggATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCgg                                                                                                                                  >  1:748585/1‑233 (MQ=255)                                                                                                        gatgTCGCTCCCGTAGACCGCATCGTACAGTGACCCCACGCGGGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATtgctg                                                                                                                           <  1:299656/235‑1 (MQ=255)                                                                                                             cGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGc                                                                                                                     <  1:700076/236‑1 (MQ=255)                                                                                                              gCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGc                                                                                                                     <  2:199281/235‑1 (MQ=255)                                                                                                              gCTGCCGTATAACGCATCGTACAGTGACCCCAAGCGGGCGTTCGCCGCGTCCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGc                                                                                                                     <  2:318162/235‑1 (MQ=255)                                                                                                                            cATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGCAAGATT‑‑TAGGCCGCTTTATCTTTTACCGGCCCAGGATTGCTGCGATGCCACTCATCAAgcg                                                                                                        >  1:76445/1‑232 (MQ=255)                                                                                                                            cATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCAACTCATCAAGCGCt                                                                                                      >  1:528656/1‑236 (MQ=255)                                                                                                                            cATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGACTGCTGCGATGCCActcctcc                                                                                                            >  2:234117/1‑229 (MQ=255)                                                                                                                             aTCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACtcatca                                                                                                            >  1:3459/1‑229 (MQ=255)                                                                                                                               cGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCAc                                                                                                                  <  1:62104/221‑1 (MQ=255)                                                                                                                               cGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCAc                                                                                                                  >  2:62104/1‑221 (MQ=255)                                                                                                                                gTACGGTGAGCCCCGGGGGGCGTTCGGCGGTTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGc                                                                                                    <  2:600883/234‑1 (MQ=255)                                                                                                                                     gTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAAt                                                                                              >  1:388060/1‑235 (MQ=255)                                                                                                                                       gAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCg                                                                                           >  2:494409/1‑236 (MQ=255)                                                                                                                                        aGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGaaa                                                                                                                                                                     <  2:282030/161‑1 (MQ=255)                                                                                                                                        aGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGaaa                                                                                                                                                                     >  1:282030/1‑161 (MQ=255)                                                                                                                                        aGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATg                                                                                             >  1:600883/1‑233 (MQ=255)                                                                                                                                        aGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCg                                                                                           >  1:37123/1‑235 (MQ=255)                                                                                                                                        aGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGa                                                                                          >  2:498389/1‑236 (MQ=255)                                                                                                                                                cgagcgTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATcg                                                                                       >  1:392719/1‑231 (MQ=255)                                                                                                                                                cgagcgTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGGGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATcg                                                                                       >  1:348692/1‑231 (MQ=255)                                                                                                                                                  agcgTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGGGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTg                                                                                >  1:283688/1‑236 (MQ=255)                                                                                                                                                    cgTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGcc                                                                              <  1:498389/236‑1 (MQ=255)                                                                                                                                                                   gcggttggCGGTCTTTCATTCCCGGAACCCCCGGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGAttt                                                               <  2:360997/229‑1 (MQ=255)                                                                                                                                                                         aGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGCAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTCAATCCGACCTCGGTCTTCAGTCAACTGACACTTTTCTgcg                                                         >  2:387449/1‑234 (MQ=255)                                                                                                                                                                                     ttGCCGACACCCCCGGCTGGGGCCCGGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGa                                             <  2:3459/236‑1 (MQ=255)                                                                                                                                                                                          ggAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGAt                                        <  1:494409/236‑1 (MQ=255)                                                                                                                                                                                          ggAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGAt                                        >  1:232639/1‑236 (MQ=255)                                                                                                                                                                                                        cGGCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGCAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGTTCATGAACGATACCATCCAACTTgt                          >  2:444152/1‑235 (MQ=255)                                                                                                                                                                                                          gCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGtt                                                                                   <  1:508013/177‑1 (MQ=255)                                                                                                                                                                                                          gCCCCGCCTGGCTGGTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGtt                                                                                   >  2:508013/1‑177 (MQ=255)                                                                                                                                                                                                                    gctggTGATTTCGCTGGCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAAcgccg              <  2:206029/236‑1 (MQ=255)                                                                                                                                                                                                                      tggTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTTCCGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCa            <  1:74952/236‑1 (MQ=255)                                                                                                                                                                                                                      tggTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCa            <  1:47736/236‑1 (MQ=255)                                                                                                                                                                                                                      tggTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCa            <  2:37123/236‑1 (MQ=255)                                                                                                                                                                                                                      tggTGATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCa            <  1:678375/236‑1 (MQ=255)                                                                                                                                                                                                                          gATTTCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGt        >  2:325578/1‑236 (MQ=255)                                                                                                                                                                                                                            tttCGCTGTCCATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTc       <  1:387449/235‑1 (MQ=255)                                                                                                                                                                                                                                ccTGTCCAGGCCCGTGGTTTCCACCGTCACGCTCTCCGGTTGCGGCACCAGGTAGCCAAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGcc  <  2:481069/235‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      TTTCCAGCGGTAGAGATTCATCGAGGAAACGCCGAACCCAGGCGATAACCTGCTCACCGCGTTGCGGATCGTAGCCGCTGACCATCGCCCCTTCCTGCGGGATGATGTCGCTGCCGTATAACGCATCGTACAGTGAGCCCCAGCGAGCGTTCGCCGCGTTCAGCGCGTAGCGGGCGTTCATTGCCGGAACCACCAGCTGCGGCCCCGCCTGGCTGGTGATTTCGCTGTCAATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCC  >  NC_000913/3123262‑3123720

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