Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 844,358 C→T 25.9% G375S (GGT→AGT)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913844,3580CT25.9% 60.7 / 16.8 38G375S (GGT→AGT) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base T (6/4);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTCAACCGTTTTTGCAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGCGAATGCTCCG  >  NC_000913/844219‑844495
                                                                                                                                           |                                                                                                                                         
ccaTTCAACCGTTTTTGCAGCTCCGGATTGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATaaa                                                                                                                                  <  2:517857/149‑1 (MQ=255)
     cAACCGTTTTTGCAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCa                                                                                                                             <  1:397594/149‑1 (MQ=255)
       aCCGTTTTTGCAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATTGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAgg                                                                                                                           <  1:310692/149‑1 (MQ=255)
       aCCGTTTTTGCAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAATCAAACGCTGCGATACCCATCCTCGCCGCGCCGCGCACCGTTGCACCCATCATAAATTTCAgg                                                                                                                           >  1:531406/1‑149 (MQ=255)
       aCCGTTTTTGCAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCCTAAAGTTCAga                                                                                                                           >  2:161025/1‑148 (MQ=255)
         cGTTTTTGCAGCTCCGGATAGTTACGCTTAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCGCCCATAATAAATTTCAGGCt                                                                                                                         >  1:8189/1‑149 (MQ=255)
                 cAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCg                                                                                                                 >  2:155943/1‑149 (MQ=255)
                  aGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAACGTGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCgg                                                                                                                >  1:153479/1‑149 (MQ=255)
                        tgacaGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGaa                                                                                                          <  2:616913/145‑1 (MQ=255)
                        ggATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGaa                                                                                                          <  1:436245/149‑1 (MQ=255)
                        ggATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGaa                                                                                                          <  1:236228/149‑1 (MQ=255)
                        ggATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGaa                                                                                                          <  1:219013/149‑1 (MQ=255)
                              ttACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACga                                                                                                    >  2:103713/1‑149 (MQ=255)
                                 cGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACgaaaa                                                                                                 <  2:201616/149‑1 (MQ=255)
                                            tgGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCa                                                                                      <  1:435823/149‑1 (MQ=255)
                                                gagaggGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGAt                                                                                  <  1:354330/149‑1 (MQ=255)
                                                    gagTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATaaaa                                                                              >  1:39786/1‑149 (MQ=255)
                                                      gTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATa                                                                            <  1:746822/149‑1 (MQ=255)
                                                         aTGGTTTTCGCCAGCCAGCGGGACACACTTCCGGGAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGgc                                                                         >  2:755985/1‑149 (MQ=255)
                                                            gTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCAGGATAAAAAAGGGGCt                                                                      >  2:165356/1‑149 (MQ=255)
                                                                    cagccagcGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCAGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGACAAAATAGGCGCCTgccaccc                                                              >  2:778860/1‑149 (MQ=255)
                                                                          gcGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATTGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTgac                                                        >  2:581893/1‑147 (MQ=255)
                                                                              gAAATCATACCGGAAACAAACGCTGCGATACCAATAAACGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCACGGATCGAACAACGAAAAGAAAAACAGCACGATACAATAGGCGCTTGCCAGCCAGTGCCacac                                                    <  1:510372/149‑1 (MQ=255)
                                                                                    aTACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAgcgc                                              >  1:643406/1‑149 (MQ=255)
                                                                                            aacaaacGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGgcgc                                      <  2:99303/149‑1 (MQ=255)
                                                                                                       cGATACCAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAgg                           <  1:80940/149‑1 (MQ=255)
                                                                                                       cGATACCAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAgg                           <  1:676724/149‑1 (MQ=255)
                                                                                                       cGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAgg                           >  1:630477/1‑149 (MQ=255)
                                                                                                             cAATAATCGCCAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGaaa                     <  2:153479/149‑1 (MQ=255)
                                                                                                                  cgcGCCAGGCTGCGACCCGTTGCACCCATCATAAATTCCGGGCTGTTGCCCGTATCGAACAACGAAAAGAAAAACAGCACGATACAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGc                <  2:432168/147‑1 (MQ=255)
                                                                                                                      ccAGGCTGCGCACCGTTGCACTCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGc            >  1:586116/1‑149 (MQ=255)
                                                                                                                          gccgcgtACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGCGAAt        <  1:161025/146‑1 (MQ=255)
                                                                                                                          gCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGTCTGAAAAAGGCCAGCGAAt        <  2:531406/149‑1 (MQ=255)
                                                                                                                          gCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGACAAAGGCAAGCGaac        >  2:250001/1‑148 (MQ=255)
                                                                                                                          gCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGCGAAt        >  2:272198/1‑149 (MQ=255)
                                                                                                                             gcgcACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGCGAATGCt     >  1:563920/1‑149 (MQ=255)
                                                                                                                             gcgcACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAACACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAGACCGGCTGAAAAAGGCCAGGGACTGCt     >  1:76929/1‑149 (MQ=255)
                                                                                                                                cACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGCGAATGCTCcg  >  2:612627/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                         
CCATTCAACCGTTTTTGCAGCTCCGGATAGTTACGCTGAGTATGTGGCGAGAGAGTGATGGTTTTCGCCAGCCAGCGGGAAAACATACCGGAAACAAACGCTGCGATACCAATAATCGCCAGGCTGCGCACCGTTGCACCCATCATAAATTTCAGGCTGTTGCCCGGATCGAACAACGAAAAGAAAAACAGCACGATAAAATAGGCGCTTGCCAGCCAGTGCCACACCAGCGCAAAGGCGCGGATAAACAGGCTGAAAAAGGCCAGCGAATGCTCCG  >  NC_000913/844219‑844495

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: ‑