Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	299,490	T→C	E482G (GAA→GGA)	paoC ←	aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	299,490	0	T	C	97.7%	140.8 / NA	44	E482G (GAA→GGA)	paoC	aldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/18);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGA  >  NC_000913/299243‑299730                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  ccctccTGCGGTGGCTTCATGCTGCGTGGCGTGTCGGTGATGGTTGGTAATCTCTCCTTCTGCAAACGGCGACGGCGCAGGAACAAACCCGACTGCCGAGGCTATCAGTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGGGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGTCATCGCCGAGGTGAACCGCAATCTGCTCCAGCGGTACGCCAAGCATTTCCg                                                                                                                                                                                                                                                <  2:221834/247‑1 (MQ=255)          ggTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTGGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGGGGAGGGATTCGCGCCCCCATGTCCACC‑CGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGGACTCCAAGCATTCCCGCTGCCGTcg                                                                                                                                                                                                                                       >  2:135507/1‑250 (MQ=255)                  caggTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACGGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAg                                                                                                                                                                                                                                <  2:151009/247‑1 (MQ=255)                                cTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGTAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGcc                                                                                                                                                                                                                 >  1:408131/1‑251 (MQ=255)                                   cGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCgg                                                                                                                                                                                                               <  1:297582/250‑1 (MQ=255)                                     ggTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGc                                                                                                                                                                                                            >  1:310386/1‑251 (MQ=255)                                      gTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGcc                                                                                                                                                                                                           >  1:109916/1‑251 (MQ=255)                                           gTTGGTAACCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGt                                                                                                                                                                                                      <  1:674/251‑1 (MQ=255)                                               gTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCggg                                                                                                                                                                                                  >  1:96688/1‑250 (MQ=255)                                                  aTCTTGCCGTCGGCAAACGGCGAATGCTAAGGATCAAACACGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGAGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAAACGCAACCGGATCACGCGGTACGCCAAGCATTCCCGCTGGCGACTGGGCCAGAAGGGTGTAGCTGCCGGTGCCAATGACGGGCAt                                                                                                                                                                                               >  2:362602/1‑251 (MQ=255)                                                      tGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATg                                                                                                                                                                                              <  1:49605/248‑1 (MQ=255)                                                      tGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATg                                                                                                                                                                                              <  1:209612/248‑1 (MQ=255)                                                      tGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCGTCATACAGGCGGCGTAAATGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCTCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATg                                                                                                                                                                                              <  1:223530/248‑1 (MQ=255)                                                       gCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTcc                                                                                                                                                                                           <  1:365061/250‑1 (MQ=255)                                                          gTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGCCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGttt                                                                                                                                                                                       >  1:390190/1‑251 (MQ=255)                                                                                 gATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGt                                                                                                                                                                <  1:209951/251‑1 (MQ=255)                                                                                                ccGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGCGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTATCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCgg                                                                                                                                                 >  1:405508/1‑250 (MQ=255)                                                                                                  gagccAATCATTTCGCGAAGTTTCATTCGGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATGGTCCACC‑AGAACCCGCAGAAACCGGGAAACGGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTGTTGTTCGAGGTGAACCCGAGc                                                                                                                                               <  2:88931/247‑1 (MQ=255)                                                                                                    gggaaatctgcaggtgaccTTCATACAGGGGGTGTAAACGCCGGAGGGGGAGGTAGTCGCGCCCCATTG‑CCACCCAGAACCCGCGGAATACGGGATACGGGAATCGCCGTGGTGACCCGCAAACTGCGCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGGAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCAc                                                                                                                                             <  2:258194/233‑1 (MQ=255)                                                                                                      caataagtgcgcgAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGGCTCCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGATGCGCCGAGGTGAACCGCAACCTGCTCCAGGGGTACGCCATGCATCCCCGCTGCCTTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCTccc                                                                                                                                           <  2:91502/242‑1 (MQ=255)                                                                                                          cATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGAtt                                                                                                                                       >  1:114757/1‑251 (MQ=255)                                                                                                             ttggcgAAGCTCTATACAGGGGGGGTAAACGCCGGAGGTGTGGGTATCCGCTCCCCATTGCCCACC‑AGACCCCGCAGAAACCGGGAAACGGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGAtttt                                                                                                                                     <  2:152386/247‑1 (MQ=255)                                                                                                               ccacccgacccctttccgtcctgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATGGTCCACC‑AGAACCCTCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTcc                                                                                                                                  <  2:29832/229‑1 (MQ=255)                                                                                                                          aTACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGattatt                                                                                                                       >  1:223512/1‑251 (MQ=255)                                                                                                                            aCAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGG‑CCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCg                                                                                                                    <  1:344806/251‑1 (MQ=255)                                                                                                                             cAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCg                                                                                                                    <  1:362602/251‑1 (MQ=255)                                                                                                                                gcggcgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGaaa                                                                                                                 <  1:334065/251‑1 (MQ=255)                                                                                                                                  ggcgTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGc                                                                                                               >  1:220817/1‑251 (MQ=255)                                                                                                                                     gTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTcc                                                                                                                                  >  1:29832/1‑229 (MQ=255)                                                                                                                                               gaggtggaggtATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGCGGTGAACCCGAGCCCCCGATTTTTCCAGCAGATTATTGCGAACGCCCGCCGCA‑ACAcc                                                                                                 >  2:372312/1‑251 (MQ=255)                                                                                                                                                                 ccccATTATCCACC‑AGAACACGCAGAAACCGGGAAGCTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAACGGTACGACACGCATTCACGCTGCCGTCTGGGCCAGAATGGCGTAGCTGCCGGGGCCACTGTCGGACATGTCTGTTTCTACGGTAACGGGTCCATTTGGTTCGAGGGTAACCAGAGCACCCGATTTTTCCCGCGGATTATTGAGAAAGCCCGCCACA‑ACACCgatgccgaactgccgccc                                                                               >  2:422489/1‑234 (MQ=255)                                                                                                                                                                    cATTGTCAACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCTCAACCTGCTCCAGCGGTACGCAAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCAGCTTCTACGGTAGCGGGACCGTTTTGTTCGAGGGGAACCCAAGCACCCGATTTTTCCAGCAGATTATAGCGAAATCCCGCC‑CACACACCGTGGCCGACTATCC‑CTccccc                                                                           >  2:361181/1‑250 (MQ=255)                                                                                                                                                                                 gAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTg                                                                >  1:456600/1‑251 (MQ=255)                                                                                                                                                                                   aCCCGCAGAAACCGGAAAACTGGCATCGCCGGGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGACAGAATGGTGTAGCTGCCGGTGCCAATGGCGACCATGTCCGTTTCTACGGTAACGGTGACGGTTTGATCGAGGTGAACCAGAGCACCCGATTTATCCGGCAACTTATTGCGTATGCCCGCCACA‑ACACCGGGGCCGACtcacccctcccagccgcccccctgcc                                                              >  2:132160/1‑226 (MQ=255)                                                                                                                                                                                         aGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGCTTATTGCGAAAGCCCGCCGCA‑ACAACGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCggggg                                                        >  2:171876/1‑249 (MQ=255)                                                                                                                                                                                                     aCTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGGGCCGTTTTGTTCGAGGAGAACCCGAGCACCCGATTTTTCCAGCAGGTTATTGCGAAAGCCCGCCGCA‑ACACCGCGGCCGACTCGCCAATCCCCGTCGCGCCCCTCTCCGGGTGGGGCGTTGCgct                                            >  2:107948/1‑251 (MQ=255)                                                                                                                                                                                                         gAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTGCGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCtt                                        <  1:105124/251‑1 (MQ=255)                                                                                                                                                                                                                              aCCTGCTCCAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTcc                   >  1:319102/1‑251 (MQ=255)                                                                                                                                                                                                                              aCCTGCCCAAGCGGTACGCAACGCATTCCCGCTGCCGGTGGGGCCAGAATGGGGTAGCTGCCGGGGCCAATGCCGGTCATGTCCTTTTTTACGGTAACGGGGCCTTTTTGCTCGGGGTAACCCCGACCACCCGTTTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTcc                   <  1:361181/251‑1 (MQ=255)                                                                                                                                                                                                                                     ccAGCGGTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTgcgc            >  1:100818/1‑251 (MQ=255)                                                                                                                                                                                                                                           gTACGCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCAc      >  1:196024/1‑251 (MQ=255)                                                                                                                                                                                                                                               gCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTACCGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCg   >  1:286254/1‑250 (MQ=255)                                                                                                                                                                                                                                               gCCAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCCCCGGTGCGCAAGCACTCgg  >  2:396414/1‑250 (MQ=255)                                                                                                                                                                                                                                                ccAAGCATTCCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCCATGTCGGTCATGTCCGTTTCTAAGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGCGCACCCGATTTTTCCAGACGATTATTGCGACAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCAGATCCGGTGCGCAAGGACtaga  >  2:444086/1‑247 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  GCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATACAGGCGGCGTAAACGCCGGAGGTGGAGGTATTCGCGCCCCATTGTCCACC‑AGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCACCCGATTTTTCCAGCAGATTATTGCGAAAGCCCGCCGCA‑ACACCGTGGCCGACTAGCCACTCCCCGTCGCGCACCTGTCCGGGTGTGGCGTTGCGCTGCTTCCAGCCAAATTTATCCGCTCCGGTGCGCAAGCACTCGA  >  NC_000913/299243‑299730

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