Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,050,472	G→A	L66L (CTG→TTG)	gcvT ←	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,050,472	0	G	A	95.2%	128.3 / NA	43	L66L (CTG→TTG)	gcvT	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/23);  minor base C (2/0);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATT  >  NC_000913/3050229‑3050716                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     tggCAATAATGGAAAGGTCATCACGAACGGAACTTTCGATGCCGAAAGGTTCAGAGTGTTGGGTAATCCCGGAGAGGTCTTTTTCGCGGGTGGCGGAGTACACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGGGGACGCTGAGGTCAT‑CTTTCACACAGCCAGAGGCATTCAAAATCCCAGAGTAGAGAGCTTTGCCGCTTTTGGCGAGCTTCGCCACCTCGTTCGCCAGCAGATAACGCAAAAACTcc                                                                                                                                                                                                                                               >  2:56249/2‑251 (MQ=255)        cagtgAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGCTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCGTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCTGGTg                                                                                                                                                                                                                                         <  2:499939/246‑1 (MQ=255)            gaaaggtactcccGAACGGGTATTTCGATGCCGAAAGGGGCAGCGTGTGGGGTAATCCAGGGGAGGTCGTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCAGTCAACAGCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCt                                                                                                                                                                                                                                    <  2:76111/239‑1 (MQ=255)             aaaGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTg                                                                                                                                                                                                                                   <  1:372936/251‑1 (MQ=255)                     tcaCGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCCGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGAGGAAGAAATCTTCAGTACAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGAAATGGGCTTTGCCGCTTTTGGTGCGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAGACTCCCGGGTGCGGCTGCCGCGagg                                                                                                                                                                                                                           >  2:65535/1‑249 (MQ=255)                     acaCGAACGGTAATTGTGATGCCGAAAGGTTCAGGGTGTTGGGTAATCCAGGAGAGGTCTTTTGCGCGGGTGGCAGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCTT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCTGCTGCCGCGAAg                                                                                                                                                                                                                           <  2:11807/250‑1 (MQ=255)                           aCGGTACTTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCCCAAAATCCCGGGAGCGGCTGCCGCGaacgccgac                                                                                                                                                                                                                     >  2:231867/1‑244 (MQ=255)                               tAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGGAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGAAAAGGGCTTTGCCGCTGTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGGGCGGCTGCCGCGAAGATCGACGATg                                                                                                                                                                                                                 >  2:311585/1‑251 (MQ=255)                                        tGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGCCGATGGTCATAtgt                                                                                                                                                                                                        >  1:90139/1‑251 (MQ=255)                                         ggcGAAAGGTGCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTTGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTATTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACAGCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCTTCTGCCGCGAAGATCGACGATGGTCATAtgt                                                                                                                                                                                                        <  1:174665/248‑1 (MQ=255)                                         gCCGAAAGGTTCAGCGTGTTGGGTAATCCGGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTCGACGAGGAGGTCAT‑CAATCACCCCGCCAGAGGCATTCCACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCACCAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATAtgt                                                                                                                                                                                                        >  2:282213/1‑250 (MQ=255)                                                       cGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCAAACCGCCAGAGGCATTCAACTTCCCCGAGTAAAGGGCTTTGCCGCTTTGGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGCTCGACGTTGGTCATATAGGACACATCAAACAtt                                                                                                                                                                                         >  2:416525/1‑251 (MQ=255)                                                                                       cgttggtcggatttcacaccGAGGCGGAAGAAACCTTCAGTAAAGTAGTTGTCTATGAGGTAAT‑CTATCCTACCGCCACAGGCATTCAACATCCCCGAGTAAAGGGCTTGTCCGTTTTTGGTGAGCTTCTCCACATCGTCCGTCAGCAGATAACGCAAAAGCTCCCGGGTGCGGCTGTCGCGAATATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGCTCGTACCGCATTATGTtc                                                                                                                                                            <  2:233245/229‑1 (MQ=255)                                                                                       cgGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTATTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTtcgtc                                                                                                                                                         <  1:189399/251‑1 (MQ=255)                                                                                       cgGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCAAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGAGAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACTTACCGCATGATGTtcgtc                                                                                                                                                         <  1:296732/251‑1 (MQ=255)                                                                                                 aGTTAAAAACGAGGGGGAAGAAATCTCTAGTAAAGTAGTGGACGATGAGGT‑ATGCTAGCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGaa                                                                                                                                               <  2:457455/251‑1 (MQ=255)                                                                                                    tAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCg                                                                                                                                            >  2:208167/1‑251 (MQ=255)                                                                                                     gacgtttttccgatcaaaaTTTTCAGTAAATTAGTATACGATGAGTTCAT‑CTATCACCCCGCCAGAGGCCTTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCATAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGGCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGt                                                                                                                                           <  2:356406/235‑1 (MQ=255)                                                                                                       caacGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACAGCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTaa                                                                                                                                         <  1:327229/251‑1 (MQ=255)                                                                                                          cGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGc                                                                                                                                      <  1:411452/251‑1 (MQ=255)                                                                                                                     aaaTCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGt                                                                                                                                           >  1:356406/1‑235 (MQ=255)                                                                                                                      aaTCTTCATGAAAGTAGTGGCCGATGAGTCCAA‑CTATCCCGCCGCCAGAGGAATGCAACATACCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACAGCGTTCGCCAGCAGATAACGCAAAAACGCCCGGGTGCGGATGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTACGGCATCGGTACGTACCGCATGATGTGCGTAGATTTGCGAACCGTAATGCAGCGGcatcata                                                                                                                          <  2:140880/251‑2 (MQ=255)                                                                                                                                            gATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGc                                                                                                    >  1:459090/1‑251 (MQ=255)                                                                                                                                               gAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCg                                                                                                   >  1:466779/1‑249 (MQ=255)                                                                                                                                                 ggTCAT‑CTATCACACCGCCAGAGGCCTTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAGACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCCTGCGAgcg                                                                                               >  1:401406/1‑251 (MQ=255)                                                                                                                                                    gtacctttccccccGCCAGGGGCATTCAACCCCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGGGCGTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGTCATCGGTACGTACCGCATGACGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCTACCATGCGAGcgcc                                                                                             <  2:477020/239‑1 (MQ=255)                                                                                                                                                               acCGCCCGAGGCCTTCAACATCCCCGAGTACCGGGCTTTGCCGCTTTTGGTGAGCTTCGCCCCATCGTTCGCCAGCCGATAACGCAAACCCTCCCGGGTGCGGCTGCCGCGACGATCGACGATGGTCATATGTGACCCATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCTCCATGCGAGCGCCGCACAGcttgg                                                                                  >  2:530280/1‑247 (MQ=255)                                                                                                                                                                    cAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGtgttgt                                                                             <  1:416525/251‑1 (MQ=255)                                                                                                                                                                           attcaaactCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGATCTTCGCCACACCGTTCGCCAGCAGATCACGCAAAAACTCCCGGGTGCGGCTGCCGCGCAGATCGCCGATGGTCATCTGTGACACATCAACCATTCCGGCATCGGTACGTACCGCCTGATGTTCGGCGATTTGCGAACCGTACTGCAGCGGCATCATCCAGCCGGGGAAATCCCACATGCGAGCGCCCCCACCCGGGTGTTGTTCGTAcc                                                                      >  2:337487/9‑250 (MQ=255)                                                                                                                                                                              cAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTGGGTCGTaaaa                                                                     >  2:55813/1‑246 (MQ=255)                                                                                                                                                                                 cATCCCCGAGTAACGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAgg                                                                  >  1:438998/1‑249 (MQ=255)                                                                                                                                                                                        gAGTAAAGGGCTTTGCCGCTTTTGGGGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGGGCGTCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCAGCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCtgttg                                                         <  2:242861/251‑1 (MQ=255)                                                                                                                                                                                         aGTAAAGGGCTTGGCCGCTTTTGGTGAGCCCCGCCACATCGCTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGTCGCGAAGATCGACGATGGTCATAGGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGAGGTCCGTCGATTTGCGAACCGTAATGCAGCGGCACCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCtgttgt                                                        <  1:337487/251‑1 (MQ=255)                                                                                                                                                                                                ggCTTTACCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGAGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCAAGCCGTGGCAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTAAAAGGGAGTCTGTTGTGCCATCt                                                 >  2:64926/1‑251 (MQ=255)                                                                                                                                                                                                                gTGAGCTTCGCCCCATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCAGCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATaa                                 <  1:548066/251‑1 (MQ=255)                                                                                                                                                                                                                   aGCTTCGCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCg                              <  1:89079/251‑1 (MQ=255)                                                                                                                                                                                                                         gCCACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTg                        <  1:89059/251‑1 (MQ=255)                                                                                                                                                                                                                          ccACATCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAACGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTg                        >  1:217618/1‑250 (MQ=255)                                                                                                                                                                                                                              aTCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCgt                                                                                    <  2:244307/186‑1 (MQ=255)                                                                                                                                                                                                                              aTCGTTCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCgt                                                                                    >  1:244307/1‑186 (MQ=255)                                                                                                                                                                                                                                   tCGCCAGCAGATAACGCAAAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACtttt               <  1:208167/250‑1 (MQ=255)                                                                                                                                                                                                                                           agATAACGCAAAAACTCCCGGGGGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGa      <  2:466779/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaCGCAAAAACCCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAAtt  <  2:407458/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CGGCAATCATGGAAAGGTCATCACGAACGGTAATTTCGATGCCGAAAGGTTCAGCGTGTTGGGTAATCCAGGAGAGGTCTTTTTCGCGGGTGGCGGAGTTAACAACGAGGCGGAAGAAATCTTCAGTAAAGTAGTAGACGATGAGGTCAT‑CTATCACACCGCCAGAGGCATTCAACATCCCCGAGTAAAGGGCTTTGCCGCTTTTGGTGAGCTTCGCCACATCGTTCGCCAGCAGATAACGCAGAAACTCCCGGGTGCGGCTGCCGCGAAGATCGACGATGGTCATATGTGACACATCAAACATTCCGGCATCGGTACGTACCGCATGATGTTCGTCGATTTGCGAACCGTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATT  >  NC_000913/3050229‑3050716

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