Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,712,670 T→C S89G (AGC→GGC)  yhjY ← conserved protein YhjY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,712,6700TC97.3% 118.0 / ‑4.0 37S89G (AGC→GGC) yhjYconserved protein YhjY
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (17/19);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

GCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTATCCCAGGTGTAACGCTCTGT  >  NC_000913/3712544‑3712802
                                                                                                                              |                                                                                                                                    
gCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTcac                                                                                                                                <  1:37630/133‑1 (MQ=255)
     gtagtaCTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCca                                                                                                                           <  1:471751/133‑1 (MQ=255)
              tCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGttt                                                                                                                  >  2:127315/1‑133 (MQ=255)
               cGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTttc                                                                                                                 >  2:65497/1‑133 (MQ=255)
                gCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTttct                                                                                                                >  1:265987/1‑133 (MQ=255)
                gCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTttct                                                                                                                >  1:274261/1‑133 (MQ=255)
                       tCGTTATACATCGAGGTGGTTGTGTCATCGTATTGTAATGCCGCGCTTGGGCCGGTTGTACCCTGTTCATATAAATGAAAATTCCAGCCCATTCGCATCCTGCCCGTCACATCCACATGGGTTTGTTCCGCCa                                                                                                         <  2:296998/133‑1 (MQ=255)
                                aTCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACtt                                                                                                >  1:33611/1‑133 (MQ=255)
                                  cGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTcc                                                                                              >  1:113733/1‑133 (MQ=255)
                                               ccATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGccc                                                                                 >  1:45309/1‑133 (MQ=255)
                                                    gtaatgtaatGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTgg                                                                            <  2:33611/133‑1 (MQ=255)
                                                            atgCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTg                                                                    >  1:302825/1‑133 (MQ=255)
                                                             tgCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGa                                                                   >  1:363130/1‑133 (MQ=255)
                                                             tgCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGa                                                                   >  2:475050/1‑133 (MQ=255)
                                                             tgCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGa                                                                   >  1:171372/1‑133 (MQ=255)
                                                             tgCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGCTTCTTCCGCCAGATTGACTTCCAGTCCCAGGGCCCTTTGGCTACGTTGa                                                                   >  2:428552/1‑133 (MQ=255)
                                                              gCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGaa                                                                  >  2:264492/1‑133 (MQ=255)
                                                              gCCGCGACCGAGCCGGTTGTAACCTGTTCATACAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGaa                                                                  <  2:196477/133‑1 (MQ=255)
                                                                   gACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGt                                                             >  1:323463/1‑133 (MQ=255)
                                                                               tGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAata                                                 <  2:45309/133‑1 (MQ=255)
                                                                                   aCCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGt                                             >  2:340344/1‑133 (MQ=255)
                                                                                    ccTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGtt                                            <  2:274261/133‑1 (MQ=255)
                                                                                        ttCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGt                                             <  1:262306/128‑1 (MQ=255)
                                                                                        ttCATATAAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGt                                             >  2:262306/1‑128 (MQ=255)
                                                                                             ataAAGGACAATTCCTGCTCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAgg                                   <  1:171390/133‑1 (MQ=255)
                                                                                             ataAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAgg                                   <  1:65497/133‑1 (MQ=255)
                                                                                             ataAAGGAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACATTGCGACAAAATATCGTTGTAATCAgg                                   >  2:325622/1‑133 (MQ=255)
                                                                                                  ggAAAATTCCAGCCCTTAATCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATATGTTGCGACAAAATATCGTTGTAATCAGGTtgat                              <  1:168384/133‑1 (MQ=255)
                                                                                                  ggAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATATGTTGGGACAAAATATCGGTGTACTCAGGTtgat                              >  2:168384/1‑133 (MQ=255)
                                                                                                  ggAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTtgat                              <  2:323463/133‑1 (MQ=255)
                                                                                                  ggAAAATTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTtgat                              <  2:113733/133‑1 (MQ=255)
                                                                                                       aTTCCAGCCCATACTCATACTGTCGGTCACGTCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCa                         <  1:230618/133‑1 (MQ=255)
                                                                                                       aTTCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCa                         <  1:428552/133‑1 (MQ=255)
                                                                                                        ttCCAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCAc                        <  2:380704/133‑1 (MQ=255)
                                                                                                           cAGCCCATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTAt                     <  2:171372/133‑1 (MQ=255)
                                                                                                              cccATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTATccc                  <  2:302825/133‑1 (MQ=255)
                                                                                                              cccATACTCATACTGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTATccc                  <  1:264492/133‑1 (MQ=255)
                                                                                                                        tactGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCATGTTGATGATCACGATCCAAGGTGTAACg        >  2:212441/1‑133 (MQ=255)
                                                                                                                        tactGCCGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTATCCCAGGTGTAACg        <  1:340344/133‑1 (MQ=255)
                                                                                                                              tgGTCACAGCCACAGAGGCTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTATCCCAGGTGTAACGCTCtgt  <  1:212441/133‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
GCGTGGTAGTACTGTCGCCAAACTCGTTATACATCGAGGTGGTTGTGCCATCGTAATGTAATGCCGCGACCGGGCCGGTTGTAACCTGTTCATATAAAGGAAAATTCCAGCCCATACTCATACTGCTGGTCACATCCACAGGGGTTTCTTCCGCCAGATTGACTTCCAGTCCCAGCGCCCTTTGGCTACTTTGAATACGTTGCGACAAAATATCGTTGTAATCAGGTTGATGATCACTATCCCAGGTGTAACGCTCTGT  >  NC_000913/3712544‑3712802

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