Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,028,955 C→T 28.8% V20I (GTA→ATA)  ygfS ← putative 4Fe‑4S ferredoxin‑type oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,028,9550CT28.8% 104.4 / 49.0 73V20I (GTA→ATA) ygfSputative 4Fe‑4S ferredoxin‑type oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base C (22/30);  new base T (12/9);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

CCGGCGCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTAACATATCGCGCG  >  NC_000913/3028856‑3029052
                                                                                                   |                                                                                                 
ccGGCGCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATGGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTAc                                                                                                   >  2:1334480/1‑100 (MQ=255)
 cGGCGCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACt                                                                                                  <  2:1677163/100‑1 (MQ=255)
   gcgcACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTc                                                                                                >  2:1327888/1‑100 (MQ=255)
   gcgcACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTc                                                                                                <  1:1467613/100‑1 (MQ=255)
      cACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTcaca                                                                                             <  1:2049511/100‑1 (MQ=255)
      cACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTcaca                                                                                             >  1:143674/1‑100 (MQ=255)
      cACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTcaca                                                                                             <  1:1869057/100‑1 (MQ=255)
       aCTGATGCTGTCCAGTCGCGGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACaa                                                                                            <  1:555103/100‑1 (MQ=255)
         tGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGt                                                                                          <  1:1388058/100‑1 (MQ=255)
          gATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTa                                                                                         >  1:714695/1‑100 (MQ=255)
          gATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGAAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTa                                                                                         <  2:2206021/100‑1 (MQ=255)
               tGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCa                                                                                    >  1:1182735/1‑100 (MQ=255)
               tGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCa                                                                                    <  1:43305/100‑1 (MQ=255)
               tGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCa                                                                                    <  2:1602993/100‑1 (MQ=255)
                gTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAg                                                                                   >  2:1798813/1‑100 (MQ=255)
                gTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAg                                                                                   <  1:623868/100‑1 (MQ=255)
                gTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAg                                                                                   >  1:854962/1‑100 (MQ=255)
                  ccAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGc                                                                                  >  1:335708/1‑99 (MQ=255)
                   cAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCa                                                                                <  2:63370/100‑1 (MQ=255)
                    aGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCaa                                                                               >  2:1518232/1‑100 (MQ=255)
                     gTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAAt                                                                              >  1:1197361/1‑100 (MQ=255)
                       cGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATAc                                                                            >  1:1524493/1‑100 (MQ=255)
                       cGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATAc                                                                            >  1:1659166/1‑100 (MQ=255)
                        gCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACa                                                                           <  1:1502301/100‑1 (MQ=255)
                         cTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAg                                                                          <  2:1182731/100‑1 (MQ=255)
                              ccTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGc                                                                     >  2:935272/1‑100 (MQ=255)
                                ttCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCg                                                                   >  2:499047/1‑100 (MQ=255)
                                ttCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCg                                                                   <  2:2091104/100‑1 (MQ=255)
                                 tCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCgg                                                                  <  2:860245/100‑1 (MQ=255)
                                  cAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCgg                                                                  <  1:466521/99‑1 (MQ=255)
                                    gCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGAtt                                                               >  2:105957/1‑100 (MQ=255)
                                       ggggTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAAc                                                            <  1:341605/100‑1 (MQ=255)
                                       ggggTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAAc                                                            <  2:407201/100‑1 (MQ=255)
                                         ggTAAAAAGACATCGGCATTCAACTCCTGTGCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTaa                                                             >  2:1305806/1‑97 (MQ=255)
                                            aaaaaGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGTCCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATaa                                                       <  1:2123468/100‑1 (MQ=255)
                                            aaaaaGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATCCCGTCAGCCGGATTCACGATaa                                                       >  1:967189/1‑100 (MQ=255)
                                               aaGACAACGGCATTCACCTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTa                                                    <  1:752837/100‑1 (MQ=255)
                                                aGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGCCCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTaa                                                   <  2:2000600/100‑1 (MQ=255)
                                                 gACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACAGAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAAc                                                  <  1:687294/100‑1 (MQ=255)
                                                     tCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGAtt                                              >  2:704531/1‑100 (MQ=255)
                                                     tCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGAtt                                              <  1:105957/100‑1 (MQ=255)
                                                     tCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGAtt                                              <  2:29108/100‑1 (MQ=255)
                                                          aTTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTtcatc                                         >  1:1909304/1‑100 (MQ=255)
                                                             cAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACg                                      >  2:931902/1‑100 (MQ=255)
                                                             cAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTtcatca                                        >  1:711429/1‑98 (MQ=255)
                                                             cAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACg                                      >  2:1935533/1‑100 (MQ=255)
                                                               aCTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGAt                                    >  2:1229264/1‑100 (MQ=255)
                                                               aCTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGAt                                    <  1:902929/100‑1 (MQ=255)
                                                               aCTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGAt                                    <  2:2167203/100‑1 (MQ=255)
                                                                cTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATg                                   >  2:631015/1‑100 (MQ=255)
                                                                cTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATg                                   >  1:787940/1‑100 (MQ=255)
                                                                   cTGTTCTGAAGGGTGAGCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCt                                >  2:1481608/1‑100 (MQ=255)
                                                                   cTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCt                                <  1:396345/100‑1 (MQ=255)
                                                                   cTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCACCACGATGCCt                                >  2:191093/1‑100 (MQ=255)
                                                                     gTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTc                              >  2:2037027/1‑100 (MQ=255)
                                                                         tGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACTATTTCATCACGATGCCTTCGTAt                          <  2:494391/100‑1 (MQ=255)
                                                                         tGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTa                           <  1:219213/99‑1 (MQ=255)
                                                                         tGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTAt                          >  2:1191453/1‑100 (MQ=255)
                                                                          gAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATc                         >  2:732126/1‑100 (MQ=255)
                                                                          gAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATc                         <  1:1965114/100‑1 (MQ=255)
                                                                          gAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATc                         <  1:1964945/100‑1 (MQ=255)
                                                                            aGGGTGAGCGACCACACAGGCTATTTCACAAGTAAGGCAGCCAATACAGTCAGCCGGATTAACGAGAATTAACGATTTCATTACGATGCCTTCGTATCaa                       <  1:467164/100‑1 (MQ=255)
                                                                              ggTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCaa                       >  2:1400822/1‑98 (MQ=255)
                                                                              ggTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAAc                     <  1:582693/100‑1 (MQ=255)
                                                                               gTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACa                    >  2:1748574/1‑100 (MQ=255)
                                                                               gTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACa                    >  1:1735395/1‑100 (MQ=255)
                                                                                 gAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACaga                  <  2:2099004/100‑1 (MQ=255)
                                                                                   gCGACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACagag                 >  1:2167794/1‑99 (MQ=255)
                                                                                     gACCACACAGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTa              <  2:854962/100‑1 (MQ=255)
                                                                                       ccACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACCATTTCATCACGATGCCTTCGTATCAAACAGAGTTAAc            <  2:1197357/100‑1 (MQ=255)
                                                                                          cacaGGCTATTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTAACata         <  1:1639183/100‑1 (MQ=255)
                                                                                           acaGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTAACatat        <  1:2028859/100‑1 (MQ=255)
                                                                                              ggCTATTTCAAAAGTCCGGCAGAAAATACAGTCAGCAGGATCAACCATAATTAACGATTTCAGCACGATGCCTTCGTATCAAACAGAGTTAACatacctc     >  2:411776/1‑96 (MQ=255)
                                                                                               gCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTAACATATcg      <  2:1822896/98‑1 (MQ=255)
                                                                                                 tccTTCACAAGTAGGGCAGCCAATACAGTCAGCCGGATTAACAATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTAACATATcgcgcg  >  1:1516925/3‑100 (MQ=255)
                                                                                                   cTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATAAAACAGAGTTAATATATcgag    >  1:654707/1‑94 (MQ=255)
                                                                                                   |                                                                                                 
CCGGCGCACTGATGCTGTCCAGTCGCTGTACCTTCAGCCGGGGTAAAAAGACATCGGCATTCAACTCCTGTTCTGAAGGGTGAGCGACCACACAGGCTACTTCACAAGTACGGCAGCCAATACAGTCAGCCGGATTAACGATAATTAACGATTTCATCACGATGCCTTCGTATCAAACAGAGTTAACATATCGCGCG  >  NC_000913/3028856‑3029052

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