Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,652 G→C 8.8% A578A (GCG→GCC rsxC → SoxR [2Fe‑2S] reducing system protein RsxC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,6520GC 8.8% 154.2 / 3.0 54A578A (GCG→GCCrsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base G (16/33);  new base C (2/3);  total (18/36)
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

GAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGC  >  NC_000913/1708503‑1708775
                                                                                                                                                      |                                                                                                                           
gAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                                                                                             <  2:1326332/150‑1 (MQ=255)
    cccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACTCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                                         <  1:678061/150‑1 (MQ=255)
    cccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                                                         >  1:607028/1‑150 (MQ=255)
     ccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAACAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAcc                                                                                                                        <  1:395250/150‑1 (MQ=255)
       gTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAg                                                                                                                      >  1:863896/1‑150 (MQ=255)
         aaaGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaa                                                                                                                    >  1:212483/1‑150 (MQ=255)
          aaGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaac                                                                                                                   >  1:10774/1‑150 (MQ=255)
           aGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAACCAGaaca                                                                                                                  >  1:1407350/1‑150 (MQ=255)
                        gCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcg                                                                                                     >  1:787007/1‑150 (MQ=255)
                            aaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCa                                                                                                 <  1:557868/150‑1 (MQ=255)
                                ggCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGACAAATGAGGACCCAGAACAACAGGTCGATACGCGCAAAgc                                                                                             >  2:1400118/1‑150 (MQ=25)
                                      aCTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATAGCGGAACAGAACACCACGGCGCTCCGCTCACCGcccccgg                                                                                       >  1:408643/1‑149 (MQ=25)
                                           agcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                  <  2:61883/150‑1 (MQ=255)
                                            gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc                                                                                 <  2:607028/150‑1 (MQ=255)
                                             caAACTTAGGACGCAGCAACCGTTGCCGATCCTCGTAAATCTGCCGTAGAAGCAGCTATTGCCCGTGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                <  2:343923/150‑1 (MQ=14)
                                              aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGACCACCGGGGCCCTCAGCGCACCGCAGCCGTCGCAGCCg                                                                               >  1:607563/1‑150 (MQ=21)
                                               aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGCTGCGCGCAAGCCCGCCGGCGAAGCCGc                                                                              >  1:701476/1‑150 (MQ=21)
                                                aCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                             <  2:10774/150‑1 (MQ=255)
                                                 ctcacgaccCAGCACCCGATGCTGATCCACGAAAACCTGCCGTCGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                            <  1:859694/147‑1 (MQ=17)
                                                 cTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTa                                                                            <  1:1529733/150‑1 (MQ=255)
                                                  tcccgacgCAGAAACAGAGGATGATCCTCGAAAATCTGCCGTTGACGCAGCTATTTACCGGGCCAAACCTCGCAATCTGGAACATCAACAGGATAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                                           <  2:1022511/147‑1 (MQ=2)
                                                   aaCGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCAGCGCCAAAGCGCGCAAGCTGGAACAGCAGCACGGTAATGCGGACCGCGAACATCCGGTCCCGCACCGCAAACCGGGCGGCGCAGCCGCTAta                                                                          >  1:832539/1‑149 (MQ=19)
                                                     cgacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                        <  1:922338/150‑1 (MQ=255)
                                                      gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                        <  1:286100/149‑1 (MQ=255)
                                                      gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                        >  2:286100/1‑149 (MQ=255)
                                                            gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                 <  1:533957/150‑1 (MQ=255)
                                                             caACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                                                <  2:787007/150‑1 (MQ=255)
                                                                ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAg                                                             <  2:1146645/150‑1 (MQ=255)
                                                                  gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcg                                                           <  2:1327078/150‑1 (MQ=255)
                                                                    tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                         <  1:141580/150‑1 (MQ=255)
                                                                      cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCTCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                       <  2:147906/150‑1 (MQ=255)
                                                                      cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCa                                                       <  1:666879/150‑1 (MQ=255)
                                                                       tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                                                     >  1:621466/1‑119 (MQ=255)
                                                                       tGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCg                                                                                     <  2:621466/119‑1 (MQ=255)
                                                                               gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                              <  1:1433910/150‑1 (MQ=255)
                                                                               gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                              <  2:212483/150‑1 (MQ=255)
                                                                                tAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa                                             <  2:1519846/150‑1 (MQ=255)
                                                                                 aaaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGACGTCGAAGCAGCGATTGCCCGTCCCAAAGCGCGCAAGCAGGACCac                                            >  1:758225/1‑149 (MQ=17)
                                                                                   aaCTGCCGTTGAAGCAGCTATCGCCCGCGACAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                          <  1:1531142/150‑1 (MQ=255)
                                                                                      tGCCGTTGAAGAAGCTATCGCCCGCGCCAAATCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                       <  1:1276053/150‑1 (MQ=255)
                                                                                          gTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                  <  2:947402/103‑1 (MQ=255)
                                                                                          gTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                  >  1:947402/1‑103 (MQ=255)
                                                                                            tGAAGCAGCTATCCCCCCCGCAAAATCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                 <  1:1414262/150‑1 (MQ=255)
                                                                                             gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                                         <  1:833952/109‑1 (MQ=255)
                                                                                             gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                                         >  2:833952/1‑109 (MQ=255)
                                                                                                gcagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                             <  1:218267/150‑1 (MQ=255)
                                                                                                         gCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                       <  2:278866/99‑1 (MQ=37)
                                                                                                         gCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                       >  1:278866/1‑99 (MQ=37)
                                                                                                          cccGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCagaag                   <  1:919967/150‑1 (MQ=255)
                                                                                                                ccACAGCGCGCAAGCTGGAACACCAACAGGCTAATGCGGAACCAGAACCCGAGGGCGATCCCCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg             <  2:816266/150‑1 (MQ=14)
                                                                                                                  aaatctcgcAAGCTGGAAGAGCAACAGGCTAATGAGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc           <  1:1016860/144‑1 (MQ=255)
                                                                                                                          cAAGCTTGAACATCAACATGCTAATGGGTATCCATAACATCAGGTCTATCGTCGTAAATCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg   <  2:758225/150‑1 (MQ=14)
                                                                                                                           aaccTGGAACATCTTCAGCATACTTCGGAACCAGCACCACAGGCCGATCCCCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCgcgc  <  1:1051659/147‑1 (MQ=14)
                                                                                                                              cTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                             <  2:16200/72‑1 (MQ=34)
                                                                                                                              cTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                             >  1:16200/1‑72 (MQ=34)
                                                                                                                                                      |                                                                                                                           
GAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGC  >  NC_000913/1708503‑1708775

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