Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,269,012 A→G S200G (AGT→GGT)  yhaC → pentapetide repeats‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,0120AG98.4% 208.7 / ‑4.7 61S200G (AGT→GGT) yhaCpentapetide repeats‑related protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (34/26);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

TCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGAC  >  NC_000913/3268877‑3269143
                                                                                                                                       |                                                                                                                                   
tctctGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGtt                                                                                                                                  <  2:727661/139‑1 (MQ=255)
 ctctGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTc                                                                                                                                 <  1:170727/139‑1 (MQ=255)
           aCGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAAc                                                                                                                       >  2:290777/1‑139 (MQ=255)
            cGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACg                                                                                                                      >  2:735379/1‑139 (MQ=255)
            cGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCT‑ATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGc                                                                                                                     >  1:283262/1‑139 (MQ=255)
                    cGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAAtc                                                                                                              <  1:290777/139‑1 (MQ=255)
                      aTTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAAtctc                                                                                                            >  2:6871/1‑139 (MQ=255)
                                 tatGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattatt                                                                                                 >  2:499802/1‑139 (MQ=255)
                                 tatGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                  <  2:1004878/138‑1 (MQ=255)
                                 tatGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                  >  1:1004878/1‑138 (MQ=255)
                                     aTTTCACAGATAAACATATGGCT‑ATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTa                                                                                            <  2:283262/139‑1 (MQ=255)
                                      tttCACAGATAAACATATGCCTAATACTATTTTAAATCCTGTTAAATTGAATTCAACAAATTTTACTAATGCCAACACGTTTCAGGTTAGTTTTGTTGGTTCAATACAAACCGCCACAATCTCCTGGGATTCTTTACTa                                                                                            <  2:695028/139‑1 (MQ=255)
                                          acaGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGtt                                                                                                                                  >  1:166837/1‑97 (MQ=255)
                                          acaGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGtt                                                                                                                                  <  2:166837/97‑1 (MQ=255)
                                          acaGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTaaaaa                                                                                        >  2:730158/1‑139 (MQ=255)
                                          acaGATAAACATATGGCTAATACTATTGTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATGATTTACTaaaaa                                                                                        >  2:640695/1‑139 (MQ=255)
                                                aaaCATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGc                                                                                  <  1:6871/139‑1 (MQ=255)
                                                aaaCATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGc                                                                                  <  1:735379/139‑1 (MQ=255)
                                                        ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGtt                                                                                                                                  >  2:844212/1‑83 (MQ=255)
                                                        ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGtt                                                                                                                                  <  1:844212/83‑1 (MQ=255)
                                                        ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                          >  1:783396/1‑139 (MQ=255)
                                                        ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                          >  2:103833/1‑139 (MQ=255)
                                                            aaTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaataa                                                                      >  1:473920/1‑139 (MQ=255)
                                                             aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGc                                                                                  >  2:693795/1‑126 (MQ=255)
                                                             aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGc                                                                                  <  1:693795/126‑1 (MQ=255)
                                                             aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAg                                                                     >  2:909529/1‑139 (MQ=255)
                                                             aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAg                                                                     >  1:471322/1‑139 (MQ=255)
                                                             aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAg                                                                     >  2:867099/1‑139 (MQ=255)
                                                              tACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGc                                                                    >  2:302349/1‑139 (MQ=255)
                                                                   ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACaa                                                                                                                >  1:774557/1‑90 (MQ=255)
                                                                   ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACaa                                                                                                                <  2:774557/90‑1 (MQ=255)
                                                                   ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTggg                                                                                                       >  2:826330/1‑99 (MQ=255)
                                                                   ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTggg                                                                                                       <  1:826330/99‑1 (MQ=255)
                                                                   ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAg                                                               >  2:544663/1‑139 (MQ=255)
                                                                    tttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGc                                                              <  2:943808/139‑1 (MQ=255)
                                                                    tttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGc                                                              <  1:544663/139‑1 (MQ=255)
                                                                       aaaTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAAACTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGat                                                           <  1:909529/139‑1 (MQ=255)
                                                                        aaTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGata                                                          >  1:486269/1‑139 (MQ=255)
                                                                        aaTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGata                                                          >  2:820096/1‑139 (MQ=255)
                                                                             ttttAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTAt                                                     >  2:600337/1‑139 (MQ=255)
                                                                             ttttAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTAt                                                     >  1:554347/1‑139 (MQ=255)
                                                                                                aaTTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGt                                  <  2:471322/139‑1 (MQ=255)
                                                                                                 aTTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGtt                                 >  1:15958/1‑139 (MQ=255)
                                                                                                   tttACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaat                                                                        <  2:396899/98‑1 (MQ=255)
                                                                                                   tttACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaat                                                                        >  1:396899/1‑98 (MQ=255)
                                                                                                      aCTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaat                                                                        >  1:94950/1‑95 (MQ=255)
                                                                                                      aCTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaat                                                                        <  2:94950/95‑1 (MQ=255)
                                                                                                        tAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGaa                          <  2:514161/139‑1 (MQ=255)
                                                                                                               aaCATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                  <  2:617350/60‑1 (MQ=255)
                                                                                                               aaCATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                  >  1:617350/1‑60 (MQ=255)
                                                                                                                     tttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTcc                                                                                                           <  1:698811/45‑1 (MQ=255)
                                                                                                                     tttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTcc                                                                                                           >  2:698811/1‑45 (MQ=255)
                                                                                                                      ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTaaaact                                                                                       <  2:645238/64‑3 (MQ=255)
                                                                                                                      ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAAt                                                                                       >  1:645238/1‑64 (MQ=255)
                                                                                                                      ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAg                                   <  1:184374/116‑1 (MQ=255)
                                                                                                                      ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAg                                   >  2:184374/1‑116 (MQ=255)
                                                                                                                             tAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGga     <  2:15958/139‑1 (MQ=255)
                                                                                                                               aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTaaa               <  2:678304/127‑1 (MQ=255)
                                                                                                                               aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTaaa               >  1:678304/1‑127 (MQ=255)
                                                                                                                               aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGgaga   <  1:730158/139‑1 (MQ=255)
                                                                                                                                ttttGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGAc  >  1:155921/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                   
TCTCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGAC  >  NC_000913/3268877‑3269143

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