Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,064,986 C→T 7.4% N409N (AAC→AAT rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,064,9860CT 7.4% 288.0 / 9.4 94N409N (AAC→AATrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base C (53/33);  new base T (4/3);  total (58/36)
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.09e-01

Only 100 of 101 total aligned reads displayed.

GTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCT  >  minE/2064918‑2065056
                                                                    |                                                                      
gTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGTCAAc                                                                        <  1:2013513/69‑1 (MQ=255)
gTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAAc                                                                        <  1:4406459/69‑1 (MQ=255)
gTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAAc                                                                        <  1:1138291/69‑1 (MQ=255)
gTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAAc                                                                        <  1:3933078/69‑1 (MQ=255)
   aTTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    <  1:3020128/70‑1 (MQ=255)
   aTTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    <  1:1839200/70‑1 (MQ=255)
    tttCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGt                                                                  >  1:1417124/1‑71 (MQ=255)
     ttCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGt                                                                  >  1:2054521/1‑70 (MQ=255)
     ttCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATCCAGGAAGGCAACATCGGt                                                                  >  1:4196186/1‑70 (MQ=255)
          tCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    >  1:1066766/1‑63 (MQ=255)
          tCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    >  1:4376074/1‑63 (MQ=255)
          tCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    >  1:3201409/1‑63 (MQ=255)
          tCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    >  1:5035159/1‑63 (MQ=255)
          tCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    >  1:976847/1‑63 (MQ=255)
          tCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCg                                                                    >  1:1156612/1‑63 (MQ=255)
                    gtACACCAACCGCGGCCTGCAATTCCTCGACCTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGtgg                                                  <  1:1573557/70‑3 (MQ=19)
                       cacCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGAt                                                <  1:2039917/70‑1 (MQ=255)
                       cacCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                               <  1:1647461/71‑1 (MQ=255)
                        acCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGAt                                                <  1:4790281/69‑1 (MQ=255)
                        acCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaa                                              <  1:4754065/71‑1 (MQ=255)
                            aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                               <  1:3518716/66‑1 (MQ=255)
                            aCCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                               <  1:1577864/66‑1 (MQ=255)
                             ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                               <  1:3399149/65‑1 (MQ=255)
                             ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                               <  1:2718564/65‑1 (MQ=255)
                             ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                               <  1:2650053/65‑1 (MQ=255)
                               gTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:1726146/1‑53 (MQ=255)
                               gTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:2044388/1‑53 (MQ=255)
                               gTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:2652513/1‑53 (MQ=255)
                                tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                      <  1:2571706/71‑1 (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:3085226/71‑1 (MQ=255)
                                tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                      <  1:1216528/71‑1 (MQ=255)
                                  gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                    >  1:1420719/1‑71 (MQ=255)
                                  gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                    >  1:3951664/1‑71 (MQ=255)
                                  gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                    >  1:1575573/1‑71 (MQ=255)
                                  gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                    >  1:4756893/1‑71 (MQ=255)
                                  gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                    >  1:1623571/1‑71 (MQ=255)
                                  gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                    >  1:1426723/1‑71 (MQ=255)
                                    ttGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCtga                                                              >  1:1399795/1‑43 (MQ=255)
                                    ttGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCtga                                                              >  1:627335/1‑43 (MQ=255)
                                    ttGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCtga                                                              >  1:1012276/1‑43 (MQ=255)
                                           tCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                           >  1:2204377/1‑71 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:1433436/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:3760790/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:4069536/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:3265514/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:1786136/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:2682547/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:2631084/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:4511602/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                         >  1:2003139/1‑40 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                           >  1:4955241/1‑70 (MQ=255)
                                            ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                           >  1:425082/1‑70 (MQ=255)
                                             cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                         <  1:924103/55‑1 (MQ=255)
                                             cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                         <  1:4709308/55‑1 (MQ=255)
                                              ttGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTa                         >  1:5009297/1‑70 (MQ=255)
                                               tGACCTGATTCAGGAAGGCAAAATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTAc                        >  1:1820771/1‑70 (MQ=255)
                                                 aCCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                     <  1:2093895/71‑1 (MQ=255)
                                                  ccTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                      <  1:1237831/53‑1 (MQ=255)
                                                  ccTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAg                     <  1:2491365/70‑1 (MQ=255)
                                                  ccTGATCCAGGAAGGCAATATCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                    <  1:1108804/71‑3 (MQ=12)
                                                   cTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGtt                   <  1:618086/71‑1 (MQ=255)
                                                    tGATCCAGGAAGGCAATATCGGCCTGATGAAGGCTGTGGACAAGTTCGAATACCGTCGCGGCTATaaat                    <  1:2569484/69‑3 (MQ=9)
                                                     gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGt                              >  1:2818843/1‑58 (MQ=255)
                                                     gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                      >  1:491083/1‑66 (MQ=255)
                                                     gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                      >  1:2829745/1‑66 (MQ=255)
                                                       ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                           <  1:1495372/59‑1 (MQ=255)
                                                       ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                           <  1:4595197/59‑1 (MQ=255)
                                                       ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc               <  1:4621225/71‑1 (MQ=255)
                                                       ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc               <  1:1625803/71‑1 (MQ=255)
                                                       ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTcc               <  1:948051/71‑1 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:5027191/1‑40 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:60001/1‑40 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:2681223/1‑40 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:393144/1‑40 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:1587890/1‑40 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:1851902/1‑40 (MQ=255)
                                                        tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                             >  1:4147094/1‑40 (MQ=255)
                                                        ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga              >  1:2428536/2‑69 (MQ=255)
                                                        ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga              >  1:3170143/2‑69 (MQ=255)
                                                        ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga              >  1:3173542/2‑69 (MQ=255)
                                                        acAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGt                           >  1:1878456/2‑58 (MQ=255)
                                                            gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                       <  1:371870/58‑1 (MQ=255)
                                                            gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                       <  1:2127057/58‑1 (MQ=255)
                                                            gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                       <  1:2026326/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:1291949/1‑69 (MQ=255)
                                                              aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa          >  1:2254085/1‑69 (MQ=255)
                                                              aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa          >  1:5043968/1‑69 (MQ=255)
                                                              aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTATCCACCTa          >  1:1578098/1‑69 (MQ=255)
                                                                   aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc   >  1:991692/1‑71 (MQ=255)
                                                                   aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc   >  1:4898665/1‑71 (MQ=255)
                                                                   aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc   <  1:3107784/71‑1 (MQ=255)
                                                                    cATCGTTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt  <  1:483579/71‑1 (MQ=255)
                                                                    cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt  <  1:4866391/71‑1 (MQ=255)
                                                                    cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt  <  1:483580/71‑1 (MQ=255)
                                                                    |                                                                      
GTTATTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCT  >  minE/2064918‑2065056

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: