Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,388,041	T→G	T50P (ACG→CCG)	aaeB ←	aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,388,041	0	T	G	98.2%	162.9 / ‑5.2	55	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (26/28);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCAC  >  NC_000913/3387807‑3388280                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  cgGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGAACGGCGGCAACAATc                                                                                                                                                                                                                                   >  1:185144/1‑251 (MQ=255)       aaCCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCt                                                                                                                                                                                                                             <  1:166254/251‑1 (MQ=255)       aaCCGGCCAGCCCCCAAGCATACTGGTTTTCTTTTCGTCCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACCCCATAATCATCAATATGGGGGCGCGAAGCATCGAAATGATGATCACCAGTCCGGCAATCCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGTCCGGCGGCAACATTCGCCGCt                                                                                                                                                                                                                             <  2:293966/251‑1 (MQ=255)                gCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCTGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTAAGTAcc                                                                                                                                                                                                                    <  1:121177/251‑1 (MQ=255)                    ccACGCATACGAGTTTTCAATTCGTACCAGCGAGGATATCCAGGTACAAAAACCTGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGccc                                                                                                                                                                                                                <  2:284978/251‑1 (MQ=255)                        gCATACGAGTTTTCTATTCGTACCCGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGCTCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGCTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCAACTCCCTCAGCAATGGCCTGACCGGCGGCAACAATCCCCCCTGTCAAGACCGccccg                                                                                                                                                                                                              >  2:373578/1‑247 (MQ=255)                         cATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCgg                                                                                                                                                                                                           >  1:293165/1‑250 (MQ=255)                           tACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGATCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTgg                                                                                                                                                                                                         >  1:27035/1‑251 (MQ=255)                                  tttCGATTCGTACCAGCGAGGATAGCCAGGTACAAAACCCGGCCCAGATACAGCAAACTAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATAACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAAGAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTTAGTACCGCCCAGCGTGGCGTTTc                                                                                                                                                                                                   <  2:54229/250‑1 (MQ=255)                                    tCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGACTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAcaccccccagtggcgcttcccg                                                                                                                                                                                                >  2:254959/1‑230 (MQ=255)                                        ttCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTgg                                                                                                                                                                                            >  1:340548/1‑251 (MQ=255)                                         tCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTgg                                                                                                                                                                                            <  1:335659/250‑1 (MQ=255)                                          cGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGa                                                                                                                                                                                           >  1:321638/1‑250 (MQ=255)                                                aGCGAGGATATCCAGGTACAAAAACCGGCCCAGATCAAGCACACCAGAATCATCAATAATGGAGCGCAAATCATCGCAATGATGATAACAAGTCCGGCAATACAGCCAATAACTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATCGCGCAAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCACAGCTGTCAGTAAAGCCCGACGAGGCGTTTCCAGCCGGCAGTGaa                                                                                                                                                                                     >  2:191074/1‑250 (MQ=255)                                                     ggATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTTCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaa                                                                                                                                                                               <  2:394387/251‑1 (MQ=255)                                                     ggATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaa                                                                                                                                                                               <  1:89696/251‑1 (MQ=255)                                                              ggTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcg                                                                                                                                                                      <  1:99784/251‑1 (MQ=255)                                                                  cAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                                          <  2:439922/195‑1 (MQ=255)                                                                  cAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                                          >  1:439922/1‑195 (MQ=255)                                                                      acccggccaatatccAGCACACCAGAATCATCAATAATGGGGCGCGGATCACCGCAATGAGGACCACCAGCCCGGCAATCCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGAGAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACa                                                                                                                                                              <  2:478417/237‑1 (MQ=255)                                                                             caaGATACATCACCCCAAAATCATCAATAATGGTGCGCGGAACATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCgg                                                                                                                                                       <  2:450852/249‑1 (MQ=255)                                                                                gATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGCTGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCCCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCCACAATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGGAGGGAAAGCCAACAAACAAGGCCAGTACAATGGCgggcg                                                                                                                                                    >  2:272435/1‑248 (MQ=255)                                                                                  tACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGCGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGCCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCACCGTGGCGTTTCCAGCTGGACATGAAAGCCAACAAACAGCGCCAGTACAATGGcagtactc                                                                                                                                                  >  2:203795/1‑244 (MQ=255)                                                                                    cAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCgggcccccg                                                                                                                                                >  2:101573/1‑244 (MQ=255)                                                                                            cAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCg                                                                                                                                       <  1:101573/251‑1 (MQ=255)                                                                                             aGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAACTGTGCCGATGATGCGCAAAAAGCCACGATAGCGCATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAGGGCGGTCGCCCG‑TTATACCGc                                                                                                                                      >  2:241838/1‑251 (MQ=255)                                                                                                  catcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAg                                                                                                                                  <  1:261351/250‑1 (MQ=255)                                                                                                   ataaataatGGTGCGCGGACCATCGCAATGATGATCACCAGCCCGGCAATACAGCCAATAAATTTGCCGTTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCGCCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCg                                                                                                                                <  2:47500/248‑1 (MQ=255)                                                                                                         aatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATg                                                                                                                          <  1:458285/251‑1 (MQ=255)                                                                                                                  gggaacatcGCAATGATGAACACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGc                                                                                                                 <  2:340879/246‑1 (MQ=255)                                                                                                                     atcatcGCAAGGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAAt                                                                                                              <  1:272435/251‑1 (MQ=255)                                                                                                                     atcaccGCAATGTTGCTCACCAGTCCATCAATACAGCCAATACATGGGCCCATTATGCTCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCGCCCGCAGCAAAGGCCGGACCTGCGGCAACAATCGCCGCTGTCAGTACCTCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATTGCGGTCGCCAG‑GGTTACCGCAATGCTGATATGTTGGTTAGCAAt                                                                                                              <  2:465526/251‑1 (MQ=255)                                                                                                                                   gatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTTCCGCCAAGCGGATATGTTGGTTAGCAATGGGGAAAAAAccc                                                                                                 >  1:202506/1‑250 (MQ=255)                                                                                                                                     tCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCATc                                                                                              >  1:22665/1‑251 (MQ=255)                                                                                                                                       aCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGt                                                                                            <  1:382244/251‑1 (MQ=255)                                                                                                                                          aGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAG‑TTTTACCCCAAAACGAAGATGGTGGTTTGTAATgggggaaaaacaccacggggt                                                                                         >  2:51157/1‑232 (MQ=255)                                                                                                                                                gaaaTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCTAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCa                                                                                   <  1:275819/249‑1 (MQ=255)                                                                                                                                                           aaTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACTGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑g                                                                       <  1:22559/251‑1 (MQ=255)                                                                                                                                                           aaTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑g                                                                       <  2:293165/251‑1 (MQ=255)                                                                                                                                                                     ccGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAAGACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGcc                                                             <  2:144459/251‑1 (MQ=255)                                                                                                                                                                        atgatgCGCAAAAAGCCACGATAGCGAGTAGCGCCAGAATACGGTCCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCAtt                                                          <  2:328072/251‑1 (MQ=255)                                                                                                                                                                                  aaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAG‑TTTTTCCGCCAAGCGAATATGTTGGGTAGCAATGGCGAAAATACCCATCGTGATTAACCAAACTC‑CCCACGGGGATGGGCCCTTTTACGga                                                   >  1:272137/1‑248 (MQ=255)                                                                                                                                                                                                 gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTAGGCAG‑GTTTAACGCAAAGCGAATATGTTGGTTAGCCATGGAGAAAATACCCAACGTGGTTTAACAAAATCACACA‑GGCGAATGACCATTTTACGGACAAACGGAACCTGGCttt                                 >  1:476548/1‑251 (MQ=255)                                                                                                                                                                                                 gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAACGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAAATTTTTCGGCAA‑GGGAATATGTTTGTTTGCAATTGAGAAAATAACCCTCGTTGTTAAACAAACTTACGCAAGGCGA‑GGGCCCTTTTACTGAAGAACGAAACCTGGAttt                                 >  1:312792/1‑251 (MQ=255)                                                                                                                                                                                                  aagaGCGCCAGAATTCGGTTCACCTCCCGCAGCAAAAGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCACCCAGCGTGGTGTTGCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCTAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                <  2:312792/248‑1 (MQ=255)                                                                                                                                                                                                     aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                             <  1:51157/251‑1 (MQ=255)                                                                                                                                                                                                      gcgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAACGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTAAGGAAGAACGAATCCTTGCTTtcgtcg                            >  1:427421/1‑251 (MQ=255)                                                                                                                                                                                                       cgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAAACAAACTCACGCA‑GGCGATGGGGCATTTTACGGAAGAACGACACCCGGCTTtcgcctc                           >  1:424930/1‑247 (MQ=255)                                                                                                                                                                                                        gcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTTCCGCAAAACGAATATGTTGGTTAGCAATGGGGAAAATACCCATCGTGATTAAACAAACTCACGCA‑GGCGGTGGGCCATTTTACGGAAAAACGAAACCTGGCTTTCGTcgcg                          >  1:347268/1‑251 (MQ=255)                                                                                                                                                                                                        gcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTcgcg                          <  1:425546/251‑1 (MQ=255)                                                                                                                                                                                                          cAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTAACGCCAAACGAATATGTTGGTTTGCAATGGAGAAAATACCCCTCGTGATTAAACCAACTCACGCA‑GGCGGTGGGCCATTTTAAGGGAGAACGGATCCTGGCTTTCGTAGCGAt                        >  1:371289/1‑251 (MQ=255)                                                                                                                                                                                                                  ggTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTAACGCAAAGCGGATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑AGCGATGGGCCCTTTTACGGGAAAACGAATACTGGCTTTAGTCGCGATATTGTTTg                >  1:228809/1‑251 (MQ=255)                                                                                                                                                                                                                    ttCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCACCCAACAGCGCCAGTACAATGGCGGTCGCCCG‑TTTTACCGCCAAACGAATATGTTGGTTAGCAATTGGGAAAATAACCATCCTGATTAACCAAACTCACACA‑GGCGCTGGGCCATTTTACGGAAGAAAGACTCCTTGCTTTCGTCGCGATCTTGTTTGcc              >  1:310368/1‑251 (MQ=255)                                                                                                                                                                                                                         cTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAACGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTgg         >  1:30405/1‑250 (MQ=255)                                                                                                                                                                                                                                agcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACcaccac  <  2:22665/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  CGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAG‑TTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCA‑GGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCAC  >  NC_000913/3387807‑3388280

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