Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,804,133	C→T	S211F (TCT→TTT)	nrdF →	ribonucleoside‑diphosphate reductase 2 subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,804,133	0	C	T	97.9%	145.0 / ‑6.5	47	S211F (TCT→TTT)	nrdF	ribonucleoside‑diphosphate reductase 2 subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/18);  minor base A (0/1);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATCTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGAT  >  NC_000913/2803907‑2804374                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  cAGATTATTCAGCAACATTATCGCGGTGATGTTCCGCTGAAAAAGAAAATCGCCATTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGtgaagagttg                                                                                                                                                                                                                             <  2:125233/249‑1 (MQ=255) cAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGtgaagagttga                                                                                                                                                                                                                            <  1:6729/250‑1 (MQ=255)           agtaaCATTATCGCGGTGATTATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTATGGAATCTTTTTTGTGCTATTCCGGTTTCTGGGTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATATCAGGTACAACGTGAAGAGTTGAAGAGTTTCGcc                                                                                                                                                                                                                 <  2:48299/248‑1 (MQ=255)            gCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCt                                                                                                                                                                                                                >  1:278425/1‑251 (MQ=255)              aaCATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTc                                                                                                                                                                                                              >  1:206251/1‑251 (MQ=255)                  ttATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACAAATACCGCGGACCTGATACGTCTGATTATCCGCGATGAAGCAGTCCACTGTTACTACATCGGCTATAAATATCAGAAAAACATGGAAAAATTATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGAtt                                                                                                                                                                                                          >  2:441734/1‑251 (MQ=255)                             atgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGAATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGa                                                                                                                                                                                                >  1:150454/1‑250 (MQ=255)                              tgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATttgtt                                                                                                                                                                                                      >  2:173435/1‑243 (MQ=255)                              tgatCCGCTGAAAAAGAAAAGCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTAGCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATttgtt                                                                                                                                                                                                      <  1:173435/243‑1 (MQ=255)                                   cGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACg                                                                                                                                                                                                                                       <  2:364898/205‑1 (MQ=255)                                   cGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACg                                                                                                                                                                                                                                       >  1:364898/1‑205 (MQ=255)                                     cTGAAAAAGCAACTCGCCTGTGTGATTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCAGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATACGTCTGATTATCCGCGATGAAGCAGACCACGGTTACTACATAGGATATAAATATCAGAAAAACATGAAAAAGATATTTCTGGGACAACGTGAAAAGTTGAAGAGTTTCGCATTCGATTTGTTGCTGGAACTCTACg                                                                                                                                                                                        >  2:109750/1‑250 (MQ=255)                                         aaaaGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGAt                                                                                                                                                                                                           >  2:85125/1‑227 (MQ=255)                                         aaaaGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGAt                                                                                                                                                                                                           <  1:85125/227‑1 (MQ=255)                                                      cAGTGTGCTTCTTGAATCTTTTTTGTTCTATTCCGGTCTCTGGCTGCCGTTGTATTTTTCCAGACGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGCTTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGCCAAACATGGAAAAGATATTTCTGGGACAACGTAAAGAGTTGAAGAGTTTCGTCTTCGATTTGTTGCTGGAACTCTCCGCCAACGAGTTGCAATac                                                                                                                                                                       >  2:214232/1‑250 (MQ=255)                                                            gtTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATg                                                                                                                                                                <  2:206251/251‑1 (MQ=255)                                                             tttCTTGAGTCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGCTGTATTTTTCCAGCCTCGGAAAGCTGCCCCATACCGCGGACCTGATCCGTCTGATTATCCGCGATTAAGCAGTCCACGGTTACTACATAGTCTATAAGTATCAGAACAACATGGAACAGCTATTTCTGGGCCAACGTGAAGCGTTGATGAGTTTCGCCTTCGATTTGTTGCTGGAACTATACGCCAGCGATTTTCAATacacagctg                                                                                                                                                                >  2:469499/1‑245 (MQ=255)                                                                 ttGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAACTATCAGAAAAACATGGAACAGATATTTCTGGGCCAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTg                                                                                                                                                           >  2:37434/1‑251 (MQ=255)                                                                     aTCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACg                                                                                                                                                       >  1:68425/1‑251 (MQ=255)                                                                     aTCTTTTTTGTTCTATTCCGGTTTCTGGCTGACGATGTATTTTTCCAGCCGCGGAAAGCTGACCACTACCGCGGACCTGATCCGTCTGATTCTCCGCGATGAAGCAGTCCACGGGTACTACATAGGCTATAAATATCAGAAAAACATGAAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGATTCGCCTTCGCTTTGTTGCAGAAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTAc                                                                                                                                                        >  1:384200/1‑250 (MQ=255)                                                                           tttGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTACAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGaa                                                                                                                                                  >  1:195948/1‑250 (MQ=255)                                                                             tGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAAcc                                                                                                                                               >  1:324542/1‑251 (MQ=255)                                                                                        ggatgcTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGAAGAAGAAGTCCACGGTTACTACATAGGCTCTAAATATCAGAAAAACATGGAAAACATATTTCTGGGACAACGTGAAGAGATTAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATAAACCGATGCGTTGTACGCCGAAAACCCGAGGGCTgg                                                                                                                                    >  2:420981/6‑250 (MQ=255)                                                                                         gTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTTAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGAc                                                                                                                                   >  1:451206/1‑251 (MQ=255)                                                                                                    ccGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGc                                                                                                                        <  1:285143/251‑1 (MQ=255)                                                                                                         gTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACCCCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACTATGTGAAAGCGTTtc                                                                                                                   <  1:370135/251‑1 (MQ=255)                                                                                                                tCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGAAGTCCACGGTTACTCCATAAGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGCGTTGCACTACACCGGTGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGtt                                                                                                             >  2:263039/1‑250 (MQ=255)                                                                                                                 ccAGCCGCGGAAAGCTGACCAATTCCGCGGACCTGATCCGTCTGAGTATCCGCGATGAAGCAGTCCACGGTTACTATATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACGCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCTAAACCCCGTGGGCTGACGATGTTAAAGCGTTTCTCTGTTAc                                                                                                           <  1:372584/251‑1 (MQ=255)                                                                                                                   aGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACa                                                                                                          >  1:281803/1‑250 (MQ=255)                                                                                                                                             ggACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTgg                                                                               >  1:316953/1‑251 (MQ=255)                                                                                                                                             ggACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTgg                                                                               >  2:45369/1‑251 (MQ=255)                                                                                                                                                   gATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCACTACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGCATCTGGGCTACg                                                                         >  1:384083/1‑251 (MQ=255)                                                                                                                                                          cggATTATCCGCGATTAAGCAGTCCACGGTTACGACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGACTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACTCCGAAACCCCGTGGGCTGACGATGTGAAAGCTTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGt                                                                   <  2:446581/248‑1 (MQ=255)                                                                                                                                                                  ccGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCg                                                          <  1:7307/251‑1 (MQ=255)                                                                                                                                                                       aTGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGCTGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGa                                                      >  2:40712/1‑250 (MQ=255)                                                                                                                                                                        tGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAATAGTTTCGCCTTCGACTTGTTGCTGGAACTCTACGACATCGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGCAATCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACTAACCGTTATTTCCCGCAGaa                                                     <  2:163520/250‑1 (MQ=255)                                                                                                                                                                           agtaatCCACGGTTTCTACATAGGCTATAAATAGCGGAAAAACATGGAAAAGAAATTTCTGGGACAACGGGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACCCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTTACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTGATTTCCCGCAGAAATgg                                                 <  1:78602/246‑1 (MQ=255)                                                                                                                                                                             cAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCAGGGCTACGAACCGTTATTTCCCGCAGAAATGGCg                                               >  2:439619/1‑251 (MQ=255)                                                                                                                                                                                  cACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGt                                          <  1:330222/251‑1 (MQ=255)                                                                                                                                                                                             cATAGGCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAATGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCaa                               <  1:127776/251‑1 (MQ=255)                                                                                                                                                                                                  gCTATAAATATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTc                          <  1:45369/251‑1 (MQ=255)                                                                                                                                                                                                     ataAATATCAGAACAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTcgcc                       >  1:135011/1‑251 (MQ=255)                                                                                                                                                                                                        aaTATCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATACGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCcgcg                    >  1:501329/1‑251 (MQ=255)                                                                                                                                                                                                         atatCAGAAAAACATGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGACTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCgcgc                   >  1:255034/1‑251 (MQ=255)                                                                                                                                                                                                                    aCATTGAAAAGATATTTCTGGGACAACGGGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCTCCGAAt        <  2:195948/251‑1 (MQ=255)                                                                                                                                                                                                                      aTGGAAAAGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGc      >  2:309302/1‑251 (MQ=255)                                                                                                                                                                                                                          aaaaGATATTTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGAt  <  1:40712/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CAGATTATTCAGCAACATTATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAGGCTATAAATATCAGAAAAACATGGAAAAGATATCTCTGGGACAACGTGAAGAGTTGAAGAGTTTCGCCTTCGATTTGTTGCTGGAACTCTACGACAACGAGTTGCAATACACCGATGAGCTGTACGCCGAAACCCCGTGGGCTGACGATGTGAAAGCGTTTCTCTGTTACAACGCCAATAAGGCTTTGATGAATCTGGGCTACGAACCGTTATTTCCCGCAGAAATGGCGGAAGTGAATCCGGCAATCCTCGCCGCGCTTTCGCCGAATGCCGAT  >  NC_000913/2803907‑2804374

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