Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,421,120	Δ1 bp	noncoding (60/121 nt)	rrfF ←	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,421,119	0	G	.	94.1%	125.7 / 4.6	51	noncoding (61/121 nt)	rrfF	5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (26/22);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421028‑3421207                                                                                            |                                                                                         aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  1:652979/1‑92 (MQ=255) aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  2:2254355/1‑92 (MQ=255)  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:2975630/1‑91 (MQ=255)  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  1:3419141/1‑91 (MQ=255)   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  1:2094162/1‑93 (MQ=255)   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  1:2433565/1‑93 (MQ=255)   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       <  2:1615664/93‑1 (MQ=37)     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCa                                                                                    >  1:872931/1‑93 (MQ=255)     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCa                                                                                    >  1:4312309/1‑93 (MQ=255)      tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCAc                                                                                   <  1:2450653/93‑1 (MQ=16)            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             <  2:905160/93‑1 (MQ=9)             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            <  1:1948681/93‑1 (MQ=9)              aGCCTTTTGTTTTATTTTATGTCTGGCAGTTCCCTACTCTCGCAT‑‑GGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  2:1256796/93‑1 (MQ=11)              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGGCCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                           >  2:1027996/1‑93 (MQ=255)              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACGCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                           <  1:3346032/93‑1 (MQ=255)              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                           >  2:4110919/1‑93 (MQ=16)               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                          <  2:4381072/93‑1 (MQ=9)                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACGCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  2:2534805/93‑1 (MQ=9)                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  2:4288688/93‑1 (MQ=9)                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGTGTAGACCCCACACTACCATCGGCGCTACGGC‑GTGTCACTTCTGAGTTCg                                                                        >  2:4209244/1‑93 (MQ=11)                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACGCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        <  2:1580297/93‑1 (MQ=9)                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACGCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:333757/1‑93 (MQ=9)                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:4477545/1‑93 (MQ=12)                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:4505142/1‑93 (MQ=12)                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:4032718/1‑93 (MQ=12)                                                                    ccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     >  2:4559710/1‑93 (MQ=9)                                                                    ccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     >  2:1067291/1‑93 (MQ=9)                                                                     cacgCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGAGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:2242071/93‑1 (MQ=11)                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    >  2:1485850/1‑93 (MQ=12)                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:907389/93‑1 (MQ=16)                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    >  1:4650190/1‑93 (MQ=12)                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:3176845/93‑1 (MQ=16)                                                                      acgcTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   >  1:2703431/4‑93 (MQ=255)                                                                      acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   >  1:755066/1‑93 (MQ=16)                                                                      acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   <  1:4646584/93‑1 (MQ=16)                                                                      acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   <  2:2637772/93‑1 (MQ=16)                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  <  1:1729213/93‑1 (MQ=19)                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  >  2:857033/1‑93 (MQ=19)                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  <  2:1691945/93‑1 (MQ=19)                                                                        acTCCCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:1735814/93‑1 (MQ=255)                                                                         cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                <  1:756426/93‑1 (MQ=255)                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               >  1:796943/1‑93 (MQ=255)                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  1:2717584/93‑1 (MQ=255)                                                                           aCCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt              >  2:3279783/1‑93 (MQ=255)                                                                             cgccGGCGCTGCGG‑GTTTTCTTTTCTGAGTTCGGCATGGGGTCAGGGGGGACCCCCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            <  2:3910756/90‑1 (MQ=17)                                                                              aTCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGGCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           >  1:4024689/1‑93 (MQ=255)                                                                               tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      >  1:2428708/1‑81 (MQ=0)                                                                               tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg          >  2:1494580/1‑93 (MQ=255)                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  1:358510/93‑1 (MQ=255)                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  1:3493670/93‑1 (MQ=255)                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  2:1387760/93‑1 (MQ=255)                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       >  1:886444/1‑93 (MQ=255)                                                                                     cTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACttt    >  2:2463271/1‑93 (MQ=255)                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   >  2:3765602/1‑93 (MQ=255)                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   >  2:4742289/1‑93 (MQ=255)                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  1:476158/1‑93 (MQ=255)                                                                                            |                                                                                         AGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421028‑3421207

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