Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,734	Δ1 bp	noncoding (409/2906 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,734	0	A	.	98.3%	103.9 / ‑1.2	60	noncoding (409/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (29/29);  minor base G (0/1);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCG  >  NC_000913/3941644‑3941814                                                                                           |                                                                                 gATACAGGGTGACAGCCGCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                >  2:1008889/1‑93 (MQ=255) gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                >  2:1419095/1‑93 (MQ=25)   tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                                             >  1:2337120/1‑93 (MQ=12)     cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa                                                                           <  1:833403/93‑1 (MQ=9)      aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                                                          >  2:259374/1‑93 (MQ=16)        ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                                                        >  1:1692389/1‑93 (MQ=9)        ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                                                        <  1:342070/93‑1 (MQ=9)        ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                                                        >  1:1816060/1‑93 (MQ=16)         gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct                                                                       <  2:4631702/93‑1 (MQ=9)           gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctcc                                                                     >  1:216530/1‑93 (MQ=16)            aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa                                                                    <  2:1915162/93‑1 (MQ=9)            aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa                                                                    >  2:4235845/1‑93 (MQ=16)             cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                   <  2:3285852/93‑1 (MQ=9)             cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                   <  1:68039/93‑1 (MQ=9)             cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                   <  1:2918840/93‑1 (MQ=9)             cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCAGAATATGGG‑GGGACCATCCTCCaa                                                                   <  1:2918839/93‑1 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                                                  <  1:752246/93‑1 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                                                  <  2:1316644/93‑1 (MQ=9)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                                                  >  1:3951287/1‑93 (MQ=16)              aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                                                  <  2:2002770/93‑1 (MQ=9)               gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCACgg                                                                 >  1:1017678/1‑93 (MQ=255)               gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                                                 >  2:3988286/1‑93 (MQ=16)                ccccGTACACAAAAATGCACATATTGTGAGCTCGATGCGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCACCCTCCAAGGc                                                                >  2:2490370/1‑93 (MQ=11)                ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                                                <  1:3420445/93‑1 (MQ=9)                 cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                                               >  1:4641993/1‑93 (MQ=16)                 cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCt                                                               >  1:1062606/1‑93 (MQ=12)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                              <  1:3828512/93‑1 (MQ=9)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                              >  1:1812181/1‑93 (MQ=16)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                              >  2:3844243/1‑93 (MQ=16)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                              >  1:1842220/1‑93 (MQ=16)                  ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                              >  1:3063362/1‑93 (MQ=16)                   cGTACACAAACATCCACCTATTCTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGGCCATCCTCCAAGGCTaa                                                             <  2:4263767/93‑1 (MQ=11)                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                            >  2:328265/1‑93 (MQ=9)                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                            <  2:4182017/93‑1 (MQ=9)                    gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                            >  1:4238343/1‑93 (MQ=16)                     tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                                                           >  1:1461690/1‑93 (MQ=16)                     tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                                                           <  2:2948037/93‑1 (MQ=9)                     tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                                                           >  2:2147577/1‑93 (MQ=16)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                          <  2:2377144/93‑1 (MQ=9)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                          >  2:370100/1‑93 (MQ=16)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                          <  2:992233/93‑1 (MQ=9)                      acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                          <  1:1947753/93‑1 (MQ=9)                       ctcaAAAATGCACATATTGTGAGCTCGATGAGTAGTGCGGGACACGTGGTATCCTGTCTAAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                                                         <  2:3153557/91‑1 (MQ=11)                        acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                                                        >  1:1103197/1‑93 (MQ=16)                         caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAGACTc                                                       >  2:1361875/1‑93 (MQ=255)                         caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGCGGGACCATCCTCCAAGGCTAAATACTc                                                       >  1:149662/1‑93 (MQ=9)                          aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                                                      >  1:1057517/1‑93 (MQ=16)                          aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                                                      <  1:3426493/93‑1 (MQ=9)                           aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                                                     <  1:3183478/93‑1 (MQ=9)                           aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                                                     <  2:3562189/93‑1 (MQ=9)                           aaaaTGCACATATTGGGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                                                     <  2:3592890/93‑1 (MQ=9)                            aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                                                    >  1:365286/1‑93 (MQ=9)                            aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                                                    >  1:498244/1‑93 (MQ=16)                            aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctg                                                    <  2:4492547/93‑1 (MQ=9)                              aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                                                  >  1:4315706/1‑93 (MQ=9)                              aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                                                  <  2:1813297/93‑1 (MQ=9)                               tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGGCCATCCTCCAAGGCTAAATACTCctgact                                                 >  1:2886788/1‑93 (MQ=255)                               tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                                                 >  2:1286667/1‑93 (MQ=16)                                gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactg                                                >  1:65703/1‑93 (MQ=16)                                 cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga                                               >  2:3397992/1‑93 (MQ=9)                                  acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                                              <  2:1871211/93‑1 (MQ=9)                                  acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                                              <  1:2059524/93‑1 (MQ=9)                                  acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                                              >  2:57809/1‑93 (MQ=9)                                   caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                             <  2:628544/93‑1 (MQ=0)                                                                               tCTGAATATGGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg  <  2:1572547/93‑1 (MQ=16)                                                                                           |                                                                                 GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCG  >  NC_000913/3941644‑3941814

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