Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,122,984 C→T R119H (CGT→CAT)  ftsN ← essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,122,9840CT90.7% 113.3 / 0.2 43R119H (CGT→CAT) ftsNessential cell division protein
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (19/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGT  >  NC_000913_3_bme_pgi/4122846‑4123122
                                                                                                                                          |                                                                                                                                          
tGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGCCTGTTCTTTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTgg                                                                                                                                  >  1:370757/1‑149 (MQ=255)
  ttCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGgtg                                                                                                                                >  1:298774/1‑149 (MQ=255)
        tgagtgctggGCCTCACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGTATACCAGCCTGGATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGtt                                                                                                                          <  2:827764/145‑1 (MQ=255)
               tgcGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCgg                                                                                                                   >  1:69273/1‑149 (MQ=255)
                       aCGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttca                                                                                                           >  1:717844/1‑149 (MQ=255)
                           tGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTttcacttc                                                                                                       <  1:509671/149‑1 (MQ=255)
                                cTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTcaccac                                                                                                  <  1:442444/149‑1 (MQ=255)
                                             tgccGCTGCTCCGGCGTCTGCTCATTCCACGGAACTTCAACCAGCCGCGGTGGCTGCTGTCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTc                                                                                     <  2:717844/146‑1 (MQ=255)
                                               aCGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGTTGGCGCATATCAGCCTGCATTTGTTCAGGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTGCGTTCTCTCCTGGCAGAAGTTTCtg                                                                                   >  2:249346/1‑149 (MQ=255)
                                                 gCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCTCCACCGGCAGAAGGTCCcgtg                                                                                 >  2:373745/1‑149 (MQ=255)
                                                 gCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCt                                                                                    >  2:451437/1‑146 (MQ=255)
                                                 gCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCt                                                                                    <  1:451437/146‑1 (MQ=255)
                                                 gCTGTTCCGGCGGCTGTTAATTCCACAGAACTTCAATCAGCTGCGTTGGCTGCTTGTTCATATCAGGCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCtgtg                                                                                 <  1:602879/149‑1 (MQ=255)
                                                  cTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTcacca                                                                                                   <  1:731279/130‑1 (MQ=255)
                                                  cTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTcacca                                                                                                   >  2:731279/1‑130 (MQ=255)
                                                   tGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGGGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTTTggg                                                                               >  2:597789/1‑149 (MQ=255)
                                                        cGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATGAGCCTGCATTTGTTCAAGAAGCTGACGTTTTTCTGGTCTCATTCGCGCCGGCGTTTTCACTTCACCGCCGGCGGACGGTTCTGTGGGcgcac                                                                          >  1:323684/1‑149 (MQ=255)
                                                         ggCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGcgcacg                                                                         >  2:530640/1‑149 (MQ=255)
                                                               tGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGCTGTTGTTCTGGTGTCAGTTGCGCCGGCGTTTTCACTTCACCACCGCCAGAAGGTGCTGTTGACGCActccctcc                                                                   >  2:10870/1‑149 (MQ=255)
                                                                 ttCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCCCGCACTCCCg                                                                 >  1:596710/1‑149 (MQ=255)
                                                                   cATCCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGc                                                               <  1:346993/149‑1 (MQ=255)
                                                                     ttCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTg                                                             >  1:432913/1‑149 (MQ=255)
                                                                         aCGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCCCCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGa                                                         <  1:664788/149‑1 (MQ=255)
                                                                         aCGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGa                                                         >  2:522458/1‑149 (MQ=255)
                                                                                 tCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAg                                                 <  2:432913/149‑1 (MQ=255)
                                                                                   aaCCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGCGCGCACGCACTCCCGGCTGGCGACTTTCCAGct                                               >  2:482429/1‑149 (MQ=255)
                                                                                              ttGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAg                                    <  1:126935/149‑1 (MQ=255)
                                                                                               tGGCTGCTGGCGCATATCATCATGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAgc                                   <  1:101732/149‑1 (MQ=255)
                                                                                               tGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAgc                                   >  2:288637/1‑149 (MQ=255)
                                                                                                ggCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGCTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGttttg                                  <  2:770618/149‑5 (MQ=255)
                                                                                                   tgctgGCGCATATCAGCCTGCATTTGTTCAAGAAGATGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGCCTGGCGACTTTCCAGCTCTTTATCGTAGCGCCa                               >  2:516945/1‑149 (MQ=255)
                                                                                                          gcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCaccacc                                                                                                 <  2:373564/76‑1 (MQ=255)
                                                                                                          gcATATCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCaccacc                                                                                                 >  1:373564/1‑76 (MQ=255)
                                                                                                             tatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCg                                                                                                >  2:418576/1‑74 (MQ=255)
                                                                                                             tatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCg                                                                                                <  1:418576/74‑1 (MQ=255)
                                                                                                             tatCAGCCTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTtcttct                     <  1:392312/149‑1 (MQ=255)
                                                                                                             tatCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTtcttct                     <  1:530640/149‑1 (MQ=255)
                                                                                                                    cTGCATTTGTTCAAGAAGCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCTCGCGCTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtg              <  1:373745/149‑1 (MQ=255)
                                                                                                                      gCATTTGTTCAAGAAGCTGATGTTGTTTTGGTGTCTGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGGAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTTATGTAGCGCCAGCGTTCTTCTGGTTTtggt            <  1:767686/149‑1 (MQ=255)
                                                                                                                       taTTTGTTTAAGAAGCTGATGGTGTTCTGGTGTCAGTTGCACCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTtggtg           <  1:795505/148‑1 (MQ=255)
                                                                                                                                aagaagCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACGGGGGTCGGAGTGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGt  >  2:557449/1‑149 (MQ=255)
                                                                                                                                  gaagCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGtttt               >  2:37065/1‑134 (MQ=255)
                                                                                                                                  gaagCTGATGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGtttt               <  1:37065/134‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TGTTCCGCTAATTGCTGCGCCTGACGTTGGCGCTGTAGCGTTTGCTGACGCTGTTCCGGCGTCTGTTCATTCCACGGAACTTCAACCAGCTGCGTTGGCTGCTGGCGCATATCAGCCTGCATTTGTTCAAGAAGCTGACGTTGTTCTGGTGTCAGTTGCTCCGGCGTTTTCACTTCACCACCGGCAGAAGGTTCTGTGGGCGCACGCACTCCCGGCTGGCGACTTTCCAGCTCTTTAATGTAGCGCCAGCGTTCTTCTGGTTTTGGTGGTAGTCCGT  >  NC_000913_3_bme_pgi/4122846‑4123122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑