Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,275 T→C 26.3% G911G (GGT→GGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,2750TC26.3% 122.7 / 46.5 76G911G (GGT→GGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (23/33);  new base C (10/10);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

TCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC  >  minE/834205‑834343
                                                                      |                                                                    
tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGGGAAAATCGGt                                                                      <  1:918393/71‑1 (MQ=255)
tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                                                      <  1:4618116/71‑1 (MQ=255)
tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                                                      <  1:2213556/71‑1 (MQ=255)
  gTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                                                    >  1:2423121/1‑71 (MQ=255)
    agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                                                  >  1:1718139/1‑71 (MQ=255)
    agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                                                  >  1:2292926/1‑71 (MQ=255)
     gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                                                                 <  1:2770130/71‑1 (MQ=255)
     gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                                                                 <  1:2331647/71‑1 (MQ=255)
      agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc                                                                <  1:2732854/71‑1 (MQ=255)
          cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                                             >  1:2191912/1‑70 (MQ=255)
             ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                                                    >  1:5053776/1‑60 (MQ=255)
             ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                                             >  1:2175296/1‑67 (MQ=255)
             ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAg                                                         >  1:4055910/1‑71 (MQ=255)
               aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                                                        <  1:3738083/70‑1 (MQ=255)
                  cGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                                                        <  1:4685043/67‑1 (MQ=255)
                   ggCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTg                                                   <  1:2770132/71‑1 (MQ=255)
                       aCTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAGGGATTGCCTgg                                              <  1:419840/71‑1 (MQ=255)
                         ttACAAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAggg                                                       <  1:3928344/61‑1 (MQ=255)
                           aCGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGaa                                            <  1:291860/70‑1 (MQ=255)
                            tgAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAAt                                                                   <  1:271263/45‑1 (MQ=255)
                              aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                                              <  1:267508/65‑1 (MQ=255)
                              aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                        >  1:4518783/1‑71 (MQ=255)
                                cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                                                      <  1:3049915/39‑1 (MQ=255)
                                cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                                                      <  1:1646611/39‑1 (MQ=255)
                                cGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAggg                                                       <  1:584154/54‑1 (MQ=255)
                                 gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGa                                                      <  1:2556674/54‑1 (MQ=255)
                                 gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGa                                                      <  1:2753986/54‑1 (MQ=255)
                                 gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                                              <  1:3178416/62‑1 (MQ=255)
                                  cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                                                    >  1:2996780/1‑39 (MQ=255)
                                  cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                                                    >  1:3443090/1‑39 (MQ=255)
                                  cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                                                    >  1:4420579/1‑39 (MQ=255)
                                  cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGCTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga                                    <  1:761213/71‑1 (MQ=255)
                                   tttACCTCTATCGGCCCGCGGATGGAGAAAATCGGTa                                                                     >  1:2821064/1‑37 (MQ=255)
                                    ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                                                  <  1:2732616/55‑1 (MQ=255)
                                    ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                                                  <  1:1780629/55‑1 (MQ=255)
                                     tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                                            >  1:2957626/1‑44 (MQ=255)
                                     tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc                                  <  1:3115838/70‑1 (MQ=255)
                                     tACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATACCCaggat                                  <  1:4847850/70‑4 (MQ=255)
                                      aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCg                                 <  1:5329022/70‑1 (MQ=255)
                                       ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                <  1:3160474/70‑1 (MQ=255)
                                       ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                <  1:2660199/70‑1 (MQ=255)
                                        ctctATCGGCCCGCTGATGGAGAAAATCGGTAATgg                                                                 <  1:2364127/36‑1 (MQ=255)
                                        ctctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                              <  1:4396547/71‑1 (MQ=255)
                                          ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc                                  <  1:210801/65‑1 (MQ=255)
                                           tATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg                            >  1:1766222/1‑70 (MQ=255)
                                            aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa                          >  1:894272/1‑71 (MQ=255)
                                            aTCGGCCCGCTGATGGAGAAAATAGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa                          >  1:4372801/1‑71 (MQ=255)
                                             tCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                              <  1:4481599/50‑1 (MQ=255)
                                             tCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGAt                         >  1:3797728/1‑71 (MQ=255)
                                              cGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATgg                           <  1:802142/68‑1 (MQ=255)
                                                   cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                              <  1:2845189/44‑1 (MQ=255)
                                                   cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                              <  1:4663557/44‑1 (MQ=255)
                                                   cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                              <  1:2551366/44‑1 (MQ=255)
                                                   cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                <  1:5435845/58‑1 (MQ=255)
                                                   cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                <  1:3581008/58‑1 (MQ=255)
                                                    gCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                              <  1:1579095/43‑1 (MQ=255)
                                                     cTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgt                                              <  1:4322964/42‑2 (MQ=255)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                            >  1:935884/1‑42 (MQ=16)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                            >  1:5430540/1‑42 (MQ=16)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                            >  1:5658164/1‑42 (MQ=16)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                            >  1:2891107/1‑42 (MQ=16)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGa                                             >  1:1536115/1‑41 (MQ=16)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGa                                             >  1:198260/1‑41 (MQ=16)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                        >  1:2030043/1‑46 (MQ=255)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                        >  1:515111/1‑46 (MQ=255)
                                                       gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                        >  1:2966700/1‑46 (MQ=255)
                                                         tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                       <  1:4887448/45‑1 (MQ=255)
                                                         tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                       <  1:1775653/45‑1 (MQ=255)
                                                         tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                       <  1:2677410/45‑1 (MQ=255)
                                                            agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt           >  1:5344887/1‑70 (MQ=255)
                                                            agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt           >  1:998717/1‑70 (MQ=255)
                                                             gaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCa         <  1:4824569/71‑1 (MQ=255)
                                                               aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc       <  1:205472/71‑1 (MQ=255)
                                                               aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc       <  1:792991/71‑1 (MQ=255)
                                                               aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc       >  1:1411784/1‑71 (MQ=255)
                                                                aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                              >  1:3745241/1‑47 (MQ=255)
                                                                aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                              >  1:2346196/1‑47 (MQ=255)
                                                                aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                <  1:3094082/45‑1 (MQ=255)
                                                                 aTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTa     <  1:5024182/71‑1 (MQ=255)
                                                                 aTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCAAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTa     >  1:665690/1‑71 (MQ=255)
                                                                    ggTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacac  >  1:3989957/1‑71 (MQ=255)
                                                                      |                                                                    
TCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC  >  minE/834205‑834343

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