Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,998	Δ1 bp	noncoding (100/109 nt)	spf →	small regulatory RNA Spot 42

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,995	0	G	.	100.0%	150.1 / NA	36	noncoding (97/109 nt)	spf	small regulatory RNA Spot 42
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (22/14);  total (22/14)

TTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATT  >  NC_000913/4049920‑4050070                                                                            |                                                                            ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggcgt                                                                        >  1:748025/1‑78 (MQ=255)      tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGttttt                                                                   <  1:747979/80‑1 (MQ=255)      tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGtttt                                                                    <  1:747980/79‑1 (MQ=255)        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTAt                                                                >  1:748019/1‑81 (MQ=255)         ttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTAtt                                                               >  1:748034/1‑81 (MQ=255)          tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTg                                                              >  2:748017/1‑81 (MQ=255)           cAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTgg                                                             >  2:748034/1‑81 (MQ=255)           cAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTgg                                                             >  2:748025/1‑81 (MQ=255)           cAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTgg                                                             >  1:748042/1‑81 (MQ=255)              aCCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCg                                                          <  1:747981/81‑1 (MQ=255)               ccTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCga                                                         <  2:747981/81‑1 (MQ=255)                   ttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCgaaag                                                      >  1:748037/1‑80 (MQ=255)                   ttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCgaaag                                                      >  2:748037/1‑80 (MQ=255)                     aCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCgaaagaa                                                    >  1:748032/1‑80 (MQ=255)                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCgaaagaaa                                                   <  1:747983/80‑1 (MQ=255)                           cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                             >  2:748053/1‑81 (MQ=255)                            gTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                            >  2:748030/1‑81 (MQ=255)                               aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                         <  2:747983/81‑1 (MQ=255)                                     tttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                   <  2:747984/81‑1 (MQ=255)                                      ttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGa                                  >  2:748044/1‑81 (MQ=255)                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGAt                                 <  1:747986/80‑1 (MQ=255)                                         gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                               <  1:747987/81‑1 (MQ=255)                                          cTGAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                              >  1:748055/1‑81 (MQ=255)                                            gAATATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                             >  2:748043/1‑80 (MQ=255)                                             aaTATTTTAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                            <  1:747988/80‑1 (MQ=255)                                                  tttAGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                       >  1:748041/1‑80 (MQ=255)                                                     aGCCGCCCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAt                    >  2:748038/1‑80 (MQ=255)                                                     aGCCGCCCCAGCCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAt                    >  2:748039/1‑80 (MQ=255)                                                      gccgccCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAGAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAt                    >  2:748046/1‑79 (MQ=255)                                                       ccgccCCAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                 <  1:747989/81‑1 (MQ=255)                                                          ccccAGTCAGTAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGc               >  2:748052/1‑80 (MQ=255)                                                              agtcagtAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg          <  2:747989/81‑1 (MQ=255)                                                                tcagtAATGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt        <  1:747991/81‑1 (MQ=255)                                                                     aaTGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt   >  1:748048/1‑81 (MQ=255)                                                                     aaTGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt   >  2:748058/1‑81 (MQ=255)                                                                      aTGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt  <  2:747991/81‑1 (MQ=255)                                                                      aTGACT‑GGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt  >  2:748069/1‑81 (MQ=255)                                                                            |                                                                            TTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATT  >  NC_000913/4049920‑4050070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑