Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,713,572 A→G 30.4% I263I (ATT→ATC nirB ← nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,713,5720AG30.4% 74.0 / 37.7 82I263I (ATT→ATCnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base A (24/33);  new base G (14/11);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGCGGACG  >  minE/2713516‑2713631
                                                        |                                                           
cATTCACCGATGGCGTAGATATCCGTATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCAcg                                                  <  1:785391/68‑1 (MQ=255)
cATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCAcg                                                  <  1:1195662/68‑1 (MQ=255)
cATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCAcg                                                  <  1:1726519/68‑1 (MQ=255)
cATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCAcg                                                  <  1:552212/68‑1 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1717911/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:2643480/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:2081316/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:2914001/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1975817/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:43805/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:51841/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:674329/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1589849/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1399015/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1388240/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1531136/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1252598/1‑67 (MQ=255)
 aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:243918/1‑67 (MQ=255)
    cACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGa                                              >  1:2939817/1‑68 (MQ=255)
      ccGATTGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:84649/68‑1 (MQ=255)
      ccGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:2433336/68‑1 (MQ=255)
      ccGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:2420086/68‑1 (MQ=255)
      ccGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:1268211/68‑1 (MQ=255)
      ccGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:333502/68‑1 (MQ=255)
      ccGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:1190037/68‑1 (MQ=255)
                 gATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:2396073/57‑1 (MQ=255)
                   tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:2976942/68‑1 (MQ=255)
                   tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:2151579/68‑1 (MQ=255)
                   tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:208712/68‑1 (MQ=255)
                   tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:207206/68‑1 (MQ=255)
                   tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:3138167/68‑1 (MQ=255)
                   tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:447488/68‑1 (MQ=255)
                    atCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACcaca                             >  1:2914705/1‑69 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2616022/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:1097170/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2706985/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:515519/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:791166/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:887489/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:90896/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:1472002/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2265327/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:1648169/56‑1 (MQ=255)
                      ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2538738/56‑1 (MQ=255)
                       cGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTg                          <  1:2093811/69‑1 (MQ=255)
                                aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:403707/1‑69 (MQ=255)
                                aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:1022963/1‑69 (MQ=255)
                                aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:1421586/1‑69 (MQ=255)
                                aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:452353/1‑69 (MQ=255)
                                 gTCTTGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:3233622/41‑1 (MQ=255)
                                  tCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCt               <  1:822975/69‑1 (MQ=255)
                                  tCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCt               <  1:1878508/69‑1 (MQ=255)
                                  tCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCt               <  1:1989074/69‑1 (MQ=255)
                                   cTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCt               <  1:2497175/68‑1 (MQ=255)
                                    tGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTa             >  1:3081235/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:1992323/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:2820135/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:1827444/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:739262/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      <  1:902076/69‑1 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:632990/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:1742218/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:1373507/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:1739165/1‑69 (MQ=255)
                                           aaTCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATcgcgcg      >  1:1682144/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:71010/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:1605414/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:418466/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:524593/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:1764228/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:2722188/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:2014278/1‑69 (MQ=255)
                                             tCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcgga    >  1:2834040/1‑69 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:2609711/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:2550228/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:1206712/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:2659475/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:532712/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:2677281/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:1822482/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:3302696/69‑1 (MQ=255)
                                               aTTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGcggacg  <  1:3106898/69‑1 (MQ=255)
                                                        |                                                           
CATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGCGGACG  >  minE/2713516‑2713631

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