Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,188,572	C→T	R1075C (CGT→TGT)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,188,572	0	C	T	98.4%	215.4 / NA	63	R1075C (CGT→TGT)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/32);  minor base G (0/1);  total (30/33)
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

AAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCA  >  NC_000913/4188442‑4188695                                                                                                                                   |                                                                                                                            aaGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCa                                                                                                                     <  1:234677/139‑1 (MQ=255)   gCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCatt                                                                                                                   >  2:821214/1‑137 (MQ=255)         ttGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAAt                                                                                                             >  2:263971/1‑139 (MQ=255)         ttGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAAt                                                                                                             >  2:808142/1‑139 (MQ=255)                   tACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCt                                                                                                   >  1:939689/1‑139 (MQ=255)                   tACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCt                                                                                                   <  2:939689/139‑1 (MQ=255)                        aCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAggg                                                                                              <  1:263971/139‑1 (MQ=255)                             aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAAt                                                                                                             >  2:31789/1‑119 (MQ=255)                             aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAAt                                                                                                             <  1:31789/119‑1 (MQ=255)                             aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTg                                                                                                        <  1:149155/124‑1 (MQ=255)                             aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTg                                                                                                        >  2:149155/1‑124 (MQ=255)                             aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacg                                                                                         >  2:198428/1‑139 (MQ=255)                             aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacg                                                                                         >  2:359741/1‑139 (MQ=255)                               cGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacgac                                                                                       <  1:70906/139‑1 (MQ=255)                                  cGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGtt                                                                                    >  1:365684/1‑139 (MQ=255)                                        gACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATc                                                                              <  2:365684/139‑1 (MQ=255)                                         aCCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATcc                                                                             <  1:883538/139‑1 (MQ=255)                                         aCCATTACGCGCCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATcc                                                                             <  1:808142/139‑1 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGtt                                                                                                         <  2:679639/106‑1 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGtt                                                                                                         >  1:679639/1‑106 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATc                                                                              <  2:217398/133‑1 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATc                                                                              >  1:217398/1‑133 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGt                                                                        <  1:359741/139‑1 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGt                                                                        <  2:942515/139‑1 (MQ=255)                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGt                                                                        >  1:1143522/1‑139 (MQ=255)                                                   gTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGTTCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGa                                                                   >  2:795603/1‑139 (MQ=255)                                                        aCCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGc                                                              >  2:602320/1‑139 (MQ=255)                                                        aCCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGc                                                              <  2:1143522/139‑1 (MQ=255)                                                           gacgaACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAg                                                           <  1:1156075/139‑1 (MQ=255)                                                               aaCTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCa                                                       >  1:352415/1‑139 (MQ=255)                                                                          ctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGc                                            <  1:989900/139‑1 (MQ=255)                                                                           tgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGcc                                           >  1:134081/1‑139 (MQ=255)                                                                                       tggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCt                                              >  2:572519/1‑124 (MQ=255)                                                                                       tggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCt                                              <  1:572519/124‑1 (MQ=255)                                                                                       tggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAt                               >  1:269969/1‑139 (MQ=255)                                                                                       tggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAt                               >  2:506471/1‑139 (MQ=255)                                                                                            cTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTc                          <  1:1214547/139‑1 (MQ=255)                                                                                            cTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTc                          <  1:564117/139‑1 (MQ=255)                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCt                                                                                                   >  1:713343/1‑58 (MQ=255)                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCt                                                                                                   <  2:713343/58‑1 (MQ=255)                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATcc                                                                             >  2:1189325/1‑80 (MQ=255)                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATcc                                                                             <  1:1189325/80‑1 (MQ=255)                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGt                                                                        >  2:637630/1‑85 (MQ=255)                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGt                                                                        <  1:637630/85‑1 (MQ=255)                                                                                                      cAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGa                <  2:134081/139‑1 (MQ=255)                                                                                                      cAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGa                >  1:8807/1‑139 (MQ=255)                                                                                                      cAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGa                >  2:186406/1‑139 (MQ=255)                                                                                                        gAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATg              <  1:795603/139‑1 (MQ=255)                                                                                                        gAACGTACCGCAGGTGGTAAAGATCTGGGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGTCGTTCTGATCCCACGTACCGATATGCCAGGGCTGTAATTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATg              <  2:156354/139‑1 (MQ=255)                                                                                                           cGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCg                                          <  1:683797/108‑1 (MQ=255)                                                                                                           cGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCg                                          >  2:683797/1‑108 (MQ=255)                                                                                                            gTACCGCAGGTGGTAATGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGt          <  1:841602/139‑1 (MQ=255)                                                                                                              aCCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATg                                                                                                     >  1:677980/1‑46 (MQ=255)                                                                                                              aCCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATg                                                                                                     <  2:677980/46‑1 (MQ=255)                                                                                                              aCCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacgac                                                                                       >  2:1075424/1‑60 (MQ=255)                                                                                                              aCCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacgac                                                                                       <  1:1075424/60‑1 (MQ=255)                                                                                                                  cAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGAt    <  2:352415/139‑1 (MQ=255)                                                                                                                   aGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATc   >  1:317832/1‑139 (MQ=255)                                                                                                                   aGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATc   >  2:1221669/1‑139 (MQ=255)                                                                                                                    ggtggtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCa  <  2:317832/139‑1 (MQ=255)                                                                                                                     gtggtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAgc                                                          >  1:949874/1‑82 (MQ=255)                                                                                                                     gtggtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAgc                                                          <  2:949874/82‑1 (MQ=255)                                                                                                                        gtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAtt                              >  1:209854/1‑107 (MQ=255)                                                                                                                        gtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAtt                              <  2:209854/107‑1 (MQ=255)                                                                                                                                   |                                                                                                                            AAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCA  >  NC_000913/4188442‑4188695

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