Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,721,339 T→C W26R (TGG→CGG)  ydhJ → putative membrane fusion protein YdhJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,721,3390TC97.5% 133.4 / ‑5.8 40W26R (TGG→CGG) ydhJputative membrane fusion protein YdhJ
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/18);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGG  >  NC_000913/1721206‑1721471
                                                                                                                                     |                                                                                                                                    
tATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAatt                                                                                                                                >  2:768619/1‑140 (MQ=255)
        tttGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATg                                                                                                                        >  1:557948/1‑140 (MQ=255)
            cGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAAt                                                                                                                    <  2:361178/140‑1 (MQ=255)
            cGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAAt                                                                                                                    <  1:708774/140‑1 (MQ=255)
                  gCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGt                                                                                                              >  2:142876/1‑140 (MQ=255)
                           ttGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACAcgcg                                                                                                     <  1:768619/140‑1 (MQ=255)
                                   gCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGaaa                                                                                             >  2:133956/1‑140 (MQ=255)
                                            aTTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGa                                                                                                    <  2:694678/124‑1 (MQ=255)
                                            aTTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGa                                                                                                    >  1:694678/1‑124 (MQ=255)
                                             ttGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAg                                                                                   >  1:235717/1‑140 (MQ=255)
                                             ttGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAg                                                                                   >  1:254184/1‑140 (MQ=255)
                                                 gTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACa                                                                               >  1:670818/1‑140 (MQ=255)
                                                   ggTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAgg                                                                             <  1:133956/140‑1 (MQ=255)
                                                          aTGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTAt                                                                      <  1:43962/140‑1 (MQ=255)
                                                            gTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAgg                                                                             >  1:821363/1‑131 (MQ=255)
                                                            gTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAgg                                                                             <  2:821363/131‑1 (MQ=255)
                                                               aaTTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCAccc                                                                 >  2:69965/1‑140 (MQ=255)
                                                               aaTTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCAccc                                                                 >  2:384453/1‑140 (MQ=255)
                                                                  tAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCcac                                                              <  1:667283/140‑1 (MQ=255)
                                                                               tATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACg                                                 >  1:808313/1‑140 (MQ=255)
                                                                                 tttCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTa                                               <  2:580160/140‑1 (MQ=255)
                                                                                       aacaaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTAttgttg                                         >  1:402949/1‑140 (MQ=255)
                                                                                       aacaaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTAttgttg                                         >  1:172242/1‑140 (MQ=255)
                                                                                         caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCggg                                                    <  1:29959/127‑1 (MQ=255)
                                                                                         caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCggg                                                    >  2:29959/1‑127 (MQ=255)
                                                                                         caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAg                                       <  2:808313/140‑1 (MQ=255)
                                                                                         caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAg                                       <  2:670818/140‑1 (MQ=255)
                                                                                         caaTCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAg                                       <  2:557948/140‑1 (MQ=255)
                                                                                            tCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCggg                                                    <  1:153431/124‑1 (MQ=255)
                                                                                            tCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCggg                                                    >  2:153431/1‑124 (MQ=255)
                                                                                            tCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTg                                    >  2:813104/1‑140 (MQ=255)
                                                                                                tGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAata                                >  2:548460/1‑140 (MQ=255)
                                                                                                       gTAGTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGa                         >  2:716920/1‑140 (MQ=255)
                                                                                                          gTTGCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACaa                      >  1:85722/1‑140 (MQ=255)
                                                                                                             gCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCa                   >  2:828754/1‑140 (MQ=255)
                                                                                                             gCGGTTCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGGGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCa                   >  2:251560/1‑140 (MQ=255)
                                                                                                                  tCTTGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTgg              >  1:258412/1‑140 (MQ=255)
                                                                                                                    ttGCCGGATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGtt            <  2:235717/140‑1 (MQ=255)
                                                                                                                         ggATGGTGGCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAAcgc       <  1:548460/140‑1 (MQ=255)
                                                                                                                             ggtggCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCgggg   <  1:813104/140‑1 (MQ=255)
                                                                                                                              gtggCTGCGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCggggg  <  2:258412/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCGACGGGAAAATACGCGCAGAACAGGTTTCTATCACCCCACAGGTGTCGGGACGTATTGTTGAGCTGAATATAAAAGACAACCAGCTGGTTAACGCGGGGG  >  NC_000913/1721206‑1721471

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