Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,017,118 C→A R353R (CGC→CGA ygfK → putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,017,1180CA93.8% 149.3 / 0.1 48R353R (CGC→CGAygfKputative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (25/20);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATCAACAATAAAGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGCGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTATTCGCCCTATTACT  >  NC_000913/3017023‑3017211
                                                                                               |                                                                                             
aTCAACAATAAAGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCt                                                                                           >  1:453833/1‑100 (MQ=255)
   aacaaTAAAGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTg                                                                                        >  2:256210/1‑100 (MQ=255)
   aacaaTAAAGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTg                                                                                        <  2:927823/100‑1 (MQ=255)
     caaTAAAGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGAc                                                                                      >  2:828469/1‑100 (MQ=255)
         aaaGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCa                                                                                  <  1:1230527/100‑1 (MQ=255)
              cgcACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTg                                                                             >  2:431009/1‑100 (MQ=255)
                cACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGcc                                                                           >  2:950986/1‑100 (MQ=255)
                    gCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCAtt                                                                       >  2:159905/1‑100 (MQ=255)
                        ggTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGTAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCtt                                                                   >  1:834359/1‑100 (MQ=255)
                           gAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTAtt                                                                >  1:954260/1‑100 (MQ=255)
                             agagATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGTAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTcc                                                              <  2:834353/100‑1 (MQ=255)
                                gaTGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCggt                                                           <  1:1188090/100‑1 (MQ=255)
                                     atatGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGc                                                      <  1:119766/100‑1 (MQ=255)
                                      tatGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGcc                                                     <  2:708820/100‑1 (MQ=255)
                                        tGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAg                                                   >  1:89613/1‑100 (MQ=255)
                                        tGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAg                                                   <  2:12114/100‑1 (MQ=255)
                                           cAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCa                                                <  1:76047/100‑1 (MQ=255)
                                           cAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCa                                                <  2:132991/100‑1 (MQ=255)
                                           cAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCa                                                >  2:1079611/1‑100 (MQ=255)
                                               ccGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTg                                            <  2:1173913/100‑1 (MQ=255)
                                                 gTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGAc                                          <  2:727332/100‑1 (MQ=255)
                                                     gcTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTa                                        >  1:304489/1‑98 (MQ=255)
                                                        gTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCcgc                                   <  1:1024422/100‑1 (MQ=255)
                                                        gTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCcgc                                   >  1:158998/1‑100 (MQ=255)
                                                         ttCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGCGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCcgcg                                  <  2:948904/100‑1 (MQ=255)
                                                         ttCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCcgcg                                  <  2:472472/100‑1 (MQ=255)
                                                           cccGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGat                                <  1:950994/100‑1 (MQ=255)
                                                           cccGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGat                                <  1:6365/100‑1 (MQ=255)
                                                            acGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGCGCCTTTTACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCcg                                    >  2:790972/2‑95 (MQ=255)
                                                             cGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGCGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGatat                              <  2:750348/100‑1 (MQ=255)
                                                             cGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGatat                              <  2:959205/100‑1 (MQ=255)
                                                               ctctCCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATAttt                            >  2:58073/1‑100 (MQ=255)
                                                                 ctcCATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATAttttt                          <  2:1270916/100‑1 (MQ=255)
                                                                   ccATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGa                        >  2:472104/1‑100 (MQ=255)
                                                                   ccATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGa                        <  1:803541/100‑1 (MQ=255)
                                                                    cATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGAt                       >  2:1928/1‑100 (MQ=255)
                                                                    cATCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGAt                       >  1:14703/1‑100 (MQ=255)
                                                                      tCAATGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATAc                     >  2:1205716/1‑100 (MQ=255)
                                                                        aaTGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAg                   <  2:158996/100‑1 (MQ=255)
                                                                         aTGTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAgg                  >  1:5216/1‑100 (MQ=255)
                                                                           gTTGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTa                >  2:702483/1‑100 (MQ=255)
                                                                             tGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTAtt              >  2:1145408/1‑100 (MQ=255)
                                                                             tGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTAtt              >  2:1047347/1‑100 (MQ=255)
                                                                             tGCAGCAGTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTAtt              >  1:737893/1‑100 (MQ=255)
                                                                               cagcagTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTATTCg            >  2:394539/1‑100 (MQ=255)
                                                                               cagcagTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTATTCg            >  2:557948/1‑100 (MQ=255)
                                                                                 gcagTTCTCTCTCGAGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTATTCGcc          <  1:260054/100‑1 (MQ=255)
                                                                                         ctctcGAGCCTTTGACGGCAAACTGCCTATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTATTCGCCCTATTACt  >  1:1027674/1‑100 (MQ=255)
                                                                                               |                                                                                             
ATCAACAATAAAGGCGCACTGCCTGGTGAAGAGATGTATATGTCAGGCCGTGCGCTGTTCCCGCTCTCCATCAATGTTGCAGCAGTTCTCTCTCGCGCCTTTGACGGCAAACTGCCCATTTCTTATTCCGGTGGTGCCAGTCAGCTGACTATCCGCGATATTTTTGATACAGGTATTCGCCCTATTACT  >  NC_000913/3017023‑3017211

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