Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,833 T→C 34.2% G1097G (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,8330TC34.2% 79.4 / 63.8 76G1097G (GGT→GGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base C (12/14);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCT  >  minE/834771‑834897
                                                              |                                                                
cAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactc                                                          >  1:2093182/1‑71 (MQ=255)
     aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc                                                      >  1:4653409/1‑70 (MQ=255)
        cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTa                                                    >  1:5375499/1‑69 (MQ=255)
           cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAg                                                 >  1:2368751/1‑69 (MQ=255)
           cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc                                                >  1:4637303/1‑70 (MQ=255)
            ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa                                              >  1:2051331/1‑71 (MQ=255)
              cAACCAGGAAAAAGCGCTCAACTTTCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa                                            <  1:3852092/71‑1 (MQ=255)
               aaCCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa                                            <  1:2695136/70‑1 (MQ=255)
                  ccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc                                         <  1:4368715/68‑1 (MQ=255)
                  cAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACaa                                           >  1:1402299/1‑68 (MQ=255)
                   cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                          <  1:567457/67‑1 (MQ=255)
                   cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                          <  1:189225/67‑1 (MQ=255)
                   cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                          <  1:3157486/67‑1 (MQ=255)
                   aGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctg                                       >  1:3144394/1‑71 (MQ=255)
                     gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca                                                             >  1:2401833/1‑47 (MQ=255)
                     gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                     <  1:3603640/71‑1 (MQ=255)
                     gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                     <  1:916451/71‑1 (MQ=255)
                     gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                     <  1:1125150/71‑1 (MQ=255)
                     gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTAGAGTGAAAACct                                        >  1:2600314/1‑68 (MQ=255)
                     gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAgtg                                               <  1:2542890/61‑3 (MQ=255)
                       aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga                                   >  1:4887880/1‑71 (MQ=255)
                        aaaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCAttcact                                                           >  1:2436063/1‑46 (MQ=21)
                        aaaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaat                                  >  1:4634689/1‑71 (MQ=255)
                         aaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc                                                               <  1:666464/41‑1 (MQ=255)
                         aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg                                 <  1:5477393/71‑1 (MQ=255)
                         aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg                                 <  1:3896675/71‑1 (MQ=255)
                          aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca                                                             >  1:3706156/1‑42 (MQ=255)
                          aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca                                                             >  1:1472731/1‑42 (MQ=255)
                          aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcca                                                             >  1:4597606/1‑42 (MQ=255)
                          aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc                                >  1:305566/1‑71 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                       >  1:5211438/1‑47 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                       >  1:5217717/1‑47 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                       >  1:1098225/1‑47 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactccac                                                       >  1:839167/1‑47 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccactcca                                                        >  1:1030869/1‑46 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                               <  1:778149/71‑1 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                               <  1:3647857/71‑1 (MQ=255)
                           gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                               <  1:4230928/71‑1 (MQ=255)
                            cgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                               >  1:3873271/1‑70 (MQ=255)
                              cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc                            >  1:3654334/4‑71 (MQ=255)
                              cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc                            >  1:132462/4‑71 (MQ=255)
                               ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACcttct                                     <  1:2322660/59‑1 (MQ=255)
                               ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                           >  1:2037702/3‑70 (MQ=255)
                               ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                           >  1:2673952/3‑70 (MQ=255)
                                 aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGCCtt                          >  1:456030/1‑70 (MQ=255)
                                 aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                         <  1:5208427/71‑1 (MQ=255)
                                 aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                         <  1:2071783/71‑1 (MQ=255)
                                 aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                         <  1:3100376/71‑1 (MQ=255)
                                 aaCTTCCTGACGCCGCACCAGAAATGGGGCAttcactc                                                          >  1:236988/1‑38 (MQ=25)
                                 aaCTTCCTGACGCCGCACCAGAAATGGGGCAttcactc                                                          >  1:5025560/1‑38 (MQ=25)
                                 aaCTTCCTGACGCCGCACCAGAAATGGGGCAttcactc                                                          >  1:1503412/1‑38 (MQ=25)
                                 aaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                         >  1:270191/1‑68 (MQ=255)
                                  aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                          >  1:1076048/1‑69 (MQ=255)
                                  aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                          >  1:1695595/1‑69 (MQ=255)
                                   cTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAg                       >  1:5416809/1‑71 (MQ=255)
                                     tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGGCa                                            <  1:2015197/48‑1 (MQ=255)
                                     tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGt                     >  1:3679281/1‑71 (MQ=255)
                                     tCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGaaa                                            <  1:4457032/48‑3 (MQ=19)
                                       cTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg                   <  1:1511231/71‑1 (MQ=255)
                                         gACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgctaa                                   <  1:3006596/53‑3 (MQ=255)
                                          aCGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg                <  1:4046506/71‑1 (MQ=255)
                                          aCGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                         <  1:1805862/62‑4 (MQ=255)
                                           cgccgcACCAGAAATGGGGCATTCAGTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                         <  1:1920594/61‑4 (MQ=255)
                                            gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                     <  1:1119767/48‑1 (MQ=19)
                                            gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                     <  1:1606080/48‑1 (MQ=19)
                                             ccgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtgg             <  1:4735454/71‑1 (MQ=255)
                                              cgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggt            >  1:2559416/1‑71 (MQ=255)
                                               gcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc           >  1:4439869/1‑71 (MQ=255)
                                                 aCCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg         <  1:368621/71‑1 (MQ=255)
                                                   cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                 <  1:5671764/45‑1 (MQ=255)
                                                   cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                 <  1:434018/45‑1 (MQ=255)
                                                   cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                 <  1:3101428/45‑1 (MQ=255)
                                                    aGAAGTAGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg                <  1:3939067/61‑1 (MQ=255)
                                                     gAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAgg                      >  1:2987148/1‑54 (MQ=255)
                                                        gTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCt  >  1:1825984/1‑71 (MQ=255)
                                                         tGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                               <  1:5218110/41‑1 (MQ=255)
                                                              tATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc           <  1:4076261/56‑1 (MQ=255)
                                                              |                                                                
CAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCT  >  minE/834771‑834897

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: