Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 898,949 A→G 7.1% N70S (AAC→AGC)  pspB → transcriptional regulator of psp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported898,9490AG 7.1% 205.4 / 5.5 70N70S (AAC→AGC) pspBtranscriptional regulator of psp operon
Reads supporting (aligned to +/- strand):  ref base A (27/38);  new base G (3/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCA  >  Exported/898883‑899016
                                                                  |                                                                   
atgaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACt                                                                   <  1:4454830/69‑1 (MQ=255)
atgaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTgg                                                                 <  1:2187039/71‑1 (MQ=255)
atgaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:2818099/70‑1 (MQ=255)
atgaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:2064135/70‑1 (MQ=255)
atgaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:440347/70‑1 (MQ=255)
atgaAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:1940459/70‑1 (MQ=255)
   aaGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTgg                                                                 <  1:3140463/68‑1 (MQ=255)
   aaGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTgg                                                                 <  1:2401599/68‑1 (MQ=255)
   aaGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTgg                                                                 <  1:1655957/68‑1 (MQ=255)
   aaGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAgg                                                              >  1:474326/1‑71 (MQ=255)
    aGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAGCTGGAggg                                                             <  1:3545162/71‑1 (MQ=255)
       aaaaCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATACGAACTGGAGGGATc                                                          >  1:4394312/1‑71 (MQ=255)
         aaCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:1444377/1‑71 (MQ=255)
          aCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:362264/1‑70 (MQ=255)
          aCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:3895440/71‑1 (MQ=255)
          aCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:3744327/71‑1 (MQ=255)
          aCGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:923133/71‑1 (MQ=255)
           cGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTa                                                      <  1:1163082/71‑1 (MQ=255)
               tGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATgg                                                  <  1:2760364/71‑1 (MQ=255)
               tGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATgg                                                  <  1:3468531/71‑1 (MQ=255)
                gcgcGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATgg                                                  <  1:2733132/70‑1 (MQ=255)
                  gcgAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:3534707/52‑1 (MQ=255)
                    gAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCggg                                              >  1:241169/1‑70 (MQ=255)
                     aaCGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGAt                                                           >  1:1780149/1‑56 (MQ=255)
                     aaCGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGAt                                                           >  1:1255080/1‑56 (MQ=255)
                     aaCGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGAt                                                           >  1:94700/1‑56 (MQ=255)
                     aaCGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAAt                                                    >  1:3307911/1‑63 (MQ=255)
                        gTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTa                                         <  1:2005812/71‑1 (MQ=255)
                        gTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTa                                         <  1:3062003/71‑1 (MQ=255)
                         tATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTaa                                        <  1:1666744/71‑1 (MQ=255)
                           ttCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:3129112/54‑1 (MQ=255)
                           ttCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:3876857/54‑1 (MQ=255)
                           ttCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:3973539/54‑1 (MQ=255)
                           ttCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:1981137/54‑1 (MQ=255)
                           ttCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCt                                                       <  1:1429011/54‑1 (MQ=255)
                            tCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:492963/42‑1 (MQ=255)
                             cAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAGCTGGAGGGATCGCTAATGGCGGGCATTAAtct                                     >  1:2148339/1‑70 (MQ=255)
                             cAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:3943609/1‑51 (MQ=255)
                             cAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:3166172/1‑51 (MQ=255)
                             cAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:1096544/1‑51 (MQ=255)
                             cAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:1119282/1‑51 (MQ=255)
                             cAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGc                                                        >  1:2697332/1‑51 (MQ=255)
                               ggCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCaa                                  >  1:3730194/1‑71 (MQ=255)
                                gcgcTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCaa                                  >  1:3045665/1‑70 (MQ=255)
                                gcgcTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCaa                                  >  1:2626465/1‑70 (MQ=255)
                                gcgcTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCaa                                  >  1:2539613/1‑70 (MQ=255)
                                  gcTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaa                               <  1:2825291/71‑1 (MQ=255)
                                  gcTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaa                               <  1:2045300/71‑1 (MQ=255)
                                  gcTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATa                                >  1:309350/1‑70 (MQ=255)
                                   cTGGAATCTATTCTTGATGCAGAACATCCGAACTg                                                                  <  1:2774369/35‑1 (MQ=255)
                                     ggAATCTATTCTTGATGCAGAACATCCGAACTGGa                                                                >  1:4219870/1‑35 (MQ=255)
                                     ggAATCTATTCTTGATGCAGAACATCCGAACTGGa                                                                >  1:2642631/1‑35 (MQ=255)
                                     ggAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGAt                                                           >  1:113228/1‑40 (MQ=255)
                                     ggAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaaa                             >  1:2297050/1‑70 (MQ=255)
                                      gAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaaaaa                           >  1:470531/1‑71 (MQ=255)
                                      gAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaaa                             <  1:1262735/69‑1 (MQ=255)
                                      gAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGAGCGCTAATGGCGGGCATTAATCTCAATaaaaa                            <  1:1069894/70‑1 (MQ=255)
                                       aaTCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATaaa                              >  1:4661884/1‑67 (MQ=255)
                                        aTCTATTCTTGATGCAGAACATGCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAAtt                         <  1:1179878/71‑1 (MQ=255)
                                        aTCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAAt                          >  1:2598312/1‑70 (MQ=255)
                                           tATTCTTGATGCAGAACATCCGAGCTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATg                      >  1:437494/1‑71 (MQ=255)
                                           tATTCTTGATGCAGAACATCCGAGCTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATg                      >  1:57823/1‑71 (MQ=255)
                                             ttCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTa                        <  1:1731967/67‑1 (MQ=255)
                                              tCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATg                      <  1:4431714/68‑1 (MQ=255)
                                                   aTGCAGAACATCCGAGCTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTc              <  1:3426431/71‑1 (MQ=255)
                                                   aTGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTc              <  1:3666089/71‑1 (MQ=255)
                                                         aaCATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAAt                          <  1:3026363/53‑1 (MQ=255)
                                                          aCATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCAcagca       <  1:3704146/71‑1 (MQ=255)
                                                             tCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCg                   <  1:3387210/56‑1 (MQ=255)
                                                               cGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCa  >  1:1097227/1‑71 (MQ=255)
                                                                  |                                                                   
ATGAAGCAAAACGGATGCGCGAACGTATTCAGGCGCTGGAATCTATTCTTGATGCAGAACATCCGAACTGGAGGGATCGCTAATGGCGGGCATTAATCTCAATAAAAAATTATGGCGTATTCCACAGCAGGGCA  >  Exported/898883‑899016

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: