Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,477,434 T→C intergenic (‑15/‑206) fkpA ← / → slyX FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)/phi X174 lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,4340TC98.7% 268.4 / ‑4.4 76intergenic (‑15/‑206)fkpA/slyXFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)/phi X174 lysis protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (44/31);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCTATTATCATCTA  >  NC_000913/3477339‑3477528
                                                                                               |                                                                                              
ttCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCggg                                                                                            >  2:1517810/1‑100 (MQ=255)
 tCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCgggg                                                                                           >  2:936961/1‑100 (MQ=255)
 tCAGCAGCAAAAATGGTTGGTGCATGCAGGGCCACCGGCATTTTGGGCGCCAGCAGCGTTACTTTAAACAATGATTTCATCCATATCTCCCGGACCgggg                                                                                           >  2:695686/1‑100 (MQ=255)
  cagcagCTAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGc                                                                                          <  2:849826/100‑1 (MQ=255)
   agcagcAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCa                                                                                         <  1:726309/100‑1 (MQ=255)
   agcagcAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCa                                                                                         >  1:530071/1‑100 (MQ=255)
     cagcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCAtc                                                                                       <  1:537088/100‑1 (MQ=255)
     cagcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCAtc                                                                                       <  1:666093/100‑1 (MQ=255)
     cagcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCAtc                                                                                       <  2:1106895/100‑1 (MQ=255)
       gcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCAtctc                                                                                     >  1:280757/1‑100 (MQ=255)
       gcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCAtctc                                                                                     <  2:1339796/100‑1 (MQ=255)
        caAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCa                                                                                    >  2:855707/1‑100 (MQ=255)
         aaaaGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCAc                                                                                   <  1:167710/100‑1 (MQ=255)
         aaaaGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCAc                                                                                   >  2:1176963/1‑100 (MQ=255)
              tGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAg                                                                              <  1:1144339/100‑1 (MQ=255)
              tGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCAGGGCATCTCACCCCAg                                                                              <  1:1388234/100‑1 (MQ=255)
               gATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAgg                                                                             >  1:834962/1‑100 (MQ=255)
               gATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAgg                                                                             <  1:166053/100‑1 (MQ=255)
               gATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAgg                                                                             >  1:186305/1‑100 (MQ=255)
                aTTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAggg                                                                            >  2:917289/1‑100 (MQ=255)
                aTTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAggg                                                                            <  1:834429/100‑1 (MQ=255)
                  tGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGtt                                                                          <  2:356635/100‑1 (MQ=255)
                   ggTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCGGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTa                                                                         >  1:1112744/1‑100 (MQ=255)
                   ggTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTa                                                                         >  1:168367/1‑100 (MQ=255)
                      gcatgcaGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACt                                                                      >  1:584644/1‑100 (MQ=255)
                      gcatgcaGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACt                                                                      >  2:1327909/1‑100 (MQ=255)
                      gcatgcaGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGAGTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACt                                                                      >  1:584186/1‑100 (MQ=255)
                      gcatgcaGGGCAACGGACATTGTGGTCTCCAGCGGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGtt                                                                          <  2:1063960/96‑1 (MQ=255)
                       catgaaGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACta                                                                     <  1:625606/100‑1 (MQ=255)
                          gcaGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTAGTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACtatta                                                                  >  1:1118036/1‑100 (MQ=255)
                          gcaGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACtatta                                                                  >  2:139171/1‑100 (MQ=255)
                            aGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAAGTATTATc                                                                >  1:1205823/1‑100 (MQ=255)
                                    ggCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGt                                                        >  1:917278/1‑100 (MQ=255)
                                     gCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTa                                                       >  1:1107291/1‑100 (MQ=255)
                                     gCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTa                                                       >  1:572247/1‑100 (MQ=255)
                                      ccATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTAc                                                      <  2:378063/100‑1 (MQ=255)
                                      ccATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTAc                                                      >  1:201585/1‑100 (MQ=255)
                                      ccATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTAc                                                      >  1:1182498/1‑100 (MQ=255)
                                       cATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACt                                                     <  2:953598/100‑1 (MQ=255)
                                        aTTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACta                                                    <  1:855712/100‑1 (MQ=255)
                                          tgtgGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACtata                                                  <  2:834957/100‑1 (MQ=255)
                                          tgtgGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACtata                                                  <  1:695688/100‑1 (MQ=255)
                                           gtgGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATaa                                                 <  1:917294/100‑1 (MQ=255)
                                           gtgGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATaa                                                 >  2:140731/1‑100 (MQ=255)
                                              gTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGc                                              >  1:792398/1‑100 (MQ=255)
                                                  ccAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGTGTACTATAAAGCGTTg                                          <  2:186205/100‑1 (MQ=255)
                                                  ccAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTg                                          <  1:1517819/100‑1 (MQ=255)
                                                  ccAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTg                                          >  2:1059208/1‑100 (MQ=255)
                                                  ccAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAAATATTATCAGAAGGGTACTATAAAGCGTTg                                          >  2:798297/1‑100 (MQ=255)
                                                   cagcagCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGt                                         <  2:319186/100‑1 (MQ=255)
                                                   cagcagCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGt                                         >  1:48914/1‑100 (MQ=255)
                                                     gcagcGTTACTTTAATCAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCg                                       <  1:1555191/100‑1 (MQ=255)
                                                      cagcGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGa                                      <  1:604619/100‑1 (MQ=255)
                                                        gcGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAAc                                    >  2:1039242/1‑100 (MQ=255)
                                                        gcGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAAc                                    <  2:632981/100‑1 (MQ=255)
                                                          gTTACTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACaa                                  >  1:333249/1‑100 (MQ=255)
                                                              cTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCt                              >  1:54183/1‑100 (MQ=255)
                                                              cTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCt                              >  1:873555/1‑100 (MQ=255)
                                                              cTTTAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCt                              >  2:293424/1‑100 (MQ=255)
                                                                ttAAACAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCtac                            <  2:1360031/100‑1 (MQ=255)
                                                                 tAAACAGCGATATCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTCTTATCAGAAGGGTACTATAAAGCGGTGTCGAACAAATCtaca                           >  2:831783/1‑100 (MQ=255)
                                                                    aCAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCtacatac                        >  1:399900/1‑100 (MQ=255)
                                                                    aCAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCtacatac                        <  2:1154077/100‑1 (MQ=255)
                                                                     cAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCtacataca                       >  2:35794/1‑100 (MQ=255)
                                                                     cAGTGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCtacataca                       >  1:638078/1‑100 (MQ=255)
                                                                        tGATTTCATCCATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGac                    >  1:1259709/1‑100 (MQ=255)
                                                                              catccatcTCTGCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAGGCGTTGTCGAACAAATCTACATACAGACACGccc              >  2:831511/1‑100 (MQ=255)
                                                                                 ccATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCta            >  2:1341179/1‑99 (MQ=255)
                                                                                  cATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCtatt          >  1:734820/1‑100 (MQ=255)
                                                                                  cATATCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCtatt          >  2:8099/1‑100 (MQ=255)
                                                                                   atatCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCtatta         <  1:537183/100‑1 (MQ=255)
                                                                                    tatCTCCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCtattat        <  2:1305971/100‑1 (MQ=255)
                                                                                     atCTCCAGGACCGGGGCATCTCACCCCAGTGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTTGAACAAATCTACATACAGACACGCCCTATTatc       <  2:948366/100‑1 (MQ=255)
                                                                                       ctcCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCTATTatcat     >  2:1400479/1‑100 (MQ=255)
                                                                                        tcCAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCTATTatcatc    >  1:366959/1‑100 (MQ=255)
                                                                                          cAGGACCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCTATTATCATCTa  <  1:1092039/100‑1 (MQ=255)
                                                                                               |                                                                                              
TTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAACAAATCTACATACAGACACGCCCTATTATCATCTA  >  NC_000913/3477339‑3477528

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