Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,065,001 A→G 9.3% K414K (AAA→AAG rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,065,0010AG 9.3% 335.2 / 19.0 118K414K (AAA→AAGrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (65/42);  new base G (8/3);  total (73/45)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 133 total aligned reads displayed.

TACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGT  >  Exported/2064939‑2065067
                                                              |                                                                  
  cacCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGAt                                                           <  1:2039917/70‑1 (MQ=255)
  cacCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                          <  1:1647461/71‑1 (MQ=255)
   acCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGAt                                                           <  1:4790281/69‑1 (MQ=255)
       aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                          <  1:3518716/66‑1 (MQ=255)
       aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                          <  1:1577864/66‑1 (MQ=255)
        ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                          <  1:2650053/65‑1 (MQ=255)
        ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                          <  1:2718564/65‑1 (MQ=255)
          gTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:2044388/1‑53 (MQ=255)
           tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                                 <  1:1285275/71‑1 (MQ=255)
           tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                                 <  1:2268039/71‑1 (MQ=255)
           tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                                 <  1:2272023/71‑1 (MQ=255)
           tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                                 <  1:2571706/71‑1 (MQ=255)
           tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                                 <  1:1216528/71‑1 (MQ=255)
             gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                               >  1:1575573/1‑71 (MQ=255)
             gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                               >  1:1623571/1‑71 (MQ=255)
             gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                               >  1:3951664/1‑71 (MQ=255)
             gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                               >  1:1420719/1‑71 (MQ=255)
             gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                               >  1:4756893/1‑71 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:2682547/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:4069536/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:2631084/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:4511602/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:3265514/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:3760790/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:2003139/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:1433436/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                    >  1:1786136/1‑40 (MQ=255)
                       ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                      >  1:4955241/1‑70 (MQ=255)
                        cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                    <  1:924103/55‑1 (MQ=255)
                        cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                    <  1:4709308/55‑1 (MQ=255)
                          tGACCTGATTCAGGAAGGCAAAATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTAc                                   >  1:1820771/1‑70 (MQ=255)
                            aCCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                <  1:2093895/71‑1 (MQ=255)
                             ccTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                <  1:2491365/70‑1 (MQ=255)
                             ccTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                               <  1:1108804/71‑3 (MQ=12)
                                gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGt                                         >  1:2818843/1‑58 (MQ=255)
                                gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                 >  1:491083/1‑66 (MQ=255)
                                  ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                      <  1:1495372/59‑1 (MQ=255)
                                  ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                          <  1:4621225/71‑1 (MQ=255)
                                  ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                          <  1:948051/71‑1 (MQ=255)
                                  ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                          <  1:1625803/71‑1 (MQ=255)
                                   tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                        >  1:60001/1‑40 (MQ=255)
                                   tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                        >  1:5027191/1‑40 (MQ=255)
                                   tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                        >  1:4147094/1‑40 (MQ=255)
                                   tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                        >  1:2681223/1‑40 (MQ=255)
                                   ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                         >  1:3173542/2‑69 (MQ=255)
                                   ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                         >  1:2428536/2‑69 (MQ=255)
                                   ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                         >  1:3170143/2‑69 (MQ=255)
                                   acAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGt                                      >  1:1878456/2‑58 (MQ=255)
                                       gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                  <  1:2127057/58‑1 (MQ=255)
                                        aaGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                  <  1:3690460/57‑1 (MQ=255)
                                        aaGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTAc                    >  1:2511348/1‑71 (MQ=255)
                                         aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCtg                     >  1:1743491/1‑68 (MQ=255)
                                         aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                     >  1:5043968/1‑69 (MQ=255)
                                         aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                     >  1:2254085/1‑69 (MQ=255)
                                         aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTATCCACCTa                     >  1:1578098/1‑69 (MQ=255)
                                              aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc              <  1:3107784/71‑1 (MQ=255)
                                              aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc              >  1:991692/1‑71 (MQ=255)
                                               taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                              >  1:75893/2‑54 (MQ=255)
                                               taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                              >  1:4970356/2‑54 (MQ=255)
                                               taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                              >  1:2656642/2‑54 (MQ=255)
                                               cATCGTTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt             <  1:483579/71‑1 (MQ=255)
                                               cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt             <  1:4866391/71‑1 (MQ=255)
                                               cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt             <  1:483580/71‑1 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                 >  1:742241/1‑50 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                 >  1:2881537/1‑50 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                 >  1:3728995/1‑50 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                 >  1:3117681/1‑50 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                                 >  1:2774363/1‑50 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                  >  1:72336/1‑49 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                           >  1:3423525/1‑56 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                           >  1:1387133/1‑56 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                           >  1:3295095/1‑56 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc                  >  1:1390675/1‑65 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc                  >  1:2513981/1‑65 (MQ=255)
                                                aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATAACGCCGTGGTTACAAGTtctc                           >  1:5133367/1‑56 (MQ=255)
                                                 tCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCtgg           <  1:3119972/71‑1 (MQ=255)
                                                  cGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                                  <  1:126012/47‑1 (MQ=255)
                                                      cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                      >  1:4496882/1‑55 (MQ=255)
                                                      cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGAt      >  1:1189773/1‑71 (MQ=255)
                                                      cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTATCCACCTACGCAACCTGGTGGAt      >  1:5067196/1‑71 (MQ=255)
                                                       tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                      <  1:1532944/38‑1 (MQ=255)
                                                       tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                      <  1:3878127/38‑1 (MQ=255)
                                                       tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                      <  1:2316836/38‑1 (MQ=255)
                                                       tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATc     >  1:3301009/1‑71 (MQ=255)
                                                         atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGtt                              >  1:4831013/1‑44 (MQ=255)
                                                         atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGt                               >  1:4011548/1‑43 (MQ=255)
                                                         atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATcc    <  1:2333211/70‑1 (MQ=255)
                                                         atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATcc    >  1:3176972/1‑70 (MQ=255)
                                                         atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATcc    <  1:2991830/70‑1 (MQ=255)
                                                         atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCg   >  1:523219/1‑71 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:2834245/1‑40 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:4143497/1‑40 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:4428889/1‑40 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:2001162/1‑40 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:2045447/1‑40 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:2190888/1‑40 (MQ=255)
                                                           gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                >  1:4738799/1‑40 (MQ=255)
                                                            aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGt  >  1:1267602/1‑69 (MQ=255)
                                                            aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGt  >  1:3969095/1‑69 (MQ=255)
                                                             aaGCGGTTGATAAATTCGAATACCGCCGTTGTTACAAGTTCTcc                          <  1:1189966/44‑1 (MQ=255)
                                                              |                                                                  
TACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGT  >  Exported/2064939‑2065067

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: