Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,315,705	T→C	K22K (AAA→AAG)	yciF ←	putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,315,705	0	T	C	98.7%	269.8 / ‑4.4	75	K22K (AAA→AAG)	yciF	putative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (27/47);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCC  >  NC_000913/1315608‑1315799                                                                                                  |                                                                                               tCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGctt                                                                                              >  1:2009752/1‑100 (MQ=255)   tGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGctttt                                                                                            >  2:661000/1‑100 (MQ=255)      ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGCTAATTGCTTTTCTg                                                                                         >  1:1692710/1‑100 (MQ=255)       cATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgc                                                                                        >  2:1003164/1‑100 (MQ=255)         tGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcgc                                                                                      >  1:1771999/1‑100 (MQ=255)         tGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcgc                                                                                      >  1:1205890/1‑100 (MQ=255)          gAGTTTCCTCGAGGTGCGCATGAAAAGCCCGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTCTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCt                                                                                     <  2:1230280/100‑1 (MQ=255)               tCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAgg                                                                                <  2:980529/100‑1 (MQ=255)                ccTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGt                                                                               >  1:1936621/1‑100 (MQ=255)                 cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTa                                                                              <  1:1287018/100‑1 (MQ=255)                 cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTa                                                                              >  1:528651/1‑100 (MQ=255)                  tCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTAt                                                                             >  1:630403/1‑100 (MQ=255)                   tccgatctgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATc                                                                            <  1:1672706/92‑1 (MQ=255)                   cGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATc                                                                            >  2:1840232/1‑100 (MQ=255)                      ggTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGa                                                                         <  1:1091436/100‑1 (MQ=255)                       gTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGaa                                                                        <  1:661002/100‑1 (MQ=255)                       gTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGaa                                                                        <  1:2018752/100‑1 (MQ=255)                        tGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGaaa                                                                       <  1:1115074/100‑1 (MQ=255)                         gcgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAg                                                                      >  2:2005285/1‑100 (MQ=255)                           gcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATcagatcgga                                                                    >  2:1672693/1‑92 (MQ=255)                             aTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                  >  1:210112/1‑100 (MQ=255)                             aTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                  <  2:210111/100‑1 (MQ=255)                               gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                <  2:1031588/100‑1 (MQ=255)                                aaaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGa                                                               >  1:57401/1‑100 (MQ=255)                                 aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                              <  2:2009734/100‑1 (MQ=255)                                 aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                              <  1:378433/100‑1 (MQ=255)                                    gCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataa                                                           <  1:1422085/100‑1 (MQ=255)                                      cTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaat                                                         >  1:1095787/1‑100 (MQ=255)                                      cTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaat                                                         <  1:1453576/100‑1 (MQ=255)                                      cTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaat                                                         <  1:955386/100‑1 (MQ=255)                                       tGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                        >  1:402018/1‑100 (MQ=255)                                       tGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                        <  2:1058868/100‑1 (MQ=255)                                        gACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATAc                                                       <  2:1575807/100‑1 (MQ=255)                                         aCTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACa                                                      <  1:589415/100‑1 (MQ=255)                                           ttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATc                                                    <  1:1793857/100‑1 (MQ=255)                                           ttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATc                                                    <  2:1225266/100‑1 (MQ=255)                                              aTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTc                                                 <  2:970674/100‑1 (MQ=255)                                              aTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTc                                                 <  2:970610/100‑1 (MQ=255)                                               tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCa                                                <  1:1840248/100‑1 (MQ=255)                                                  ttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATg                                             <  2:891923/100‑1 (MQ=255)                                                  ttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATg                                             <  1:87897/100‑1 (MQ=255)                                                    cATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGt                                           >  2:986606/1‑100 (MQ=255)                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                          >  1:1366638/1‑100 (MQ=255)                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                          >  2:2010395/1‑100 (MQ=255)                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                          >  2:1288176/1‑100 (MQ=255)                                                     aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                          <  2:817629/100‑1 (MQ=255)                                                      tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCt                                         <  1:1726848/100‑1 (MQ=255)                                                       ttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCtt                                        <  2:662357/100‑1 (MQ=255)                                                         gATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttca                                      <  2:1083315/100‑1 (MQ=255)                                                         gATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttca                                      >  1:1422819/1‑100 (MQ=255)                                                             ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatatt                                  <  2:201608/100‑1 (MQ=255)                                                             ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatatt                                  <  1:1395324/100‑1 (MQ=255)                                                              tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                 <  2:1771985/100‑1 (MQ=255)                                                               gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattca                                <  2:528650/100‑1 (MQ=255)                                                                 tctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcata                              <  1:632575/100‑1 (MQ=255)                                                                  ctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTtcatattcatat                             >  1:744156/1‑100 (MQ=255)                                                                  ctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTtcatattcatat                             >  2:858516/1‑100 (MQ=255)                                                                  ctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTtcatattcatat                             <  2:1211798/100‑1 (MQ=255)                                                                     gCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                          >  2:1463277/1‑100 (MQ=255)                                                                     gCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                          >  1:796287/1‑100 (MQ=255)                                                                     gCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                          >  1:1648080/1‑100 (MQ=255)                                                                       gAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                        >  1:1480893/1‑100 (MQ=255)                                                                       gAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                        >  2:1428514/1‑100 (MQ=255)                                                                       gAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                        <  2:1205884/100‑1 (MQ=255)                                                                       gAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                        <  1:916750/100‑1 (MQ=255)                                                                         gTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtctc                      <  2:1646069/100‑1 (MQ=255)                                                                          ttttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTcc                     <  1:2005303/100‑1 (MQ=255)                                                                                 aGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGaaa              <  1:1003169/100‑1 (MQ=255)                                                                                   tGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATc            <  1:1053336/100‑1 (MQ=255)                                                                                       cGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTg        <  1:43261/100‑1 (MQ=255)                                                                                        gggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGc       <  1:877513/100‑1 (MQ=255)                                                                                         ggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGACTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg      <  2:1095782/100‑1 (MQ=255)                                                                                           ttAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAg    <  2:1615578/100‑1 (MQ=255)                                                                                           ttAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAg    <  2:1366629/100‑1 (MQ=255)                                                                                            tAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc   <  1:1602636/100‑1 (MQ=255)                                                                                             aaTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc  >  2:269885/1‑100 (MQ=255)                                                                                             aaTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc  <  2:796284/100‑1 (MQ=255)                                                                                                  |                                                                                               TCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCC  >  NC_000913/1315608‑1315799

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