Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	989,976	C→T	A127A (GCG→GCA)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,976	0	C	T	97.8%	142.3 / NA	46	A127A (GCG→GCA)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/21);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCATAAATACCCAGCAG  >  NC_000913/989737‑990221                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      gcCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGAGGGGTTGCGC‑AACGTCGGCCCGCGGTGAGAGATAGCGATGTACCATTTCCCTGATAAATGCCAGCAAAgg                                                                                                                                                                                                                                            >  1:32161/1‑251 (MQ=255)     gATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAAACGACTAATTTTCTTTCCAGCGGGTGGGGGGCCTTCGCCAACG‑CGGCCTGCGGGGAGGGACAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACt                                                                                                                                                                                                                                        >  1:177728/1‑251 (MQ=255)             ttGCTGGTGCGCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGCGAAGTAGACTAATTTGCTTTTCCGCGGGTTGAGGGGGTGTGC‑AACGTCGGGCTTCGGTGGGCGATAGCGCTGTACCATTTCACTTATAACCGCCCGCAAACGAAATTCCCaa                                                                                                                                                                                                                                  >  2:52264/1‑248 (MQ=255)              tGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCagag                                                                                                                                                                                                                                <  1:328501/250‑1 (MQ=255)                 tGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGAATACTTTGCTTTTCAGCGGGTCGAGGGCTTGAGC‑CACGTCGGCCTGCGGGGAGCGATAACGATGTACCCTTTTACTGATAACTGCCAGCAAAGGCACTTCCCAGAgggtg                                                                                                                                                                                                                            >  1:214444/1‑248 (MQ=255)                        cacaAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCAc                                                                                                                                                                                                                     <  2:1252/251‑1 (MQ=255)                         acaAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACg                                                                                                                                                                                                                    <  1:215979/251‑1 (MQ=255)                             aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGAGGACTTCACGCCAc                                                                                                                                                                                                                <  1:72787/251‑1 (MQ=255)                                  gATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AAAGTCGGCCTGCGGTGAGCGATAGCGATGTTCCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGcc                                                                                                                                                                                                           >  1:196616/1‑251 (MQ=255)                                           tGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAAc                                                                                                                                                                                                   <  2:13383/250‑1 (MQ=255)                                            gCAGACGCTTAACGATGGTTTCTTGTATTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGGGAAGTCGACTAATTTGCTCTCCCGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAAATGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACCCCCCGGGCCGCTTAAACg                                                                                                                                                                                                  >  2:259734/1‑250 (MQ=255)                                                     tAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCTGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAg                                                                                                                                                                                         <  1:216151/250‑1 (MQ=255)                                                      aaCGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCACACTCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGGGAAATCGACAAATTTGCTTTCCAGCGGGTCGAGGGCTTGCGC‑AAAGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCCGCCAAGGAACTTCCCAGAGGGTGACTTCACGCCACGGGCCGCTTAAACGAAAATcc                                                                                                                                                                                           >  1:336421/1‑247 (MQ=255)                                                               ttCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGGCGACTAATTTGCTTTCCCGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCCTTTCACTTATAACTGCCCGCAAAGGAACTTCCCCGAGGATGACTTCACCCCACGGGCCGCTTAAACGAAAATCCAGCTTGccctttt                                                                                                                                                                              >  2:127340/1‑246 (MQ=255)                                                                 cTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCCGAGGATGACTTCACGCCCCGGGCCGCTTa                                                                                                                                                                                                      >  1:273307/1‑225 (MQ=255)                                                                 cTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                                                                                                                                      <  2:273307/225‑1 (MQ=255)                                                                         tCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGaca                                                                                                                                                                     <  2:55894/250‑1 (MQ=255)                                                                             gagaAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGCCTAATTTGCTTTTCAGCGTGTCGAGGGCTTGCGC‑CAAGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGaa                                                                                                                                                                                                                          >  1:213800/1‑192 (MQ=255)                                                                             gagaAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGAc                                                                                                                                                                                                                          <  2:213800/193‑1 (MQ=255)                                                                               gaAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAc                                                                                                                                                              <  1:188593/251‑1 (MQ=255)                                                                                              gTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTaaa                                                                                                                                               <  2:274345/251‑1 (MQ=255)                                                                                                  cAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAATCGAATAATTTGCTTTCCAACGTGTCGAGGGCTTGCGC‑AAAGTCGGCCTGCGGGGAGCGGTAACGATGTACCATTTCACTGATAACTGCCCGCAAAGGAACTTCCCAGAGGATGAATTCACGCCCCCGGCCGCTTAAACGAATATCCAGCTTGCCCTTATCGGTGGACACGGTGACTTTTTCCGGGTTAAAGCgg                                                                                                                                           >  1:69383/1‑251 (MQ=255)                                                                                                        ccATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGTATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTAGCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTc                                                                                                                                      <  1:359114/250‑1 (MQ=255)                                                                                                                           cATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                                                                  <  1:233086/251‑1 (MQ=255)                                                                                                                            atatCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                                                                 <  1:101787/251‑1 (MQ=255)                                                                                                                                    aCCGGCGGGTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaa                                                                                                         <  2:143930/251‑1 (MQ=255)                                                                                                                                       ggcggTTAACGCCGAGCAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGAAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaa                                                                                                         >  1:309659/1‑248 (MQ=255)                                                                                                                                       ggcggTTAACGCCGAGAAGTCGACTAATTTTCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCGTaaaga                                                                                                      <  2:329023/251‑1 (MQ=255)                                                                                                                                       ggcggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaa                                                                                                         >  1:328610/1‑248 (MQ=255)                                                                                                                                       ggcggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaa                                                                                                         >  1:116518/1‑248 (MQ=255)                                                                                                                                         cggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaa                                                                                                     >  1:24825/1‑250 (MQ=255)                                                                                                                                                   cGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGACCTCGCGTCACCAGTTAAGATAGTCGGCCTTAAAGACAggcagggc                                                                                            >  2:91604/1‑249 (MQ=255)                                                                                                                                                     agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTGTCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACAGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaa                                                                                        <  2:74631/251‑1 (MQ=255)                                                                                                                                                       aaGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaa                                                                                        >  1:104427/1‑249 (MQ=255)                                                                                                                                                                aaTTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATc                                                                                                                                                                             >  1:82611/1‑155 (MQ=255)                                                                                                                                                                aaTTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATc                                                                                                                                                                             <  2:82611/155‑1 (MQ=255)                                                                                                                                                                              cGTGTCGAGGGCTTGCAC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa                                                               <  2:99589/251‑1 (MQ=255)                                                                                                                                                                                       ggCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGTCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAggg                                                      >  2:260875/1‑250 (MQ=255)                                                                                                                                                                                          ttGCGC‑AAGGTCGGCCTGCAGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACATCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGCGAGAAAGCCACCGATATTCATCATCCTGCAGGAGCAGGtgct                                                   >  2:272937/1‑251 (MQ=255)                                                                                                                                                                                                   cGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTGAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGcc                                           >  2:56048/1‑251 (MQ=255)                                                                                                                                                                                                            gCGGTGAGCGATAGCGATGTACCATTTCACGTATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCt                                    <  1:357072/249‑1 (MQ=255)                                                                                                                                                                                                                tGAGCGATAGCGATGTACCATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGAGATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCAc                              >  2:307259/1‑251 (MQ=255)                                                                                                                                                                                                                                  cATTTCACTGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCATaa             >  2:118442/1‑250 (MQ=255)                                                                                                                                                                                                                                          tGATAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCATAAATACCcagc    >  2:96442/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCATAAATACCcagcag  <  2:320961/251‑1 (MQ=255)                                                                                                                                                                                                                                             tAACTGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGGCTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCATAAATACCcagc    >  2:334412/1‑248 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      GCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGC‑AACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCATAAATACCCAGCAG  >  NC_000913/989737‑990221

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