Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,268,809 C→T *68* (TAG→TAA yjbS ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,268,8090CT97.2% 110.6 / NA 36*68* (TAG→TAAyjbSuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/15);  minor base A (1/0);  total (21/15)
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

ACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTT  >  NC_000913/4268682‑4268943
                                                                                                                               |                                                                                                                                      
aCCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATaa                                                                                                                             >  1:453595/1‑139 (MQ=255)
     atatTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCg                                                                                                                        >  1:978298/1‑139 (MQ=255)
     atatTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCg                                                                                                                        >  2:53981/1‑139 (MQ=255)
           aTTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGaaaa                                                                                                                  >  2:330767/1‑139 (MQ=255)
           aTTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTAACCCTGATTTTCTTAAATAAATCAATCAATGCGTGaaaa                                                                                                                  >  2:904694/1‑139 (MQ=255)
                        aaTAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                                                                                                     <  1:884733/139‑1 (MQ=255)
                             aaaGGTCCTGTAATGAGAACAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAAt                                                                                                                           <  1:113272/112‑1 (MQ=255)
                             aaaGGTCCTGTAATGAGAACAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAAt                                                                                                                           >  2:113272/1‑112 (MQ=255)
                                     tGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATAc                                                                                        >  1:401438/1‑139 (MQ=255)
                                           gagaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACttcttc                                                                                  >  1:885580/1‑139 (MQ=255)
                                            agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGttaattaa                                                                                                                                     <  2:648691/87‑1 (MQ=255)
                                            agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGttaattaa                                                                                                                                     >  1:648691/1‑87 (MQ=255)
                                              aataatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTcc                                                                               >  1:752844/1‑139 (MQ=255)
                                                 aatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGa                                                                                                                     >  1:425180/1‑98 (MQ=255)
                                                 aatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGa                                                                                                                     <  2:425180/98‑1 (MQ=255)
                                                 aatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAgg                                                                            >  1:663622/1‑139 (MQ=255)
                                                         ttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTg                                                                                                                      <  2:1122757/89‑1 (MQ=255)
                                                         ttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTg                                                                                                                      >  1:1122757/1‑89 (MQ=255)
                                                                 ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTg                                                            >  2:112875/1‑139 (MQ=255)
                                                                             gCAAAGCGTGCAAATATACATAAAGTCTTTCTGACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                                <  1:1107100/139‑1 (MQ=255)
                                                                              cAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATa                                                                                                   <  1:591946/87‑1 (MQ=255)
                                                                              cAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATa                                                                                                   >  2:591946/1‑87 (MQ=255)
                                                                                      gCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATActtc                                                                                     >  1:788090/1‑93 (MQ=255)
                                                                                      gCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATActtc                                                                                     <  2:788090/93‑1 (MQ=255)
                                                                                      gCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTa                                       <  2:885580/139‑1 (MQ=255)
                                                                                                 tAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAg                            <  2:453595/139‑1 (MQ=255)
                                                                                                   aaGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGaa                          <  1:904694/139‑1 (MQ=255)
                                                                                                    aGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGAt                                                                                          >  2:46101/1‑74 (MQ=255)
                                                                                                    aGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGAt                                                                                          <  1:46101/74‑1 (MQ=255)
                                                                                                    aGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGt                                                                           <  1:558066/89‑1 (MQ=255)
                                                                                                    aGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGt                                                                           >  2:558066/1‑89 (MQ=255)
                                                                                                      tctttcttACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTaa                                                  <  1:821707/112‑1 (MQ=255)
                                                                                                      tctttcttACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTaa                                                  >  2:821707/1‑112 (MQ=255)
                                                                                                        tttcttACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCtt                     >  2:194425/1‑139 (MQ=255)
                                                                                                        tttcttACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCtt                     <  2:978298/139‑1 (MQ=255)
                                                                                                                    gATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                                >  1:287260/1‑100 (MQ=255)
                                                                                                                    gATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                                <  2:287260/100‑1 (MQ=255)
                                                                                                                          gTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCtt   >  1:772964/1‑139 (MQ=255)
                                                                                                                           ttaattaaATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCttt  <  2:401438/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
ACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGAGATATCTTT  >  NC_000913/4268682‑4268943

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