Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,301,750	C→T	T190I (ACT→ATT)	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,750	0	C	T	97.5%	127.1 / NA	40	T190I (ACT→ATT)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/17);  minor base G (0/1);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGT  >  NC_000913/1301510‑1301976                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   aTGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAACAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTg                                                                                                                                                                                                                          >  2:259157/1‑251 (MQ=255)                         aCTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTc                                                                                                                                                                                                  <  2:159610/251‑1 (MQ=255)                              gtgtATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCg                                                                                                                                                                                              <  1:67790/250‑1 (MQ=255)                                         gaaaCGTTCTGTTTATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTATAAGTCACCGTAAGTTAACCCGTTCCGTACTTCTATAAATTACGTGTGCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGaaaaagtgaccc                                                                                                                                                                                  <  2:339528/249‑8 (MQ=255)                                          cAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGAccc                                                                                                                                                                                  >  1:41165/1‑250 (MQ=255)                                                  tGTTGCTCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCAGGTATTGATGAAATTCTTGAAGGGAAAACACCGATTACCGATCTCGGCGTGAAAGCTATTGATGAACACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATAGAGAAATTCGGCGAAAAATGGACCCAGCCTGGt                                                                                                                                                                          >  1:274226/1‑250 (MQ=255)                                                      gATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTCTTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCCGCCTGGAAatat                                                                                                                                                                     >  2:326711/1‑251 (MQ=255)                                                         ccGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGt                                                                                                                                                                  >  1:165507/1‑251 (MQ=255)                                                           gAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTc                                                                                                                                                                 >  1:368089/1‑250 (MQ=255)                                                                 tGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACAAATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTTCCAAAAGCCGCTAACGTGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACg                                                                                                                                                          >  2:365769/1‑251 (MQ=255)                                                                                        tGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGACTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTgg                                                                                                                                   >  1:158271/1‑251 (MQ=255)                                                                                          cAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGAGTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTAGAAAGTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCGAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGt                                                                                                                                 <  1:537890/251‑1 (MQ=255)                                                                                            atatGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCg                                                                                                                               <  2:221128/251‑1 (MQ=255)                                                                                             tatGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGt                                                                                                                              >  1:31430/1‑251 (MQ=255)                                                                                               tGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACgg                                                                                                                                                         >  1:355774/1‑222 (MQ=255)                                                                                                      atCGCCGGTATTGATGAAAATCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAAAACGTTCCGTACGTCTATAAATTACTTGTTCAACCATCAAGTTCACCGGTGCTAAAAGCCGCTATCGAGAAATTCTGCGAAAAATGGTCCCAGCCTTGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAaacgaac                                                                                                                      <  1:525877/250‑1 (MQ=255)                                                                                                          ccGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGAGCTCGGCGTGAAATCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATGACTGGTTCCCCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGGCACCAACGGTGCCGATACCTTAAAAGATTGGGTCGTAAACGAACGAATc                                                                                                                 <  2:295545/251‑1 (MQ=255)                                                                                                            ggTATTGATGAACTTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAACCGAACGAATCGt                                                                                                               >  2:426543/1‑251 (MQ=255)                                                                                                             gTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCCTCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCCGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGtt                                                                                                              >  1:583136/1‑251 (MQ=255)                                                                                                                 tgatgaAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAAGCGTTCtt                                                                                                           <  2:368089/250‑1 (MQ=255)                                                                                                                       aaTTCTTGAAGGGAAAAAACCGATTACCGATCTCGGAGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCATTCTTGAACGc                                                                                                     >  2:175415/1‑250 (MQ=255)                                                                                                                                  ggAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTTATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGagctcctcccgtcctac                                                                                         <  1:582960/251‑17 (MQ=255)                                                                                                                                          ccGATTACCGATCTCGGCGTGAAAGCTATTTATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGaaca                                                                                  >  2:203661/1‑250 (MQ=255)                                                                                                                                             aTTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGCACAAcgc                                                                              >  1:425499/1‑251 (MQ=255)                                                                                                                                                aCCGATCACGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATGCCTTAAAAGATTGGGTCGTAAAGGAAGGAATCGTTCTTGAACGCAGCCCGATCTACTGGAACAACGCGaa                                                                           <  1:510210/251‑1 (MQ=255)                                                                                                                                                   gATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAAAGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAAc                                                                        >  2:441297/1‑251 (MQ=255)                                                                                                                                                        cGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGtt                                                                    >  1:503838/1‑250 (MQ=255)                                                                                                                                                               aaaGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCa                                                            <  2:41165/251‑1 (MQ=255)                                                                                                                                                                aaGCTATTGATGATCACACATTAGAAGTCACCCTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCAGCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAAGTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGGGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCGTGAACGCAGCCCGACCGACGGGAACAACGCGGAAACGGTTATTAACCAg                                                           <  2:583136/251‑1 (MQ=255)                                                                                                                                                                  gCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCATCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGt                                                         <  1:51884/251‑1 (MQ=255)                                                                                                                                                                  gCTAGTGATGATCACACATTAGAAGTCTCCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCTCCCATCAATTTCACCGTTGCCAAAAGCCGCTATCGAGAAAGTCGGCGAAAAATGGACCCAGCCTGGTAATATCGGCACCAACGGTGCCTATACCGTAAAAGATTGGGTCTTAAACGAACGAAGCGTTCTTAAACGCAGCCCGTCCTACAGGAATAACGTGAAAACCGTTATTAACCAGGt                                                         <  2:524029/251‑1 (MQ=255)                                                                                                                                                                   cTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGt                                                         >  1:77523/1‑250 (MQ=255)                                                                                                                                                                     aTTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAAc                                                      <  1:38307/251‑1 (MQ=255)                                                                                                                                                                     aTTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATGCGGCGAAAAATTGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAAc                                                      <  2:287478/251‑1 (MQ=255)                                                                                                                                                                           gatCACACATTAGAAGTGACCTTAAGTGAACCCGCTCCCTACTTATATCACTTACTTGTTCACCCATCAATTTCACCCGTGCAAAAAGCCGCTATCGAGCAATTAGGCGAAACATGGACCCAACCTGGTAATATCGACACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGCACGAATCGTTCTTGCACGCAGCCCGACCTCCTGGAACCACGCGAAAACCGTTATCAACCAGGTAACCTAttt                                                >  2:548093/1‑251 (MQ=255)                                                                                                                                                                                      aGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCt                                      <  1:260050/250‑1 (MQ=255)                                                                                                                                                                                                 tAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTAcc                          >  1:197664/1‑251 (MQ=255)                                                                                                                                                                                                      gAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTACCCTATTTGCCTATTGCTTCTGAAGTTACCGATGt                     >  1:26754/1‑251 (MQ=255)                                                                                                                                                                                                                   aCTTCTATAAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGCTGTCAACCGCTACCGt        >  2:398827/1‑251 (MQ=255)                                                                                                                                                                                                                         ataAATTACTTGTTCACCCATCAATTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCATCCCGACCTACTTGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTTCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAgtggt  <  1:548093/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   ATGGCACGCCAGTCACAGCACAAGACTTTGTGTATAGCTGGCAACGTTCTGTTGATCCGAACACTGCTTCTCCGTATGCCAGTTATCTGCAATATGGGCATATCGCCGGTATTGATGAAATTCTTGAAGGGAAAAAACCGATTACCGATCTCGGCGTGAAAGCTATTGATGATCACACATTAGAAGTCACCTTAAGTGAACCCGTTCCGTACTTCTATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGT  >  NC_000913/1301510‑1301976

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