Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,065,001 A→G 5.5% 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 5.5% 490.8 / 15.2 163K414K (AAA→AAGrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (89/65);  new base G (4/5);  total (93/70)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 170 total aligned reads displayed.

AAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAG  >  Exported/2064933‑2065070
                                                                    |                                                                     
aaGAAATAGACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAg                                                                      >  1:4731039/1‑70 (MQ=255)
aaGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGc                                                                     >  1:3295182/1‑71 (MQ=255)
aaGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGc                                                                     >  1:231214/1‑71 (MQ=255)
aaGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGc                                                                     >  1:6041457/1‑71 (MQ=255)
 aGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCg                                                                    >  1:3669903/1‑71 (MQ=255)
    aaTACACCAACCGTGGCTTGCAGTTCCTTTACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 <  1:4216176/71‑1 (MQ=255)
    aaTACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 <  1:5053871/71‑1 (MQ=255)
     aTACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 >  1:1079908/1‑70 (MQ=255)
     aTACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 >  1:3788427/1‑70 (MQ=255)
     aTACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 >  1:5684646/1‑70 (MQ=255)
         acCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaa                                                            <  1:1101925/71‑1 (MQ=255)
          ccAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAATGCAACATCGGTCTGATGAAAGCGGTTGATa                                                             <  1:1029277/69‑1 (MQ=255)
          ccAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                                       >  1:2643896/1‑59 (MQ=255)
          ccAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCTGTTGATa                                                             <  1:1355618/69‑1 (MQ=255)
          ccAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                           <  1:1062831/71‑1 (MQ=255)
          ccAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                             <  1:2005012/69‑1 (MQ=255)
           cAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGa                                                               >  1:6407636/1‑66 (MQ=255)
           cAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAAt                                                          >  1:659923/1‑71 (MQ=255)
              ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                       <  1:627763/71‑1 (MQ=255)
                 tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                             <  1:3884138/62‑1 (MQ=255)
                 tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                             <  1:650810/62‑1 (MQ=255)
                   gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATa                                                   <  1:5697636/70‑1 (MQ=255)
                   gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATa                                                   <  1:928443/70‑1 (MQ=255)
                       gCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAccgcc                                              >  1:1719058/1‑71 (MQ=255)
                       gCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAccgcc                                              >  1:486590/1‑71 (MQ=255)
                        caattCCTCGACCTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATAccgtcg                                             >  1:1226268/4‑71 (MQ=16)
                        caattCCTCGACCTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATAccgtcg                                             >  1:593885/4‑71 (MQ=16)
                        caattCCTCGACCTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATAccgtcg                                             >  1:181924/4‑71 (MQ=16)
                        cAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAccgccg                                             <  1:1642067/71‑1 (MQ=255)
                          gTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTg                                           >  1:2339213/1‑71 (MQ=255)
                           ttCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTgg                                          >  1:2477734/1‑71 (MQ=255)
                            tCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGt                                                                  >  1:1931985/1‑46 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTg                                                                >  1:4948457/1‑47 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTg                                                                >  1:4288257/1‑47 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTg                                                                >  1:3479955/1‑47 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTg                                                                >  1:3773531/1‑47 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                           >  1:5387391/1‑52 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                           >  1:1324349/1‑52 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                           >  1:353584/1‑52 (MQ=255)
                             ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                         >  1:2903113/1‑70 (MQ=255)
                              cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGtt                                        <  1:4126821/70‑1 (MQ=255)
                               ttGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTa                                       >  1:5458637/1‑70 (MQ=255)
                                tGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTAc                                      >  1:2791913/1‑70 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGa                                                               >  1:5637550/1‑44 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGa                                                               >  1:691130/1‑44 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGAt                                                              >  1:6305556/1‑45 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGAt                                                              >  1:5023907/1‑45 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTgg                                          >  1:4675642/1‑65 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTgg                                          >  1:296902/1‑65 (MQ=255)
                                 gACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTg                                           >  1:849608/1‑64 (MQ=255)
                                  aCCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                                   <  1:6198136/71‑1 (MQ=255)
                                   ccTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                                  <  1:5919113/71‑3 (MQ=12)
                                    cTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGt                                  <  1:3069935/70‑1 (MQ=255)
                                    cTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGt                                  <  1:4970678/70‑1 (MQ=255)
                                    cTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                                  <  1:6131862/70‑3 (MQ=12)
                                      gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAACGCCGTGGTTACAAGTtct                               >  1:3737742/1‑71 (MQ=255)
                                       aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 >  1:2879053/1‑36 (MQ=255)
                                       aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 >  1:2594903/1‑36 (MQ=255)
                                       aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGtt                                                                 >  1:2981297/1‑36 (MQ=255)
                                       aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTct                                                       <  1:6185135/46‑2 (MQ=255)
                                       aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtc                                <  1:5334154/69‑1 (MQ=255)
                                       aTTCAGGAAGGCAACATCGGTCTGATGAAAGCGGCTGATAAATTCg                                                       <  1:4995951/46‑1 (MQ=255)
                                         ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcg                             >  1:2251713/2‑69 (MQ=255)
                                           aGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCAc                           <  1:5898393/70‑1 (MQ=255)
                                            ggAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                         <  1:1301026/71‑1 (MQ=255)
                                            ggAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                         <  1:3936255/71‑1 (MQ=255)
                                            ggAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                         <  1:5725593/71‑1 (MQ=255)
                                             gAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTa                        <  1:5372516/71‑1 (MQ=255)
                                                ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                       <  1:6151418/37‑1 (MQ=255)
                                                ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                       <  1:470093/37‑1 (MQ=255)
                                                ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                       <  1:452616/37‑1 (MQ=255)
                                                ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                       <  1:555968/37‑1 (MQ=255)
                                                ggCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                                       <  1:799211/37‑1 (MQ=255)
                                                  cAACATCGGTCTAATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCaa                   <  1:3961382/71‑1 (MQ=255)
                                                    aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                             <  1:4809046/59‑1 (MQ=255)
                                                    aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc                             <  1:1507474/59‑1 (MQ=255)
                                                    aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc                 <  1:1828410/71‑1 (MQ=255)
                                                    aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc                 >  1:3930994/1‑71 (MQ=255)
                                                      aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCtg               >  1:1432374/1‑71 (MQ=255)
                                                       tCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt                <  1:4849767/69‑1 (MQ=255)
                                                         ggTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCaa                   >  1:2589102/1‑64 (MQ=255)
                                                         ggTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCtggtg            >  1:6804802/1‑71 (MQ=255)
                                                           tCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCtggtg            >  1:663752/1‑69 (MQ=255)
                                                           tCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCtggtg            >  1:6143565/1‑69 (MQ=255)
                                                            cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGAt         >  1:3210071/1‑71 (MQ=255)
                                                            cTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGAt         >  1:276128/1‑71 (MQ=255)
                                                             tgatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCa                    <  1:199756/59‑1 (MQ=255)
                                                              gatgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCt                         <  1:4035737/53‑1 (MQ=255)
                                                               atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                              >  1:3319673/1‑47 (MQ=255)
                                                               atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCa                            >  1:5670959/1‑49 (MQ=255)
                                                               atgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAc                  >  1:4415569/1‑59 (MQ=255)
                                                                tgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc                     >  1:4294136/1‑55 (MQ=255)
                                                                tgaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTTGATCCGt     <  1:4784723/71‑1 (MQ=255)
                                                                 gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCCCCTACGCAACCTGGTGGATCCGt     >  1:5502235/1‑70 (MQ=255)
                                                                 gaAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGt     >  1:3396338/1‑70 (MQ=255)
                                                                  aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTc    <  1:1568276/70‑1 (MQ=255)
                                                                  aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTc    <  1:6107662/70‑1 (MQ=255)
                                                                  aaaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTc    <  1:2745295/70‑1 (MQ=255)
                                                                   aaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAg  >  1:5084517/1‑71 (MQ=255)
                                                                   aaGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAg  >  1:4452565/1‑71 (MQ=255)
                                                                    |                                                                     
AAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAG  >  Exported/2064933‑2065070

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: