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	97.4%	128.2 / ‑4.6	38	L12P (CTC→CCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/19);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACAT  >  NC_000913/1894600‑1895079                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 gTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAgg                                                                                                                                                                                                                                       >  1:206920/1‑251 (MQ=255)   tGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGa                                                                                                                                                                                                                                      <  1:276982/250‑1 (MQ=255)     gCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCACTCGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGAGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGAcgc                                                                                                                                                                                                                                   >  1:31273/1‑251 (MQ=255)         gcagcaCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGCGAAGGTAAAACTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTGAGGACTGATGAAGACGTTATCTCCCGCTGTGCTTAC‑TTTACTCTGGCGTCAGGACGCCGCt                                                                                                                                                                                                                               >  1:152738/1‑251 (MQ=255)                    gCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTAttt                                                                                                                                                                                                                    >  1:374068/1‑251 (MQ=255)                    gCCACCAGCGCGATTGCCTGGCCGCTGCTCAACCCTTGTGCAATCGGTTCCTGGATCCGCTCGACCGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTTGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTTATGCAGAAGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTAttt                                                                                                                                                                                                                    >  1:268136/1‑251 (MQ=255)                      caccaGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTtct                                                                                                                                                                                                                  <  1:327848/251‑1 (MQ=255)                                 ttGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGc                                                                                                                                                                                                       >  1:352000/1‑251 (MQ=255)                                          cGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCg                                                                                                                                                                                              >  1:251452/1‑251 (MQ=255)                                          cGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCg                                                                                                                                                                                              >  2:82088/1‑251 (MQ=255)                                          cGCTGCTCATCCCTTGTGCCATCAGATCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCg                                                                                                                                                                                              >  1:84148/1‑251 (MQ=255)                                                tCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGCAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCg                                                                                                                                                                                        >  1:265554/1‑251 (MQ=255)                                                                cAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTcc                                                                                                                                                                         <  2:374068/250‑1 (MQ=255)                                                                            ccGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCGGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCTCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGatc                                                                                                                                                            <  1:420177/251‑1 (MQ=255)                                                                              gCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGAAGTTATCTCACGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGCCTATGCCGatcat                                                                                                                                                          >  2:381774/1‑251 (MQ=255)                                                                                 cGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGatcatc                                                                                                                                                         <  1:376775/249‑1 (MQ=255)                                                                                    cAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGtat                                                                                                                                                    <  1:82088/251‑1 (MQ=255)                                                                                                              agtgaAGGTAAACCTGCAAACATTTGTAACTCCTGCTAAATTTTTTGCGCTAATTATTTCATTCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgccga                                                                                                                          <  2:316999/251‑1 (MQ=255)                                                                                                                             gCAAACATTGTTAACTCCTGCTAAATGGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACttt                                                                                                           <  2:77515/251‑1 (MQ=255)                                                                                                                               aaaCATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTa                                                                                                          >  2:145513/1‑250 (MQ=255)                                                                                                                                aaCATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAAc                                                                                                        >  2:325912/1‑251 (MQ=255)                                                                                                                                 aCATTGTTAACTCCTGATAAATTGTTGGCGCTAATTATTTCATGCTCCCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAAcc                                                                                                       <  2:402866/251‑1 (MQ=255)                                                                                                                                   aTTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCAc                                                                                                     <  1:298889/251‑1 (MQ=255)                                                                                                                                                    aaTTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTCAAGAAACCg                                                                                    >  1:337058/1‑251 (MQ=255)                                                                                                                                                         ttgGCGCTAATTATTTCATGATACCCGGCACATAGCCAGTAGAGGCAGGAATGATGAACACGTTATCTCCCGCTATGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGACTTTTATTTCTCCAGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATACGTATAGCCGCTTTGATATTGTGGTCGCCGCTCCGATTTGCACTTTAACCACTTTCGGTACAGAAACCgttgtt                                                                               >  2:121941/1‑251 (MQ=255)                                                                                                                                                           ggCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgttgtt                                                                               >  2:258902/1‑249 (MQ=255)                                                                                                                                                                      tttCATGCTACCCGGCACATCGCCAGTAGAGTCAGGACTGATGCAGACGTTCTCTCCCGCTGTGATTTC‑TTTCCCCTGTCGTCAGGCCGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGCTATTGTGGTCGCAGAGCCGATTTGCACTTTACCCACTTTCGGTAAAGAAACCGTTGTTAGGGAAAGCGaa                                                                   >  2:253759/1‑250 (MQ=255)                                                                                                                                                                         cATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGaaaaa                                                                <  1:29826/250‑1 (MQ=255)                                                                                                                                                                                cccGGCCCATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGGGGGCGATGCTTGTACACGCCGGCGATGCCGATCATCCGTATAGCCTCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAAc                                                        <  1:211811/251‑1 (MQ=255)                                                                                                                                                                                     cacaTAGCCAGTAGAGTCAGGACTGATGAAGACGTTTTCCCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACATCGACCa                                                   <  1:325912/251‑1 (MQ=255)                                                                                                                                                                                              aGTAGAGTCAGGACGGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTTCCCTGGCGTCAGGACGCCGCTGAATTTTATTGCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGAcc                                          <  2:251452/251‑1 (MQ=255)                                                                                                                                                                                                 agagTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGc                                       >  1:238324/1‑251 (MQ=255)                                                                                                                                                                                                   agTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTGCCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCGGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTa                                     <  1:56285/251‑1 (MQ=255)                                                                                                                                                                                                       aGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAgg                                 >  2:51225/1‑251 (MQ=255)                                                                                                                                                                                                                gaAGAGGTTATCTCCCCTTGTGAT‑ACATGTACCCTGGCGTGAGGACGCAGCTGATGTTTATTTCGCCCGTTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACACCGGCTATGCTGATCAACCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGGGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCcag                         <  2:84148/250‑1 (MQ=255)                                                                                                                                                                                                                             tgcGCTGTGATTAC‑TTTACCCTGCCGGCGGTACGCCGCTAAATGTTATTTTTCCCGTTTCAGCCACCTGCTGTGGGCGATGCGTTTACACTCCGGCTATGCCGATCGACCGTATAGCCTCTTTGATATTGTGGTCGCCGATCCGATTAGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACGACGACCACTGATGACCCGCTACCGGGGCTCCAGCAGGTGATGGATc            <  2:155842/248‑1 (MQ=255)                                                                                                                                                                                                                                     gATTAC‑TTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACa   <  2:387774/251‑1 (MQ=255)                                                                                                                                                                                                                                      aTTAC‑TTTCCCCTGGCGTCAGGACGCCGCTTAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACAt  <  1:324692/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GTTGGCACGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTAC‑TTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACAT  >  NC_000913/1894600‑1895079

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