Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,517,951 (C)6→7 coding (285/522 nt) fecI ← RNA polymerase sigma factor FecI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,9451.C90.9% 102.3 / 3.6 33coding (291/522 nt)fecIRNA polymerase sigma factor FecI
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (16/14);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

AATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑CCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGT  >  NC_000913/4517810‑4518076
                                                                                                                                        |                                                                                                                                    
aaTCTCGCTGTATGTCAGACCACCCAGTTACGAAAGCTGAAACGCTTCACGTGTTTTGCCGTTTAGTCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑cccc                                                                                                                                  <  1:235903/140‑1 (MQ=255)
aaTCTCGCTGTATGTCACACCAGCCAGTTGCTAAACCAGAAATGCTTCTCGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑cccc                                                                                                                                  <  2:333907/140‑1 (MQ=255)
       cTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑cccc‑ccctccg                                                                                                                          >  1:660982/1‑135 (MQ=255)
            tGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATa                                                                                                                      >  1:767819/1‑140 (MQ=255)
                aGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAgcg                                                                                                                  >  2:175828/1‑140 (MQ=255)
                  aCCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGcc                                                                                                                <  1:523297/140‑1 (MQ=255)
                    cATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAg                                                                                                              >  2:835414/1‑140 (MQ=255)
                     aTCCAGTTGCGAAAGCAGAAACACTACACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGc                                                                                                             <  1:292446/140‑1 (MQ=255)
                      tCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTGGCGTTCCTGAGGAGAAAGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCa                                                                                                            >  2:319861/1‑140 (MQ=255)
                      tCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCa                                                                                                            >  2:629938/1‑140 (MQ=255)
                              cGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑CCCCTCCTGCGACATAAGGGCCAGCATCTCCAGa                                                                                                    >  1:253222/1‑140 (MQ=255)
                               gAAAGCAGAAACGCTTCACGTGCTTCGCCTTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGAt                                                                                                   <  1:319861/140‑1 (MQ=255)
                                              tCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑CACC‑CCTCCGGCATAAGCACCAGCATCGCCAGGTACGCTTTTTCCAgcgc                                                                                   >  2:262213/1‑140 (MQ=255)
                                               cACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCttttt                                                                                          >  2:171414/1‑133 (MQ=255)
                                               cACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCttttt                                                                                          <  1:171414/133‑1 (MQ=255)
                                                      tttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGAAGTTGTAGGGTCTCGAGTTGGCGTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCGCCAGATACGCTTTTTCCAGCGCGTTTCgg                                                                            <  2:718011/140‑1 (MQ=255)
                                                       ttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATAGGCTTTTGCCAGCGCGTTTCGGc                                                                           >  1:378729/1‑140 (MQ=255)
                                                       ttGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGc                                                                           >  2:569866/1‑140 (MQ=255)
                                                            gTTTAGCCCGTCCAGGATGCTGTCGTGGAATTGTAGGCTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                                                      <  1:423753/140‑1 (MQ=255)
                                                            gTTTAGCCCGTCCAGCCTGCTGGCGAGGAGTTGGATGGTATCGAGTTGTCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAAc                                                                      <  2:393418/140‑1 (MQ=255)
                                                                aGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGt                                                                  <  1:175828/140‑1 (MQ=255)
                                                                  cccGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTcc                                                                >  2:450230/1‑140 (MQ=255)
                                                                         aGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCa                                                         <  1:718917/140‑1 (MQ=255)
                                                                           cATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcg                                                       <  1:569866/140‑1 (MQ=255)
                                                                             tGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCAcgcg                                                     <  2:660982/140‑1 (MQ=255)
                                                                                tGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGttt                                                  >  2:97335/1‑140 (MQ=255)
                                                                                    gaggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCg                                              >  2:18241/1‑140 (MQ=255)
                                                                                     aggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa                                             <  1:629938/140‑1 (MQ=255)
                                                                                     aggagTTGTAGGGTCTCGAGTTGGCTTTCGCGTACCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCGTCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGa                                             <  2:217889/140‑1 (MQ=255)
                                                                                       gagTTGTAGGGTCTCTAGTTGGTTTGCGCCTTCCTCATGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATa                                           <  1:268714/140‑1 (MQ=255)
                                                                                           tGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGc                                       <  1:835414/140‑1 (MQ=255)
                                                                                             tAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCGCCATCACGCGTTTGGCGATAGTGCag                                     <  2:42262/140‑1 (MQ=255)
                                                                                               gggTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCagag                                   >  1:199556/1‑140 (MQ=255)
                                                                                                ggTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCC‑CCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGagga                                 >  2:660443/1‑140 (MQ=255)
                                                                                                 gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGagga                                 >  1:78047/1‑140 (MQ=255)
                                                                                                 gTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGATTCGGCGAGACAGGTCGACCATCACACGTTTGGCGATAGTGCAGagga                                 >  2:268714/1‑140 (MQ=255)
                                                                                                            gCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGaggaa                                >  2:511921/1‑130 (MQ=255)
                                                                                                            gCTTTCGCGTTCCTCAGGTGAAGGCGCTCCCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGaggaa                                <  1:511921/130‑1 (MQ=255)
                                                                                                                                aaGGCGCTCCCCC‑CCTCCGGCATAAGCACCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGt  >  2:204393/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCAGCATGCTGTCGAGGAGTTGTAGGGTCTCGAGTTGGCTTTCGCGTTCCTCAGGTGAAGGCGCT‑CCCC‑CCTCCGGCATAAGCGCCAGCATCTCCAGATACGCTTTTTCCAGCGCGTTTCGGCGAAACAGGTCCACCATCACGCGTTTGGCGATAGTGCAGAGGAAGGAGCGAGGATCGCGGATCGTCGAGAGCGT  >  NC_000913/4517810‑4518076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑