Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,761,556 A→C L79* (TTA→TGA)  sufD ← Fe‑S cluster scaffold complex subunit SufD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,761,5560AC85.9% 262.3 / 29.2 78L79* (TTA→TGA) sufDFe‑S cluster scaffold complex subunit SufD
Reads supporting (aligned to +/- strand):  ref base A (5/6);  new base C (43/24);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCTG  >  NC_000913/1761509‑1761603
                                               |                                               
gCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCAc                                               <  1:2904935/50‑1 (MQ=255)
 cACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACg                                              >  1:4650087/1‑50 (MQ=255)
 cACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACg                                              <  1:1647260/50‑1 (MQ=255)
 cACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACg                                              <  1:1508054/50‑1 (MQ=255)
  aCGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGc                                             >  1:2320738/1‑50 (MQ=255)
   cGTAACGCCCATCGACAAACACCAGCCGCACGTAGTCTAACGTTCACGCt                                            <  1:4264452/50‑1 (MQ=255)
   cGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGc                                             >  1:3422307/1‑49 (MQ=255)
   cGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGGTCACGCt                                            <  1:4514220/50‑1 (MQ=255)
    gTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCt                                            >  1:3241480/1‑49 (MQ=255)
     tAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCt                                            >  1:7088731/1‑48 (MQ=255)
     tAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCt                                            >  1:3181119/1‑48 (MQ=255)
     tAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTaa                                          <  1:5008605/50‑1 (MQ=255)
       aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAgg                                        >  1:2903110/1‑50 (MQ=255)
       aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAgg                                        <  1:6214800/50‑1 (MQ=255)
       aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                         <  1:3804272/49‑1 (MQ=255)
        cGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAgg                                        >  1:1366318/1‑49 (MQ=255)
         gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCa                                      >  1:6912591/1‑50 (MQ=255)
          cccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCAt                                     >  1:1698185/1‑50 (MQ=255)
          cccATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCAt                                     >  1:2237672/1‑50 (MQ=255)
           ccATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATc                                    <  1:635308/50‑1 (MQ=255)
            cATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCa                                   <  1:1537849/50‑1 (MQ=255)
             aTCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCAc                                  >  1:4898691/1‑50 (MQ=255)
             aTCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCAc                                  >  1:3720364/1‑50 (MQ=255)
             aTCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCAc                                  <  1:7017179/50‑1 (MQ=255)
             aTCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCAc                                  >  1:466549/1‑50 (MQ=255)
             aTCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCAc                                  <  1:7477954/50‑1 (MQ=255)
             aTCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCAc                                  <  1:7347616/50‑1 (MQ=255)
              tCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACg                                 <  1:5146481/50‑1 (MQ=255)
              tCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACg                                 >  1:1524850/1‑50 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                               >  1:4570353/1‑50 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgc                                <  1:3504485/49‑1 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACgct                               <  1:6853985/50‑1 (MQ=255)
                 aCAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctg                              >  1:5980223/1‑50 (MQ=255)
                 aCAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctg                              <  1:195719/50‑1 (MQ=255)
                 aCAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctg                              >  1:6886394/1‑50 (MQ=255)
                  cAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgc                             <  1:4010614/50‑1 (MQ=255)
                  cAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgc                             >  1:8164955/1‑50 (MQ=255)
                   aaaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgct                            >  1:3744908/1‑50 (MQ=255)
                   aaaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgct                            >  1:5908799/1‑50 (MQ=255)
                   aaaCACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgct                            >  1:6154949/1‑50 (MQ=255)
                      cacCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCtgtg                         >  1:3045958/1‑50 (MQ=255)
                        ccAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTggg                       <  1:5411867/50‑1 (MQ=255)
                         cAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTggg                       >  1:8298637/1‑49 (MQ=255)
                          aGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGa                     >  1:388124/1‑50 (MQ=255)
                          aGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGa                     >  1:2379378/1‑50 (MQ=255)
                           gCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGat                    <  1:5253495/50‑1 (MQ=255)
                           gCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGat                    <  1:765151/50‑1 (MQ=255)
                            ccGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGata                   <  1:5411119/50‑1 (MQ=255)
                            ccGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGata                   <  1:1412906/50‑1 (MQ=255)
                             cGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGatat                  >  1:7299995/1‑50 (MQ=255)
                             cGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGatat                  >  1:7670914/1‑50 (MQ=255)
                             cGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGatat                  >  1:3518676/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                 >  1:2038911/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                 <  1:1607624/50‑1 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                 >  1:2272040/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                 >  1:1582480/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc                 >  1:1091457/1‑50 (MQ=255)
                               cACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATAtct                >  1:3980561/1‑50 (MQ=255)
                                aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc               >  1:5867793/1‑50 (MQ=255)
                                aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc               >  1:2031643/1‑50 (MQ=255)
                                 cGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctct              >  1:5191136/1‑50 (MQ=255)
                                 cGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctct              >  1:2222323/1‑50 (MQ=255)
                                  ggAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctctc             >  1:652144/1‑50 (MQ=255)
                                   gAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTcc            >  1:1670656/1‑50 (MQ=255)
                                    aGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTccc           >  1:670652/1‑50 (MQ=255)
                                     gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg          >  1:7171426/1‑50 (MQ=255)
                                     gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg          >  1:7253794/1‑50 (MQ=255)
                                       cTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGc         >  1:3283332/1‑49 (MQ=255)
                                       cTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGc         >  1:7646087/1‑49 (MQ=255)
                                         aaCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAAt      >  1:5324331/1‑50 (MQ=255)
                                          aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCaa       <  1:5295959/48‑1 (MQ=255)
                                          aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATg     <  1:1639165/50‑1 (MQ=255)
                                           cGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGc    <  1:494087/50‑1 (MQ=255)
                                           cGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGc    >  1:6275497/1‑50 (MQ=255)
                                           cGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGc    <  1:7802430/50‑1 (MQ=255)
                                           cGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGc    <  1:5172135/50‑1 (MQ=255)
                                           cGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGc    >  1:100443/1‑50 (MQ=255)
                                             ttAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCTg  <  1:6172329/50‑1 (MQ=255)
                                               |                                               
GCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGCTG  >  NC_000913/1761509‑1761603

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