Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,713,566 A→T 29.7% I265I (ATT→ATA 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,5660AT29.7% 100.9 / 53.3 84I265I (ATT→ATAnirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base A (24/35);  new base T (14/11);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GTTCCAGCTTGCGCATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGCGGACG  >  minE/2713503‑2713631
                                                               |                                                                 
gttCCAGCTTGCGCATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGAc                                                              <  1:1253037/69‑1 (MQ=255)
gttCCAGCTTGCGCATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGAc                                                              <  1:1671950/69‑1 (MQ=255)
             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:51841/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:1717911/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:674329/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1531136/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:243918/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1252598/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1589849/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1388240/1‑67 (MQ=255)
              aTTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCAcg                                                  >  1:1399015/1‑67 (MQ=255)
                 cACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGa                                              >  1:2939817/1‑68 (MQ=255)
                   ccGATTGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:84649/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:2433336/68‑1 (MQ=255)
                   ccGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:1190037/68‑1 (MQ=255)
                              gATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGc                                            <  1:2396073/57‑1 (MQ=255)
                                tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:2151579/68‑1 (MQ=255)
                                tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:208712/68‑1 (MQ=255)
                                tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:3138167/68‑1 (MQ=255)
                                tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:207206/68‑1 (MQ=255)
                                tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:2976942/68‑1 (MQ=255)
                                tatCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACca                               <  1:447488/68‑1 (MQ=255)
                                 atCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACcaca                             >  1:2914705/1‑69 (MQ=255)
                                   ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:515519/56‑1 (MQ=255)
                                   ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2616022/56‑1 (MQ=255)
                                   ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2706985/56‑1 (MQ=255)
                                   ccGGATCGGAAGTCTGGCAGGAATCATTTATGACGATACCCCCACGCGGAGCAACg                                        <  1:2538738/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:1097170/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)
                                    cGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTg                          <  1:2093811/69‑1 (MQ=255)
                                             aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:403707/1‑69 (MQ=255)
                                             aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:1421586/1‑69 (MQ=255)
                                             aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:452353/1‑69 (MQ=255)
                                             aGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAg                 >  1:1022963/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)
                                                               |                                                                 
GTTCCAGCTTGCGCATTCACCGATGGCGTAGATATCCGGATCGGAAGTCTGGCAGGAATCATTAATGACAATACCCCCACGCGGAGCAACGTCCAGACCACACTGGGTTGCCAGCTTATCGCGCGGACG  >  minE/2713503‑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: