Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,322,759	G→A	A63V (GCT→GTT)	trpE ←	anthranilate synthase subunit TrpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,322,759	0	G	A	97.8%	151.6 / NA	46	A63V (GCT→GTT)	trpE	anthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/27);  minor base C (1/0);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAGCTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTA  >  NC_000913/1322518‑1323000                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ggCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATgc                                                                                                                                                                                                                                          <  1:90547/251‑1 (MQ=255)     tctcGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATgcgc                                                                                                                                                                                                                                        <  1:240753/249‑1 (MQ=255)        cGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGCCAGGGTCACCTAAAACTGTAATgcgcagcgc                                                                                                                                                                                                                                   >  2:558539/1‑251 (MQ=255)        cGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATgcgcagcgc                                                                                                                                                                                                                                   >  1:482830/1‑251 (MQ=255)         gTTCTTCCTTCGGGACATTCAACATAGTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGTGCGTCTTCATCCAGCAGTGGACTCACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGGTATCCAGTAGTGCCAGGAGGGCTTCGCCGGTGCCGGAAAGGGCCGGGATTGGGACAGTGTCACCTAAAACTGTAAGgcgctgcgca                                                                                                                                                                                                                                  <  2:560943/251‑1 (MQ=255)                   cGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACca                                                                                                                                                                                                                        <  1:530269/251‑1 (MQ=255)                      tACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCGGGCAGGGCGTTATCCAGTAGTGCCCGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCt                                                                                                                                                                                                                            >  1:379285/1‑244 (MQ=255)                      tACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCGGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCt                                                                                                                                                                                                                            <  2:379285/244‑1 (MQ=255)                            cAACAGATGCGGCAATAAACGGAAAGCGGCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTGGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGc                                                                                                                                                                                                               <  2:462749/251‑1 (MQ=255)                                        taaTAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAatcat                                                                                                                                                                                                   >  2:323433/2‑251 (MQ=255)                                            aaaCGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATACGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCttt                                                                                                                                                                                               >  1:558695/1‑251 (MQ=255)                                                     gggtcgaaaaCCGAAAGGGAGCATAAGCGGGCGTCTCTATCTGGCAGTGGTCCTACAGGGGGGAAGCGCATGACACGGCAGTTTAGTGCTTGTTTACTTTCCACCTCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGCGGGCTTCTCCGTTGCCGGAAATTGCCTGGATTGTGACATTGTCACCTAAAACTGTAATGCGCATCACACTGTTTAGCAGCAGCAGGCTTTTCAAATCAGCTTTGCTGTAGat                                                                                                                                                                                      <  2:432103/245‑1 (MQ=255)                                                      cGTCAAAAGCCGAAAGGGCGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGCAGAGCAGCACACGGCCGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCAAGGAGGGCTTCGCCGTTGCCGGACAGTGGCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGAGCACTGTCTACCAGCCGCAGGCTTTTTAAATCATCTTTGCTGTCGatc                                                                                                                                                                                     >  2:269649/1‑250 (MQ=255)                                                            aaaCCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCg                                                                                                                                                                               >  1:44443/1‑251 (MQ=255)                                                               ccGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTGTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCTTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGAt                                                                                                                                                                            <  2:408831/251‑1 (MQ=255)                                                                cGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGAtt                                                                                                                                                                           >  1:487488/1‑251 (MQ=255)                                                                 gAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAACCTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTc                                                                                                                                                                          >  2:140508/1‑251 (MQ=255)                                                                     gggAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCAGCTTTGCTGTCGATATCTGCGGATTCcagc                                                                                                                                                                      <  2:114269/251‑1 (MQ=255)                                                                                     gTCTTCATCCAGCAGTGGACTGACCGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTACATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGAc                                                                                                                                                      >  2:276526/1‑251 (MQ=255)                                                                                         tCATCCAGCAGTGGACTGACAGGGAGGAAGCGCAGAACACGGCAGTTTGGGGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGGCTACCAGCAGCAGGCTTTTTAAATCAGCTTTGCTGTCGATAGCTGCGGATTCCAGCAGCAGCGTTGCCGGACg                                                                                                                                                     <  2:386025/248‑1 (MQ=255)                                                                                          cATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATcc                                                                                                                                                 >  1:31506/1‑251 (MQ=255)                                                                                            tCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATcccc                                                                                                                                               >  2:473092/1‑251 (MQ=255)                                                                                             ccAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTTCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCca                                                                                                                                              <  1:225190/251‑1 (MQ=255)                                                                                                          gACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGCTTCCAGCAGCAGCGTTGCCGGCCGCTcccccccccactggcgac                                                                                                                                 >  2:475971/1‑237 (MQ=255)                                                                                                               gggggAAGCGCAGCACACGGCAGTTTGGTGATTTTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACCCAACTGGTGAAAAAgc                                                                                                                            <  2:524598/251‑1 (MQ=255)                                                                                                                ggggAAGCGCAGCACACGGCAGTTTGGTGATTGTTCCCTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCTGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAgcg                                                                                                                           <  1:512319/251‑1 (MQ=255)                                                                                                                           gcaCACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTg                                                                                                                >  1:362481/1‑251 (MQ=255)                                                                                                                              cacGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTcg                                                                                                             <  1:411084/251‑1 (MQ=255)                                                                                                                                   cAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATa                                                                                                        <  1:475971/251‑1 (MQ=255)                                                                                                                                                         cctcacCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAATCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAg                                                                                   <  1:225296/247‑1 (MQ=255)                                                                                                                                                         ccACACCCTCAGGCATGGCGTTATCCAGTATGGCCATGAGTGCTTCGCCGTTGCCGGAAAGTGCCTGGATTATTACAGTGTCACCTAAAACGGTAATCGGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCTGGACGATTCCCACACAACTGGTGCAAAAGCGCGGGGGGAGTGTCGCGATAAGCGCCTTCGCAGGTTATCAg                                                                                   <  2:356339/250‑1 (MQ=255)                                                                                                                                                               ccGCAGACAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGCTAAGCGCCTTCGCAGGTTAGCAGTTCg                                                                               >  1:455255/1‑248 (MQ=255)                                                                                                                                                                        gggCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGtttt                                                                   <  2:455255/251‑1 (MQ=255)                                                                                                                                                                         ggCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGttttt                                                                  <  2:297404/251‑1 (MQ=255)                                                                                                                                                                          gCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCAGCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGtgtgtg                                                                 <  2:336988/251‑6 (MQ=255)                                                                                                                                                                            gTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACATAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTtgtg                                                               >  1:553915/1‑251 (MQ=255)                                                                                                                                                                             ttATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTgtgt                                                              >  1:447460/1‑251 (MQ=255)                                                                                                                                                                               aTCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTtgtgct                                                             <  2:362481/250‑3 (MQ=255)                                                                                                                                                                                tCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTg                                                           >  1:576949/1‑251 (MQ=255)                                                                                                                                                                                  cAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCa                                                         <  1:22351/251‑1 (MQ=255)                                                                                                                                                                                    gtagtGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCAtt                                                       <  1:473092/251‑1 (MQ=255)                                                                                                                                                                                                     tCGCCGTTGCCTGACAGTGCCTTGATTGTGACAGTGTCACCTAAAACTGTAATGAGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCAGCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGACGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCACAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAAttt                                        <  2:5923/249‑1 (MQ=255)                                                                                                                                                                                                                        ccTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTc                   <  2:437875/251‑1 (MQ=255)                                                                                                                                                                                                                        ccTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTc                   <  2:487488/251‑1 (MQ=255)                                                                                                                                                                                                                                 tgACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGt                                     >  2:42555/1‑224 (MQ=255)                                                                                                                                                                                                                                          cACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTa  <  2:44443/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAGCTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTA  >  NC_000913/1322518‑1323000

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