Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,894,839	T→C	L12P (CTC→CCC)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,894,839	0	T	C	82.6%	106.0 / 13.6	46	L12P (CTC→CCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (6/2);  new base C (20/18);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTC  >  NC_000913/1894594‑1895081                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ggAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGc                                                                                                                                                                                                                                               >  1:5643/1‑251 (MQ=255)                          gCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTGATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTa                                                                                                                                                                                                                         <  1:285067/248‑1 (MQ=255)                           ccaccaGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTtc                                                                                                                                                                                                                     <  1:114583/251‑1 (MQ=255)                                         gCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGcc                                                                                                                                                                                                        <  1:322312/250‑1 (MQ=255)                                           cTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCAcc                                                                                                                                                                                                     >  1:170389/1‑251 (MQ=255)                                                           ccTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCt                                                                                                                                                                                     <  1:341628/251‑1 (MQ=255)                                                             ttGTGCCATCAGTTCCTGGATCAGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGNttt                                                                                                                                                                                   >  2:444823/1‑251 (MQ=255)                                                                      cAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTGTATGGGTGAGGGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTAAGGACGCCGCTGAATTTTATTTCTCCCGCGTAAGCCACCTGCCGTGGGCGATGCTTTTACACTcc                                                                                                                                                                           <  2:390681/250‑1 (MQ=255)                                                                              ggATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGcc                                                                                                                                                                  >  2:150103/1‑251 (MQ=255)                                                                                    gCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGatcat                                                                                                                                                            <  1:170635/251‑1 (MQ=255)                                                                                       cgcaagttttttGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATGGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGTCGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATcc                                                                                                                                                          <  2:309167/243‑1 (MQ=255)                                                                                                      gctgTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGat                                                                                                                                          <  1:120077/251‑1 (MQ=255)                                                                                                          ttCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATtg                                                                                                                                      >  1:251884/1‑251 (MQ=255)                                                                                                                       gaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGATACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGACACCTGCCGTGGGCGATGCTTTTACAATCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgccgaccc                                                                                                                         >  2:382914/1‑248 (MQ=255)                                                                                                                           gTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGAtt                                                                                                                     >  1:120340/1‑251 (MQ=255)                                                                                                                                    cAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACttt                                                                                                             >  1:293449/1‑250 (MQ=255)                                                                                                                                      aaCATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGACGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCGGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGCTATTGTGTTCGCCTAGCCGTTTTGCACTTTAAc                                                                                                          >  2:260164/1‑251 (MQ=255)                                                                                                                                           tGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACtt                                                                                                     >  2:388762/1‑251 (MQ=255)                                                                                                                                            gTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACtt                                                                                                     <  1:308314/250‑1 (MQ=255)                                                                                                                                             ttAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTc                                                                                                   >  2:37187/1‑251 (MQ=255)                                                                                                                                              tAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGTTCGCCGAGCCGATTTGCACTTTAACCACTTTCg                                                                                                  >  1:159955/1‑251 (MQ=255)                                                                                                                                               aaCTCCTGTGAAATTTTTGTCGCTAATTATTTCATGCTCCCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTTATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGTTTACGCCACCGGCCGTGGGCGATGCTTTGACACTCCGGCTATGCCGATCATCCGTATAGCCTCTTTGTTAGTCTGGTCGCCGAGCCGATTTGCACTTTAGCCACTTTCg                                                                                                  <  2:470564/250‑1 (MQ=255)                                                                                                                                                   ccTGCTAAATTGTTGGCGCTAATTATTTCATGCTTCCCGGCACATAGCCATTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTaa                                                                                              <  2:262007/250‑1 (MQ=255)                                                                                                                                                     tGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTaa                                                                                              >  2:148917/1‑248 (MQ=255)                                                                                                                                                           aTTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGCCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGCGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgt                                                                                     >  1:268239/1‑251 (MQ=255)                                                                                                                                                                       aaTTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAATCCGTTGTTAGTGAAAg                                                                         <  1:388762/251‑1 (MQ=255)                                                                                                                                                                                aTGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACg                                                                >  2:399546/1‑251 (MQ=255)                                                                                                                                                                                aTGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACg                                                                >  2:328768/1‑251 (MQ=255)                                                                                                                                                                                          gCACATAGCCAGTGGAGGCAGGCATGTTGACGCCGTTATCTCCCGCTGTGATTACTTTACCCGGGCGTAGGGACGCCGCTGAATTTTATTCCCCCCGCTAAAGCCACCTGCCGTGGGCGATGCTTTTACCCTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGAcc                                                      <  2:159955/251‑1 (MQ=255)                                                                                                                                                                                            ccaTAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTCCCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCAc                                                    <  2:394958/250‑1 (MQ=255)                                                                                                                                                                                             caTAGCCTGAAGAGTCAGGGATGATGAAGACGTTATCTCCCGCTGTGATTAAGTTACCCTGGCGTCAGGGCGTAGATGAATTTTATTTCTCCCGCTTTAGCCACCTGCCGTGGGCGATGCTTTTACACACCGGCTAGGCCGATCATCCGTAGAGCAGCTTTGATATTGTGGTCGCCGAGCCGACTGGCACTTTAACCACTTTCGGTAAAGAAACCGTGGTTAGTGAAAGCGCAACACGCACAACGACCACt                                                   >  2:403262/1‑251 (MQ=255)                                                                                                                                                                                                agcaggTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCGGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACtgat                                                <  2:327129/246‑1 (MQ=255)                                                                                                                                                                                                      tAGAGTCAGAACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCg                                          >  2:34598/1‑251 (MQ=255)                                                                                                                                                                                                       agagTCAGGACTGATGAAGACGTTACCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCg                                          <  1:424502/250‑1 (MQ=255)                                                                                                                                                                                                          gTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGACAGCGAAAAACGCACAACGACCACTGATGACCCGCTAc                                      >  1:174234/1‑251 (MQ=255)                                                                                                                                                                                                                  tgatgaAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGCATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTCACCACTTTCGGTACAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGCCCACTGATGACCCGCTACAAGTGCTc                              >  1:473263/1‑251 (MQ=255)                                                                                                                                                                                                                  tgatgaAGACGTTATATCCCGCTGTGATTAATTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGAGAAAAACGCACAACGACCACTGATGCCACGCTACAGGTNCTc                              >  2:435714/1‑251 (MQ=255)                                                                                                                                                                                                                    atgaAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCca                            <  1:444823/251‑1 (MQ=255)                                                                                                                                                                                                                     tgaAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTcc                             >  1:77384/1‑249 (MQ=255)                                                                                                                                                                                                                       aaGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTgg                                                                                                                                                                                             <  2:338096/87‑1 (MQ=255)                                                                                                                                                                                                                       aaGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTgg                                                                                                                                                                                             >  1:338096/1‑87 (MQ=255)                                                                                                                                                                                                                              tcnctcCCGCTGTGATTATTTTACTCTGTCGTCAGTACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGGTGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGATTTGATATTGTGGTCGCCGAGCCGAGTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCt                   <  2:440355/247‑1 (MQ=255)                                                                                                                                                                                                                                tctcCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCAGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCAGTGGGCGATGCTTTTACCCTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTgg                 >  1:208158/1‑250 (MQ=255)                                                                                                                                                                                                                                     cGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATcgcg           >  1:455554/1‑251 (MQ=255)                                                                                                                                                                                                                                         gtgATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAg        <  1:37187/250‑1 (MQ=255)                                                                                                                                                                                                                                          tgATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGAc      >  1:282671/1‑251 (MQ=255)                                                                                                                                                                                                                                               actttaccCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTc  <  2:251884/250‑1 (MQ=255)                                                                                                                                                                                                                                               actttaccCGGGCTTCAGGACGCCGCTGATTTTTATTTCCCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCGCTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTc  <  2:396678/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GGAGTGGTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTC  >  NC_000913/1894594‑1895081

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