Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,879	1	.	C	98.1%	186.0 / ‑3.4	52	coding (671/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (21/30);  total (21/31)
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.97e-01

CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGT  >  NC_000913/3815741‑3816011                                                                                                                                              |                                                                                                                                      cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgat‑t‑c                                                                                                                                      <  1:848019/140‑1 (MQ=255)   aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑ccc                                                                                                                                    >  1:993901/1‑140 (MQ=255)   aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑ccc                                                                                                                                    >  1:569523/1‑140 (MQ=255)      aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑ccccct                                                                                                                                 <  1:1196762/140‑3 (MQ=255)          ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCggg                                                                                                                              >  2:406377/1‑140 (MQ=255)           gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGc                                                                                                                             <  1:826135/140‑1 (MQ=255)              cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCag                                                                                                                          <  2:254562/140‑1 (MQ=255)                gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCagag                                                                                                                        <  1:276580/140‑1 (MQ=255)                  tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGcc                                                                                                                      <  2:279963/140‑1 (MQ=255)                     ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCcaac                                                                                                                   <  2:1209069/140‑1 (MQ=255)                      gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCcaaca                                                                                                                  <  2:578923/140‑1 (MQ=255)                         tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGACCAGAGCCAACAAGa                                                                                                               >  2:998812/1‑140 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGcc                                                                                                                            >  2:229675/1‑125 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGcc                                                                                                                            <  1:229675/125‑1 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:396687/1‑140 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:531475/1‑140 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:834673/1‑140 (MQ=255)                             gCCTTCGCTCCTCATCTTACTTTCCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAg                                                                                                           <  2:1363008/140‑1 (MQ=255)                                ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        >  2:432881/1‑140 (MQ=255)                                ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  2:1070048/140‑1 (MQ=255)                                ctCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTCCGCCACCGCCGCCTTCTGCGTCGCTACAATGGATTCGCTCC‑CCCCTCGGGCCAGAGCGAACAAGATGAGTAg                                                                                                        <  2:901746/139‑1 (MQ=255)                                  cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                                                                                      <  1:1246019/140‑1 (MQ=255)                                                   ttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGACCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                                                     <  1:998812/140‑1 (MQ=255)                                                    tCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                                                     >  1:1098337/1‑140 (MQ=255)                                                    tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                    <  1:396687/140‑1 (MQ=255)                                                     cTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                   >  1:430657/1‑140 (MQ=255)                                                        cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                >  2:889907/1‑140 (MQ=255)                                                          gAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                              <  1:432881/140‑1 (MQ=255)                                                            c‑AAAAAAAAGGCGACTCATCAGTCG‑CTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           <  2:1164286/140‑1 (MQ=255)                                                              aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           >  2:832206/1‑140 (MQ=255)                                                              aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           <  2:800902/140‑1 (MQ=255)                                                              aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           <  2:654688/140‑1 (MQ=255)                                                              aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                           >  1:172776/1‑140 (MQ=255)                                                                 aaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCg                                                                        >  2:201478/1‑140 (MQ=255)                                                                      ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                   <  1:889907/140‑1 (MQ=255)                                                                      ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                   <  1:406377/140‑1 (MQ=255)                                                                        cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑Tg                                                                <  1:832206/140‑1 (MQ=255)                                                                        cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑Tg                                                                <  1:531475/140‑1 (MQ=255)                                                                        cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑Tg                                                                >  2:349640/1‑140 (MQ=255)                                                                        cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑Tg                                                                <  1:914468/140‑1 (MQ=255)                                                                          aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCa                                                              <  2:430657/140‑1 (MQ=255)                                                                                     gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAatgatg                                                   >  2:60607/1‑140 (MQ=255)                                                                                     gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAatgatg                                                   >  1:814504/1‑140 (MQ=255)                                                                                     gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑GGCAATTCAatggtg                                                   >  2:837740/1‑140 (MQ=255)                                                                                      ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATgc                                                  >  2:639306/1‑140 (MQ=255)                                                                                              aTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCt                                          <  2:1039173/140‑1 (MQ=255)                                                                                                 aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCg                                       >  1:224725/1‑140 (MQ=255)                                                                                                     tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCATGCACTTCAATGATGCGCCCGTCTTCgg                                      <  1:24695/138‑1 (MQ=255)                                                                                                     tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCATGCACTTCAATGATGCGCCCGTCTTCgg                                      >  2:24695/1‑138 (MQ=255)                                                                                                      gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGtcat                                  <  2:993901/140‑1 (MQ=255)                                                                                                            gccgccTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAc                            >  1:1146536/1‑140 (MQ=255)                                                                                                              cgccTTCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                          <  2:1098337/140‑1 (MQ=255)                                                                                                                   tCTGCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                     <  1:834673/140‑1 (MQ=255)                                                                                                                      gCGTCGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                  <  1:349640/140‑1 (MQ=255)                                                                                                                          cGCTACAATGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct              <  2:172776/140‑1 (MQ=255)                                                                                                                                  tGGATTCGAT‑TCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa      <  1:639306/140‑1 (MQ=255)                                                                                                                                      ttcgat‑tcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt  >  1:328239/1‑140 (MQ=255)                                                                                                                                              |                                                                                                                                      CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC‑TGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGT  >  NC_000913/3815741‑3816011

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