Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 246,784 G→A A25T (GCC→ACC)  yafL → putative lipoprotein and C40 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913246,7840GA97.1% 111.3 / NA 35A25T (GCC→ACC) yafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/15);  minor base C (0/1);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGG  >  NC_000913/246647‑246927
                                                                                                                                         |                                                                                                                                               
taAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCcaccac                                                                                                                                    >  1:189195/1‑151 (MQ=255)
        ttCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACcacatcacat                                                                                                                             <  2:301543/150‑1 (MQ=255)
                  aTCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTc                                                                                                                    >  2:424390/1‑149 (MQ=255)
                  aTCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTc                                                                                                                    >  1:349276/1‑149 (MQ=255)
                       ggATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTAcg                                                                                                             >  1:349832/1‑151 (MQ=255)
                         aTACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTAcgcc                                                                                                           <  1:272038/151‑1 (MQ=255)
                                     aTTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCCCATATCTCTTTCAGCTACGCCGCCCGCCAGCg                                                                                                <  2:572729/150‑1 (MQ=255)
                                     aTTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCgg                                                                                               >  1:156152/1‑151 (MQ=255)
                                        aaaGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATg                                                                                            <  2:81813/151‑1 (MQ=255)
                                                 ccAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTCCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAAc                                                                                      <  2:189620/148‑1 (MQ=255)
                                                 ccAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAAcc                                                                                     <  2:416867/149‑1 (MQ=255)
                                                                 aTGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAAc                                                                   <  1:458822/151‑1 (MQ=255)
                                                                     ccTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGt                                                                >  1:391200/1‑150 (MQ=255)
                                                                     ccTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGt                                                                <  2:89727/150‑1 (MQ=255)
                                                                     ccTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTa                                                               <  1:554433/151‑1 (MQ=255)
                                                                         gCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCa                                                            >  1:521130/1‑150 (MQ=255)
                                                                          ccGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCaa                                                           >  2:228300/1‑150 (MQ=255)
                                                                            gTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAAc                                                         >  2:615280/1‑150 (MQ=255)
                                                                            gTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAAc                                                         >  1:119576/1‑150 (MQ=255)
                                                                             tCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACt                                                        <  2:454368/150‑1 (MQ=255)
                                                                              cGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTc                                                       >  1:160846/1‑150 (MQ=255)
                                                                                       ccTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGaaaa                                              <  2:929/150‑1 (MQ=255)
                                                                                       ccTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGaaaa                                              >  1:390470/1‑150 (MQ=255)
                                                                                       ccTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGaaa                                               <  1:444170/149‑1 (MQ=255)
                                                                                                            tGTTGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATtgtg                          <  2:516120/149‑1 (MQ=255)
                                                                                                              ttGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTgg                         <  1:317979/148‑1 (MQ=255)
                                                                                                              ttGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGa                        >  2:264730/1‑149 (MQ=255)
                                                                                                               tGATTTGTTTGCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAgg                     >  1:615527/1‑151 (MQ=255)
                                                                                                                    tgtttgCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATg                >  1:348919/1‑151 (MQ=255)
                                                                                                                      tttgCCGTTTTCTTCATTTACCAGCGCCCCCACATCCCATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGc               <  1:330253/150‑1 (MQ=255)
                                                                                                                      tttgCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGcc              >  2:261927/1‑151 (MQ=255)
                                                                                                                        tgCCGTTTTCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCCTCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGcc              >  2:77080/1‑149 (MQ=255)
                                                                                                                             ttttCTTCATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGc        >  1:227685/1‑150 (MQ=255)
                                                                                                                               ttcttcATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCaa      >  1:64607/1‑150 (MQ=255)
                                                                                                                                   tcATTTACCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAgg  <  1:128400/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                               
TAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGG  >  NC_000913/246647‑246927

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