Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,301,992	A→T	N271Y (AAT→TAT)	oppA →	oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,301,992	0	A	T	97.7%	144.9 / NA	44	N271Y (AAT→TAT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/19);  minor base G (0/1);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGC  >  NC_000913/1301754‑1302231                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                aCCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGc                                                                                                                                                                                                                                     >  2:300781/1‑251 (MQ=255)   cGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCg                                                                                                                                                                                                                                   <  1:484608/251‑1 (MQ=255)                          aaaTTCGGCGAAAAATGGACCCAGCCTTGTAATATCGTCCCCAACGGTGCCTATACCGTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACCGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTGCCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGaa                                                                                                                                                                                                            <  1:499258/251‑1 (MQ=255)                                  cGAAACATGGACCCAGCCTGGTAATATCGTCCCCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAgaga                                                                                                                                                                                                    <  2:271954/251‑1 (MQ=255)                                            aCCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACg                                                                                                                                                                                           >  1:501086/1‑250 (MQ=255)                                                    tGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACg                                                                                                                                                                                  >  1:285947/1‑251 (MQ=255)                                                             cGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCAt                                                                                                                                                                         <  1:87126/251‑1 (MQ=255)                                                             cGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCAt                                                                                                                                                                         <  1:242340/251‑1 (MQ=255)                                                              gTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATa                                                                                                                                                                        <  1:436400/251‑1 (MQ=255)                                                                  ccAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCtg                                                                                                                                                                    >  1:94575/1‑251 (MQ=255)                                                                    aaCGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AAACGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTc                                                                                                                                                                                     <  2:480308/232‑1 (MQ=255)                                                                    aaCGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AAACGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTc                                                                                                                                                                                     >  1:480308/1‑232 (MQ=255)                                                                     aCGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGc                                                                                                                                                                 <  2:372569/251‑1 (MQ=255)                                                                        gTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTCCCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGACGTTCACGTTGATCCATACCTGTGCACt                                                                                                                                                              >  2:31441/1‑251 (MQ=255)                                                                          gCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCAGGTTTATCCATACCTGTGCACTTa                                                                                                                                                            <  2:51310/251‑1 (MQ=255)                                                                              ataCCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTa                                                                                                                                                         <  2:172906/250‑1 (MQ=255)                                                                              ataCCTTAAAAGAGTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGGAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGGGCACTTACTAt                                                                                                                                                        <  2:64720/251‑1 (MQ=255)                                                                                          aTNGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACTCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTGACCGATGTC‑AACCGCTACCGTAGTTGTGAAATCGACATGACGGATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCAGACCTGTGCACTTACTATTACGAAATTaa                                                                                                                                             <  2:483593/250‑1 (MQ=255)                                                                                                  gTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGa                                                                                                                                     >  1:105872/1‑250 (MQ=255)                                                                                                     aacgaacgAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAAcc                                                                                                                                 >  2:174062/1‑251 (MQ=255)                                                                                                                cGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACg                                                                                                                      >  2:323326/1‑251 (MQ=255)                                                                                                                        aaCGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGAAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGCTCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGAtgtgcgt                                                                                                              >  2:193627/1‑251 (MQ=255)                                                                                                                               cccGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGTTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACTTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTa                                                                                                       <  1:483913/251‑1 (MQ=255)                                                                                                                               cccGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATgtgcgtgtgcgt                                                                                                        >  1:298116/1‑250 (MQ=255)                                                                                                                                               aacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAgg                                                                                       >  2:391594/1‑251 (MQ=255)                                                                                                                                                     aaaaCCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGAt                                                                                                                                                                                                   <  1:92203/137‑1 (MQ=255)                                                                                                                                                     aaaaCCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGAt                                                                                                                                                                                                   >  2:92203/1‑137 (MQ=255)                                                                                                                                                           gTTATTACCCAGGTAACCTTTTTGCTTTTTCCTTCTGAAGTTCCCG‑TGTCCACCCGCTCCCTTGGTGGTAAAATCAACAGGACGTATAACAGCATGCCTATCGATTTTTCCCAAAAGCTAAAAAAAAAGACCCCGGACGAAGTTCTCGTTGACCCATTCCTGGGCATTTTCTTTTAAAAAATAAACAACAAAAAACCGCCACTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcgcg                                                                            <  2:308013/250‑1 (MQ=255)                                                                                                                                                            ttattaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGat                                                                          >  1:171761/1‑251 (MQ=255)                                                                                                                                                            tgtttacccAGGTAACCTATGTTCCTATTGCTTCTGAAGTTACCGATTTC‑AACCGCTACCGTATTGTGGAAATCAACATGACGTATAACACCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAAAGATCCCGGACGAAGTTCACTTTGATCCATCCCTGTTCACTTACTATTCCGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTCCCGCACTGATACTAGGTATGGACCGCGa                                                                           <  2:504551/243‑1 (MQ=255)                                                                                                                                                                     aGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCACCGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACAGCGATATCATTGtt                                                                 >  1:303705/1‑251 (MQ=255)                                                                                                                                                                     aGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGtt                                                                 >  1:93662/1‑251 (MQ=255)                                                                                                                                                                       gTAACCTATTGGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAAGCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAATTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATAAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTtaa                                                               <  2:93662/251‑1 (MQ=255)                                                                                                                                                                                  gCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGaa                                                      >  2:402654/1‑249 (MQ=255)                                                                                                                                                                                           ttCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTAAAACTAGGAATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGc                                            >  2:278127/1‑250 (MQ=255)                                                                                                                                                                                                   ttACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGccc                                   >  1:113647/1‑251 (MQ=255)                                                                                                                                                                                                       cGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCt                               <  1:6132/251‑1 (MQ=255)                                                                                                                                                                                                                cccGCTCCCGTAGTGGTGAAATCGACATGACGTAAAACAGCATGCCGATCGAATTTTTCCAGAAGCTAAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATCCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTac                       <  1:174062/250‑1 (MQ=255)                                                                                                                                                                                                                  cGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTacac                     >  1:10437/1‑251 (MQ=255)                                                                                                                                                                                                                       ccGTAGGGGTGAATTCGACATGACGTATAACAGCATCCCGATCAATTTGTCCCAGAAGCTAAAAAAAGAGATCCCGGACAAAGTTCACGTGGTCCCATCCCTGTGCATTTACTATTCGGAAATTAACAACAAAAAACCGCCATTCAACTATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATGGTTAATAAAGTAAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCa                 <  1:300781/250‑1 (MQ=255)                                                                                                                                                                                                                          tAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGt             >  1:121257/1‑251 (MQ=255)                                                                                                                                                                                                                           aGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGt             >  2:217340/1‑250 (MQ=255)                                                                                                                                                                                                                                tgAAATCGACATGACGTATAACAGAATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCACGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTCTGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACACAACCGTATACTg       >  2:87876/1‑251 (MQ=255)                                                                                                                                                                                                                                     tCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTCTACTGCTGgc  >  2:450336/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                ACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTC‑AACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGC  >  NC_000913/1301754‑1302231

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