Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,282,857	C→T	R998R (CGC→CGT)	narG →	nitrate reductase A subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,282,857	0	C	T	95.6%	133.5 / NA	46	R998R (CGC→CGT)	narG	nitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/19);  minor base A (2/0);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGG  >  NC_000913/1282628‑1283090                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          agagCGAGATGGATCTGCTGCGTATGCTCAAATACACCAAAGCGGAAGGTCCGGCGCACGGCCAGACGATGCTGACCACCGCAATTGCTGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGGCAGGTAGCCGTGAAAGCCGGGGCTGCCCTGACCGAATTTAACGGTCGTGAACATACGCATCTGGCGCTGAATAAAGAAGACTCGAAGATCCGCTTCCGAGATATTCAggaccagccacgc                                                                                                                                                                                                                      >  2:334631/1‑240 (MQ=255)      gagaTGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAAt                                                                                                                                                                                                                 >  2:552050/1‑251 (MQ=255)                     cGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGCCCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGAc                                                                                                                                                                                                  <  2:348047/251‑1 (MQ=255)                             cAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGGAACCAACGTTCAGGTAGCCGTGAAAGCCTGTGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCTCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGGTATTCAGGCACATCCTCGCAAAATTATCTCCAGCCCGACctggtct                                                                                                                                                                                           <  1:261152/250‑1 (MQ=255)                              aaCTACACCAAAGCGGAAGGTACGGCGAAAGGCCAGACGATGCTGAGCACCGAAATTGATGAGGCAAAGAGGATCCTGACACTGCCAACAGGAACCAAAGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTAAGCGAATTTACCGGACGTGACCATACGCATCTGGCGCTGCATAAAGAAGACGAGAAGATACGCTTCCGGAATAATCAGGCACAACCGCACAAAAttcaatacagcccggactgatctgg                                                                                                                                                                                         >  2:591759/1‑228 (MQ=255)                               aCTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAAGTAGCTCCAGCCCGACctgggctggt                                                                                                                                                                                        <  2:564501/251‑1 (MQ=255)                                          gCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGGGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGCCCTGGTCTGGTCTGGAAGATg                                                                                                                                                                              >  2:530219/1‑250 (MQ=255)                                            ggAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGAGACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCGGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGACCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCAAGCCCGGCCTGGCCTGGTCTGAAAGATGaa                                                                                                                                                                            >  2:160125/1‑250 (MQ=255)                                             gAAGGGCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAaca                                                                                                                                                                          <  2:428939/251‑1 (MQ=255)                                              aaGGTCCGGGGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGTTCTGGTCTGGAAGATGAacac                                                                                                                                                                         <  2:290820/251‑1 (MQ=255)                                                    cGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt                                                                                                                                                                   >  1:411911/1‑251 (MQ=255)                                                        gAAAGGCCAGCCGATGCTGAACCCCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                                                                                                                                               <  2:556533/251‑1 (MQ=255)                                                        gAAAGGCCAGCCGATGCTGAACAACGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGCGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc                                                                                                                                                                >  2:575941/1‑250 (MQ=255)                                                             gCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc                                                                                                                                                          >  2:108084/1‑251 (MQ=255)                                                               cAGCCGAGGCAGAACACCGCAAGTGATGCGGCAGAGATGATCCTGACACAGGCTCCGGAACCCAACGGTCAGGGAGCCGTGAAAGCCTGGGCGGCCCTGAGCGCATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCgg                                                                                                                                                        <  2:131152/251‑1 (MQ=255)                                                                     aTGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                                                                                                                                                  >  1:598865/1‑251 (MQ=255)                                                                           aaccccGCAATTGAGGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAAGTTACCGGTTGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCGTACAACGCCGGTTACACCAACGt                                                                                                                                            <  2:384560/247‑1 (MQ=255)                                                                                    aTTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCAGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGACCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAAAGCCGGTTACACCAACGTTCACGAGCt                                                                                                                                   >  1:126757/1‑251 (MQ=255)                                                                                               cgagaTGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGc                                                                                                                        <  1:167328/250‑1 (MQ=255)                                                                                                     gatCCTGACACTGGCACCGGAAACCAACGGCCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCGTCCGTGATATACAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACg                                                                                                                   <  1:335229/250‑1 (MQ=255)                                                                                                      atCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGct                                                                                                                 <  2:19942/251‑1 (MQ=255)                                                                                                       accTGACACTGGCACCGGAACCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGACGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGCTGACCACGTTTCTTACAACGCCGGTTCCACCAACGTTCACGCGCTGCTCCCCTGGCGTAccatc                                                                                                                >  2:141329/2‑247 (MQ=255)                                                                                                          tGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctctct                                                                                                             >  2:349343/1‑251 (MQ=255)                                                                                                              acTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtc                                                                                                         <  1:389426/251‑1 (MQ=255)                                                                                                                   cACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCACATGGCGTACGCTCTCTGgtcgtc                                                                                                      >  2:386956/1‑249 (MQ=255)                                                                                                                     ccGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCGGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCa                                                                                                  <  2:452307/251‑1 (MQ=255)                                                                                                                         aacccatcggtcaggtgtgcGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCAGACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACATCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCTGTTACACCGACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTg                                                                                              <  2:102794/232‑1 (MQ=255)                                                                                                                                            cGTGAAAGCCTGGGCTGCCCTGCGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAAAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGACGTCAGCAACTGTATCAGGATCACCAGTgg                                                                            >  1:540858/1‑250 (MQ=255)                                                                                                                                             gTGAAGGCCTGGGCTGCCCTGAGAGAATTTAACGGTCGTGACCATACGCATCTGGCGCTGCATAAAGAAGACGAGAAGATCCGCTTCCGTGCTATTCAGGCACAGCAGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTCCACCAACGTTCACGAGCTGATCCCATGGCGTACGCACTCTGGTCGCCAGCAACTGTATCAGGATCACCAGTGGAt                                                                          >  2:374554/1‑251 (MQ=255)                                                                                                                                              tGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATg                                                                         >  1:172492/1‑251 (MQ=255)                                                                                                                                               gAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGc                                                                        >  1:263889/1‑251 (MQ=255)                                                                                                                                                aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCATGCACAGCCGCTCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTTTATCAGGATCACCAGTGGATGCg                                                                       <  2:328576/251‑1 (MQ=255)                                                                                                                                                  aGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATCAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTg                                                                     >  2:64553/1‑251 (MQ=255)                                                                                                                                                    ccTGGGCTGCCCTGAGCGAATTTACAGGTCGTGACCATCCGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGAGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCTACGTTCCAGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGCTCACCAGAGGATGCGTGAt                                                                   >  2:125874/1‑251 (MQ=255)                                                                                                                                                       gggCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt                                                                  <  1:551067/249‑1 (MQ=255)                                                                                                                                                            gCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCACCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTg                                                            >  1:586891/1‑250 (MQ=255)                                                                                                                                                                 gAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGc                                                       <  1:569079/250‑1 (MQ=255)                                                                                                                                                                  aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAGAGCct                                                     <  2:473703/251‑1 (MQ=255)                                                                                                                                                                                    gACCATACGCATCTGGCGCTGCATACAGCAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc                                   >  2:326909/1‑251 (MQ=255)                                                                                                                                                                                     aCCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTAATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGACAGCCTGCTGGTTTATCGTccgccg                                  >  1:66179/1‑251 (MQ=255)                                                                                                                                                                                            gCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGac                            >  1:462496/1‑250 (MQ=255)                                                                                                                                                                                             cATCTGGCGCTGAATAAAGAAGCCGAGAAGACCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATGATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACTTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGAGGCGTGATTTCGGTGAAATCCTGCTGGTTTTTCGTCCGCCGATCGacac                          <  2:390484/251‑1 (MQ=255)                                                                                                                                                                                                         aatacagcacaCGAGAAGAACCGATTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCTACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGCAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTTGGTGa               >  2:521370/11‑250 (MQ=255)                                                                                                                                                                                                                 agaCGAGAAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg       >  1:202505/1‑250 (MQ=255)                                                                                                                                                                                                                  gaCGAGAAGATCAGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGACGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCCACAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGGGAt     >  2:360307/1‑251 (MQ=255)                                                                                                                                                                                                                     gagaAGATCCGCTTCCGTGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTCCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAgg  >  2:86125/1‑251 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          AGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGG  >  NC_000913/1282628‑1283090

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