Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,064,992 T→C 9.3% G411G (GGT→GGC 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,9920TC 9.3% 306.4 / 17.6 107G411G (GGT→GGCrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (61/36);  new base C (7/3);  total (68/39)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 111 total aligned reads displayed.

TTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGA  >  minE/2064922‑2065062
                                                                      |                                                                      
 ttCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGt                                                                        >  1:2054521/1‑70 (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:1577864/66‑1 (MQ=255)
                         ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                     <  1:2650053/65‑1 (MQ=255)
                         ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                     <  1:3399149/65‑1 (MQ=255)
                         ccGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATa                                                     <  1:2718564/65‑1 (MQ=255)
                           gTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:2044388/1‑53 (MQ=255)
                           gTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:1726146/1‑53 (MQ=255)
                            tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                            <  1:1216528/71‑1 (MQ=255)
                            tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                            <  1:2268039/71‑1 (MQ=255)
                            tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                            <  1:3085226/71‑1 (MQ=255)
                            tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                            <  1:1285275/71‑1 (MQ=255)
                            tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                            <  1:2272023/71‑1 (MQ=255)
                            tGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAAt                                            <  1:2571706/71‑1 (MQ=255)
                              gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                          >  1:1420719/1‑71 (MQ=255)
                              gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                          >  1:4756893/1‑71 (MQ=255)
                              gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                          >  1:3951664/1‑71 (MQ=255)
                              gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                          >  1:1575573/1‑71 (MQ=255)
                              gCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATAc                                          >  1:1426723/1‑71 (MQ=255)
                                ttGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCtga                                                                    >  1:627335/1‑43 (MQ=255)
                                ttGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCtga                                                                    >  1:1012276/1‑43 (MQ=255)
                                ttGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCtga                                                                    >  1:1399795/1‑43 (MQ=255)
                                       tCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                 >  1:2204377/1‑71 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:2682547/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:2631084/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:3265514/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:1786136/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:2003139/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:3760790/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:1433436/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:4069536/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGaaa                                                               >  1:4511602/1‑40 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                 >  1:4955241/1‑70 (MQ=255)
                                        ccTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                 >  1:425082/1‑70 (MQ=255)
                                         cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                               <  1:4709308/55‑1 (MQ=255)
                                         cTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCg                                               <  1:924103/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:2829745/1‑66 (MQ=255)
                                                 gATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:491083/1‑66 (MQ=255)
                                                   ttCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGt                                 <  1:4595197/59‑1 (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:2681223/1‑40 (MQ=255)
                                                    tCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATaaa                                                   >  1:4147094/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:5027191/1‑40 (MQ=255)
                                                    ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                    >  1:3173542/2‑69 (MQ=255)
                                                    ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                    >  1:3170143/2‑69 (MQ=255)
                                                    ccAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGTTATAAATtctcga                    >  1:2428536/2‑69 (MQ=255)
                                                    acAGGAAGGCAATATCGGCCTGCTCAAGGCCGTGGATAAATTCGAATACCGTCGCGGt                                 >  1:1878456/2‑58 (MQ=255)
                                                        gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                             <  1:2127057/58‑1 (MQ=255)
                                                        gAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                             <  1:371870/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:2254085/1‑69 (MQ=255)
                                                          aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                >  1:1291949/1‑69 (MQ=255)
                                                          aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTa                >  1:5043968/1‑69 (MQ=255)
                                                          aGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTATCCACCTa                >  1:1578098/1‑69 (MQ=255)
                                                               aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc         >  1:4898665/1‑71 (MQ=255)
                                                               aCATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAAcc         >  1:991692/1‑71 (MQ=255)
                                                                taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                         >  1:75893/2‑54 (MQ=255)
                                                                taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                         >  1:2656642/2‑54 (MQ=255)
                                                                taTCGGCCTGATGAAGGCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAAtt                         >  1:4970356/2‑54 (MQ=255)
                                                                cATCGTTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt        <  1:483579/71‑1 (MQ=255)
                                                                cATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCt        <  1:4866391/71‑1 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:3728995/1‑50 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:742241/1‑50 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:2774363/1‑50 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:2881537/1‑50 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:3117681/1‑50 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACaa                            >  1:3230415/1‑50 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                             >  1:72336/1‑49 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACa                             >  1:4060350/1‑49 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                      >  1:3295095/1‑56 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                      >  1:3423525/1‑56 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTtctc                      >  1:1387133/1‑56 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc             >  1:2513981/1‑65 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc             >  1:813768/1‑65 (MQ=255)
                                                                 aTCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGc             >  1:1390675/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)
                                                                      tCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGa  <  1:2469570/71‑1 (MQ=255)
                                                                      |                                                                      
TTTCTATCGCTAAGAAATACACCAACCGTGGCTTGCAGTTCCTTGACCTGATTCAGGAAGGCAACATCGGTCTGATGAAAGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGA  >  minE/2064922‑2065062

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: