Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,279,466 T→G 8.0% intergenic (‑271/‑68) narX ← / → narK sensory histidine kinase in two‑component regulatory system with NarL/nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,279,4660TG 8.0% 383.7 / 29.7 151intergenic (‑271/‑68)narX/narKsensory histidine kinase in two‑component regulatory system with NarL/nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base T (66/73);  new base G (6/6);  total (72/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 166 total aligned reads displayed.

ATGAAGTCACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTCTAT  >  AP009048/1279396‑1279535
                                                                      |                                                                     
aTGAAGTCACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGAtt                                                                       >  1:21945/1‑71 (MQ=255)
     gTCACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACAt                                                                  >  1:4563984/1‑71 (MQ=255)
     gTCACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACAt                                                                  >  1:4983163/1‑71 (MQ=255)
      tCACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGAtt                                                                       <  1:2738206/65‑1 (MQ=255)
       cACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATc                                                                 >  1:1337739/1‑70 (MQ=255)
        aCTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCaa                                                               <  1:1821126/71‑1 (MQ=255)
        aCTGTACTCACTATGGGAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCaaa                                                              <  1:1071537/71‑1 (MQ=255)
           gTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAAtt                                                            >  1:5755945/1‑71 (MQ=255)
            tACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAAtt                                                            <  1:1481850/70‑1 (MQ=255)
            tACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAAtt                                                            <  1:6543877/70‑1 (MQ=255)
              cTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGcc                                                         <  1:784670/71‑1 (MQ=255)
              cTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGcc                                                         <  1:2128235/71‑1 (MQ=255)
                 aCTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCttt                                                      <  1:1917836/71‑1 (MQ=255)
                 aCTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCtt                                                       >  1:2954472/1‑70 (MQ=255)
                 aCTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCtt                                                       <  1:7335596/70‑1 (MQ=255)
                 aCTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCtt                                                       <  1:7584521/70‑1 (MQ=255)
                   tATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCaaa                                                              <  1:2074588/61‑1 (MQ=255)
                   tATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCtt                                                       >  1:7507940/1‑68 (MQ=255)
                   tATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTa                                                     >  1:4162995/1‑70 (MQ=255)
                   tATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAg                                                    >  1:4541583/1‑71 (MQ=255)
                    aTGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTa                                                                     >  1:3013539/1‑53 (MQ=255)
                          aaTGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCa                                                                <  1:112780/52‑1 (MQ=255)
                          aaTGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGa                                             >  1:1734641/1‑71 (MQ=255)
                          aaTGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGa                                             >  1:5398859/1‑71 (MQ=255)
                           aTGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGAttt                                                                      <  1:2438365/45‑1 (MQ=255)
                           aTGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGa                                             <  1:4882866/70‑1 (MQ=255)
                             gATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGacac                                          >  1:3344603/1‑71 (MQ=255)
                             gATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGa                                             >  1:3454490/1‑68 (MQ=255)
                              aTAAATATCAATGATAGATAAAGTTATCTTATCGTTTGAttt                                                                      >  1:4644464/1‑42 (MQ=255)
                              aTAAATATCAATGATAGATAAAGTTATCTTATCGTTTGAttt                                                                      >  1:694314/1‑42 (MQ=255)
                              aTAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACa                                                                   <  1:2849913/45‑1 (MQ=255)
                              aTAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCttt                                                      <  1:1916933/58‑1 (MQ=255)
                               tAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACAt                                                                  >  1:6465088/1‑45 (MQ=255)
                               tAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACAt                                                                  >  1:4349614/1‑45 (MQ=255)
                               tAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTa                                        <  1:2872055/71‑1 (MQ=255)
                                aaaTATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACt                                         >  1:6475655/1‑69 (MQ=255)
                                 aaTATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAg                                      >  1:3508035/1‑71 (MQ=255)
                                 aaTATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAg                                      >  1:761979/1‑71 (MQ=255)
                                  atatCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACt                                         >  1:4389769/1‑67 (MQ=255)
                                  atatCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAgg                                     >  1:6878939/1‑71 (MQ=255)
                                  atatCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAgg                                     <  1:4018811/71‑1 (MQ=255)
                                   tatCAATGATAGATAAAGTTATCTTATCGTTTGAtt                                                                       <  1:2357023/36‑1 (MQ=255)
                                   tatCAATGATAGATAAAGTTATCTTATCGTTTGAtt                                                                       <  1:4067562/36‑1 (MQ=255)
                                   tatCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCttt                                                      <  1:3119498/53‑1 (MQ=255)
                                   tatCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCttt                                                      <  1:2533486/53‑1 (MQ=255)
                                   tatCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGt                                    <  1:5584121/71‑1 (MQ=255)
                                    atCAATGATAGATAAAGTTATCTTATCGTTTGAtt                                                                       <  1:1369394/35‑1 (MQ=255)
                                     tCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTa                                                 >  1:5768191/1‑56 (MQ=255)
                                     tCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAgg                                     >  1:3853709/1‑68 (MQ=255)
                                     tCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTgg                                  <  1:4395294/71‑1 (MQ=255)
                                       aaTGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGc                                 <  1:97197/70‑1 (MQ=255)
                                       aaTGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGc                                 <  1:4379411/70‑1 (MQ=255)
                                       aaTGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCa                                <  1:6091976/71‑1 (MQ=255)
                                        aTGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAAt                                                             >  1:3233500/1‑41 (MQ=255)
                                           atagattAAGTTATCTTATCGTTTGATTTACATCaaa                                                              >  1:15324/1‑37 (MQ=255)
                                           atagataAAGTTATCTTATCGTTTGATTTACATCaaa                                                              >  1:3377021/1‑37 (MQ=255)
                                           atagataAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACa                            >  1:4761964/1‑71 (MQ=255)
                                           atagataAAGTTATCTTATCATTTGATGTAAATCaa                                                               <  1:888061/36‑1 (MQ=21)
                                           atagataAAGTTATCTTATCATTTGATGTAAATCaa                                                               <  1:5196008/36‑1 (MQ=21)
                                           atagataAAGTTATCTTATCATTTGATGTAAATCaa                                                               <  1:5481752/36‑1 (MQ=21)
                                           atagataAAGTTATCTTATCATTTGATGTAAATCaa                                                               <  1:3730957/36‑1 (MQ=21)
                                             agataAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGa                                             >  1:3364910/1‑52 (MQ=255)
                                             agataAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTgg                                  <  1:115360/63‑1 (MQ=255)
                                             agataAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTgg                                  <  1:5649029/63‑1 (MQ=255)
                                              gataAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATc                          <  1:1958205/70‑1 (MQ=255)
                                               ataAAGTTATCTTATCGTTTGATTTACATCAAATTGc                                                          <  1:4466135/37‑1 (MQ=255)
                                               ataAAGTTATCTTATCGTTTGATTTACATCAAATTGc                                                          <  1:2809516/37‑1 (MQ=255)
                                                 aaaGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAg                                                    >  1:2047030/1‑41 (MQ=255)
                                                 aaaGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGaaa                      >  1:1731952/1‑71 (MQ=255)
                                                  aaGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGc                                                   >  1:1735172/1‑41 (MQ=255)
                                                   aGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAg                                      >  1:4637983/1‑53 (MQ=255)
                                                   aGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGAc                             <  1:5202755/62‑1 (MQ=255)
                                                   aGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGAc                             <  1:2614641/62‑1 (MQ=255)
                                                    gTTATCTTATCGTTTGATTTACATCAAATTGCCttt                                                      <  1:480836/36‑1 (MQ=255)
                                                    gTTATCTTATCGTTTGATTTACATCAAATTGCCttt                                                      <  1:2418886/36‑1 (MQ=255)
                                                     ttatcttatcGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAAc                     <  1:2776873/68‑1 (MQ=255)
                                                           tatcGTTTGATTTACATCAAATTGCCTTTAGCTACAGaca                                           >  1:1206328/1‑40 (MQ=255)
                                                             tcGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCagag          >  1:7419902/1‑71 (MQ=255)
                                                             tcGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCagag          <  1:4539736/71‑1 (MQ=255)
                                                             tcGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCagag          <  1:7183335/71‑1 (MQ=255)
                                                              cGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCagag          <  1:7396375/70‑1 (MQ=255)
                                                               gTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACAt                           >  1:4938391/1‑52 (MQ=255)
                                                               gTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGgt        <  1:5075429/71‑1 (MQ=255)
                                                               gTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGgt        <  1:7241959/71‑1 (MQ=255)
                                                                tttGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGaa                       <  1:2927810/55‑1 (MQ=255)
                                                                tttGATGTAAATCAAATTGCCTTTAGCTACAGACACTAAGGTgg                                  >  1:6637293/1‑44 (MQ=38)
                                                                tttGATGTAAATCAAATTGCCTTTAGCTACAGACACTAAGGTgg                                  >  1:2271017/1‑44 (MQ=38)
                                                                tttGATGTAAATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGgtg       >  1:637673/1‑71 (MQ=255)
                                                                tttGATGTAAATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGgtg       >  1:7817110/1‑71 (MQ=255)
                                                                  tGATTTACATCAAATTGCCTTTAGCTACAGCCACTAAGGTGGCAGACATCGAAACGAGTAt               >  1:406718/1‑61 (MQ=255)
                                                                  tGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCa                                >  1:4712626/1‑44 (MQ=255)
                                                                  tGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGgtgt      <  1:701389/70‑1 (MQ=255)
                                                                   gATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAg                               >  1:6757438/1‑44 (MQ=255)
                                                                   gATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCagag          <  1:3406408/65‑1 (MQ=255)
                                                                   gATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTc     >  1:7163622/1‑70 (MQ=255)
                                                                   gATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTc     <  1:5923115/70‑1 (MQ=255)
                                                                    aTTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTCt    <  1:5797938/70‑1 (MQ=255)
                                                                     tttACATCAAATTGCCTTTAGCTACAGACACTAAGGTgg                                  <  1:385662/39‑1 (MQ=255)
                                                                     tttACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTCTAt  <  1:1146692/71‑1 (MQ=255)
                                                                      ttACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTCTAt  <  1:1429172/70‑1 (MQ=255)
                                                                      |                                                                     
ATGAAGTCACTGTACTCACTATGGGTAATGATAAATATCAATGATAGATAAAGTTATCTTATCGTTTGATTTACATCAAATTGCCTTTAGCTACAGACACTAAGGTGGCAGACATCGAAACGAGTATCAGAGGTGTCTAT  >  AP009048/1279396‑1279535

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