Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,405,556 T→C intergenic (‑484/+147) nuoA ← / ← lrhA NADH:quinone oxidoreductase subunit A/DNA‑binding transcriptional dual regulator LrhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,405,5560TC81.1% 87.9 / 16.6 37intergenic (‑484/+147)nuoA/lrhANADH:quinone oxidoreductase subunit A/DNA‑binding transcriptional dual regulator LrhA
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (16/14);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

AGCAGGTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAATCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTATGC  >  NC_000913/2405509‑2405602
                                               |                                              
aGCAGGTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAAtct                                              >  1:3015089/1‑50 (MQ=255)
aGCAGGTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAAcct                                              >  1:1733549/1‑50 (MQ=255)
 gCAGGTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAAcct                                              <  1:3331674/49‑1 (MQ=255)
   aGGTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCtg                                           >  1:3474211/1‑50 (MQ=255)
    ggTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAATCTCtgg                                          <  1:2903983/50‑1 (MQ=255)
        aTTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTggtggt                                      >  1:2372757/1‑50 (MQ=255)
          ttttttCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGtt                                     >  1:3799923/1‑49 (MQ=255)
           tttttCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTg                                    <  1:1192169/49‑1 (MQ=255)
            ttttCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGa                                   >  1:3496102/1‑49 (MQ=255)
              ttCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGaaa                                 <  1:2847065/49‑1 (MQ=255)
               tCAGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGaaaa                                <  1:1687606/49‑1 (MQ=255)
                 aGGCGATTATTTGTGCGTTCGGGACGTGAATCTCTGGTGGTTGAAAAAt                              <  1:1357255/49‑1 (MQ=255)
                 aGGCGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATg                             <  1:476458/50‑1 (MQ=255)
                    cGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAAc                          <  1:3623514/50‑1 (MQ=255)
                    cGATTATTTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAAc                          <  1:1283577/50‑1 (MQ=255)
                     gATTATTTGTGCGTTCGGGACGTGAATCTCTGGTGGTTGAAAAATGaa                           >  1:2829975/1‑48 (MQ=255)
                      attattTGTGCGTTCGGGACGTGAATCTCTGGTGGTTGAAAAATGAACAg                        <  1:366922/50‑1 (MQ=255)
                         attTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGtt                      <  1:243077/49‑1 (MQ=255)
                         attTGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGtt                      <  1:385939/49‑1 (MQ=255)
                          tttGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGttt                     <  1:2143762/49‑1 (MQ=255)
                           ttGTGCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTg                   >  1:1394028/1‑50 (MQ=255)
                             gtgCGTTCGGGACGTGAATCTCTGGTGGTTGAAAAATGAACAGTTTTGTa                 >  1:2592215/1‑50 (MQ=255)
                             gtgCGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTa                 >  1:2410186/1‑50 (MQ=255)
                                cGTTCGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACg               <  1:3296463/49‑1 (MQ=255)
                                 gTTCGGGACGTGAATCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGtt             <  1:1136147/50‑1 (MQ=255)
                                   tCGGGACGTGAATCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCt           >  1:2894402/1‑50 (MQ=255)
                                    cGGGACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTg          >  1:640655/1‑50 (MQ=255)
                                     gggACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTg          >  1:2105063/1‑49 (MQ=255)
                                     gggACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGc         <  1:446092/50‑1 (MQ=255)
                                     gggACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGc         <  1:540616/50‑1 (MQ=255)
                                      ggACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCa        >  1:899442/1‑50 (MQ=255)
                                       gACGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCAc       >  1:1318451/1‑50 (MQ=255)
                                        aCGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCAc       >  1:2917859/1‑49 (MQ=255)
                                        aCGTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCAc       >  1:1081016/1‑49 (MQ=255)
                                          gcgAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTAt    <  1:845245/48‑1 (MQ=38)
                                          gTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTAt    >  1:25897/1‑50 (MQ=255)
                                          gTGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTAt    >  1:2000324/1‑50 (MQ=255)
                                           tGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTATg   >  1:2824312/1‑50 (MQ=255)
                                           tGAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTATg   <  1:2279471/50‑1 (MQ=255)
                                            gAACCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTATGc  >  1:1221574/1‑50 (MQ=255)
                                               |                                              
AGCAGGTGATTTTTTTCAGGCGATTATTTGTGCGTTCGGGACGTGAATCTCTGGTGGTTGAAAAATGAACAGTTTTGTACGTTCTGCACTATGC  >  NC_000913/2405509‑2405602

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: