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,5560AC81.8% 245.6 / 38.7 77L79* (TTA→TGA) sufDFe‑S cluster scaffold complex subunit SufD
Reads supporting (aligned to +/- strand):  ref base A (10/4);  new base C (33/30);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

GCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGC  >  NC_000913/1761509‑1761601
                                               |                                             
gCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCAc                                             >  1:1073438/1‑50 (MQ=255)
gCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCAc                                             <  1:4100563/50‑1 (MQ=255)
gCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCAc                                             >  1:2266621/1‑50 (MQ=255)
 cACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACg                                            <  1:3256621/50‑1 (MQ=255)
 cACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACg                                            >  1:7205707/1‑50 (MQ=255)
  aCGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGc                                           <  1:746885/50‑1 (MQ=255)
  aCGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGc                                           >  1:5247155/1‑50 (MQ=255)
   cGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGc                                           >  1:303897/1‑49 (MQ=255)
    gTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCt                                          >  1:2005687/1‑49 (MQ=255)
     tAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCt                                          >  1:6777830/1‑48 (MQ=255)
      aaCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                       >  1:4402556/1‑50 (MQ=255)
      aaCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                       >  1:3634782/1‑50 (MQ=255)
       aCGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAgg                                      <  1:3436997/50‑1 (MQ=255)
       aCGCCCATCCACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAg                                       <  1:2327926/49‑1 (MQ=39)
        cGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAgg                                      >  1:3118957/1‑49 (MQ=255)
        cGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGc                                     <  1:4700876/50‑1 (MQ=255)
         gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCa                                    >  1:1491717/1‑50 (MQ=255)
         gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACTCTAAGGCa                                    >  1:4682845/1‑50 (MQ=255)
         gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCa                                    <  1:6254145/50‑1 (MQ=255)
         gCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCa                                    >  1:3190428/1‑50 (MQ=255)
           ccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATc                                  <  1:5459034/50‑1 (MQ=255)
           ccATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATc                                  <  1:2304826/50‑1 (MQ=255)
            cATCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCa                                 >  1:3234266/1‑50 (MQ=255)
              tCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACg                               >  1:6742575/1‑50 (MQ=255)
              tCGACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACg                               <  1:154392/50‑1 (MQ=255)
              tCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACg                               >  1:7553387/1‑50 (MQ=255)
               cGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACgc                              >  1:4628831/1‑50 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                             >  1:3993358/1‑50 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                             >  1:4388150/1‑50 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgct                             >  1:6082183/1‑50 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgc                              <  1:249703/49‑1 (MQ=255)
                gACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACgct                             <  1:4293148/50‑1 (MQ=255)
                 aCAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctg                            >  1:4736816/1‑50 (MQ=255)
                 aCAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctg                            >  1:4867007/1‑50 (MQ=255)
                  cAAACACCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACgctgc                           <  1:3109392/50‑1 (MQ=255)
                      cacCAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCtgtg                       >  1:6887855/1‑50 (MQ=255)
                      cacCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCtgtg                       >  1:1673206/1‑50 (MQ=255)
                        ccAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTggg                     >  1:4850451/1‑50 (MQ=255)
                         cAGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTggg                     >  1:7121844/1‑49 (MQ=255)
                         cAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTggg                     >  1:5504850/1‑49 (MQ=255)
                          aGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGa                   <  1:1898060/50‑1 (MQ=255)
                          aGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGa                   >  1:7580741/1‑50 (MQ=255)
                          aGCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGa                   >  1:7194499/1‑50 (MQ=255)
                           gCCGCACGGAGTCTAACGTTCACGCTAGGGCATCACGCTGCTGTGGGGat                  <  1:3431402/50‑1 (MQ=255)
                           gCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGat                  <  1:2892817/50‑1 (MQ=255)
                           gCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGat                  <  1:6212764/50‑1 (MQ=255)
                           gCCGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGat                  <  1:2071521/50‑1 (MQ=255)
                            ccGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGata                 <  1:4099344/50‑1 (MQ=255)
                             cGCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGatat                >  1:4438288/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc               >  1:1553415/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc               >  1:2897103/1‑50 (MQ=255)
                              gCACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtc               <  1:5314328/50‑1 (MQ=255)
                               cACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtct              >  1:4417041/1‑50 (MQ=255)
                               cACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtct              >  1:6781650/1‑50 (MQ=255)
                               cACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtct              <  1:5448840/50‑1 (MQ=255)
                               cACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtct              >  1:4918322/1‑50 (MQ=255)
                               cACGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtct              <  1:4228006/50‑1 (MQ=255)
                                aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc             <  1:2028644/50‑1 (MQ=255)
                                aCGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctc             <  1:6745504/50‑1 (MQ=255)
                                 cGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctct            >  1:7200965/1‑50 (MQ=255)
                                 cGGAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctct            <  1:7045363/50‑1 (MQ=255)
                                 cGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATAtctct            >  1:6080691/1‑50 (MQ=255)
                                  ggAGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATAtctctc           <  1:2191281/50‑1 (MQ=255)
                                    aGTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTccc         <  1:4831471/50‑1 (MQ=255)
                                    aGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTccc         <  1:7560893/50‑1 (MQ=255)
                                     gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg        <  1:1214781/50‑1 (MQ=255)
                                     gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg        >  1:3536656/1‑50 (MQ=255)
                                     gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg        <  1:2657051/50‑1 (MQ=255)
                                     gTCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg        <  1:2019329/50‑1 (MQ=255)
                                      tCTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCg        >  1:5752440/1‑49 (MQ=255)
                                       cTAACGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGc       >  1:724089/1‑49 (MQ=255)
                                        tAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCaa     >  1:7395120/1‑50 (MQ=255)
                                          aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCaa     <  1:3978625/48‑1 (MQ=255)
                                          aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAAt    <  1:7601290/49‑1 (MQ=255)
                                          aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATg   <  1:5295923/50‑1 (MQ=255)
                                          aCGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATg   >  1:7544182/1‑50 (MQ=255)
                                           cGTTCACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGc  >  1:5698459/1‑50 (MQ=255)
                                               |                                             
GCACGTAACGCCCATCGACAAACACCAGCCGCACGGAGTCTAACGTTAACGCTAAGGCATCACGCTGCTGTGGGGATATCTCTCCCGCAATGC  >  NC_000913/1761509‑1761601

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