Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 466,148 G→A P389S (CCA→TCA)  ybaE ← uncharacterized protein YbaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913466,1480GA97.1% 105.5 / ‑4.3 34P389S (CCA→TCA) ybaEuncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/20);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGT‑CGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCAA  >  NC_000913/465921‑466385
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
aaTTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGT‑CTTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCCCAGCCTTCCGCTGCCAGGGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCTTCAAGGTTTGCGGTAGTTTGGCTTCATCCGGTACCTGCCTATGCGGAATAGGCCAGCCTGGCAATAATGCATGACGGGCGGTGa                                                                                                                                                                                                                        <  2:42387/251‑1 (MQ=255)
     aTCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGT‑CGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATccgg                                                                                                                                                                                                                   >  1:261771/1‑248 (MQ=255)
                   atGAGGTCTTCGTGTGCCTGTAGGGT‑CGTGTCGTCCCATTTTTTTGCGTGATGAAAAATAATTGTGAGTTGACAGCCTTCCGCTGCCAGTTTCGGCTTTAGGTGTTCTG‑CCATGGTATGAAGGTCGAGCTGGAGTTGATAAACCAGCGTCAGGGTTTTCGGTAGTTGGACTTCATCCGTTACTTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACt                                                                                                                                                                                                       <  1:301747/249‑1 (MQ=255)
                   atGAGGTCTGCGTGTGCCTGCAGGGT‑CGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTACCAATGCGGAATAGTCCAGACTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTATCCGACTTc                                                                                                                                                                                                     >  1:462991/1‑251 (MQ=255)
                    ttaggttcgctgtttccTGTAGGGTGCGT‑TCGTCCCAGTTTTTTGCTGTATGAAAAATAATTTTGAGTTCACAGCCTTCCGCTGCCAGGGTCGCCTGTAGTCGTTCGT‑CCATGGTATGAATGTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTTGCGGTGATCAGGTGTTCTCCGACt                                                                                                                                                                                                       <  2:371580/233‑1 (MQ=255)
                      aGGTCTGCGTGTGCCTGTAGGGT‑CGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTaa                                                                                                                                                                                                  <  1:332413/251‑1 (MQ=255)
                        gTCTGCGTGGGCCTGTAGGGT‑CGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTGCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACg                                                                                                                                                                                                <  1:191541/251‑1 (MQ=255)
                         tCTGCGTGCGCCTGAAGGGT‑CGTATCGTCCCAGCTTTTTGCGTTATGAAAAATAATTGTCAGATCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGAGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGCAAAGGTTTTAGGTAGTTTGAATTAATCCGGTACCTGCCAATGCGGCATAGTCCAGCCTGACCGTAATGAATGACTGGCGGTAAACAGGGTTTTTCCGaatactaaagt                                                                                                                                                                                               >  2:288187/1‑241 (MQ=255)
                                   ccTGTGGGGT‑CGTGGCGTCCCAGTTTTTGGAGTTATGAAAAATAATTTTGAGTTCACAGCCTTCCTCTGCCAGTGTCGCCTGTAGTCGTTCTG‑CCATGGTATGAAGTTCGATCGGTAGGTGATAAAGCAGCGTCAAGGTTTGCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGGAATGCATGACTGGCTGTGATCAGGTGTTCTTCGACTTCTAACGTTGGTAATaaa                                                                                                                                                                                     <  2:225187/251‑1 (MQ=255)
                                                      cccAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATg                                                                                                                                                                   <  1:142391/251‑1 (MQ=255)
                                                                           aaaTAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGCCTGGCGGTGATCAGGGTTTCTCCGa                                                                                                                                                                                                         >  1:437877/1‑192 (MQ=255)
                                                                                          cacaGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGAATGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCgg                                                                                                                               >  1:377490/1‑251 (MQ=255)
                                                                                              gctttCCGCTGCCAGGTTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGTTAATCCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGAc                                                                                                                           <  2:279353/248‑1 (MQ=255)
                                                                                                cttcaacgaaaaGTGTCGCCAGTCGGCCTTCAG‑CCATGGGAAGAAGTTCTATCGGTAGGTCATACACAACCGTCAGGGTTATCGGAAGCTTGACTTCATACGGAACCTGCCAATGCGGACTAGTCCAGACTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTTCAAACTCCAACGTATGCAATAAACCGGGTTGCGGAATAATGGAGATCCATTTTCGCAGCTGCCAGGTGGAGAGTCggggggtt                                                                                                                         >  2:277984/12‑247 (MQ=255)
                                                                                                        gcaatggtcgacggTAGGCGT‑CTGGCCATGGAAGGAAGTCCTATCGTGATGTGATAAACCGACGCCAAGGATTACGGTAGTTAGACTTCATCCGCTGCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCCCCGACTCCTAACGTTTGTAATAAACCGGATTGAGGAAGAATGGAGATCACGTTTCGCCCCTCCCAGAGGGAGAGTCGGGGACTTTTGCGCa                                                                                                                  <  2:462991/237‑1 (MQ=255)
                                                                                                                 gCCTGTAGGCGTTCGG‑CCATGGTAGTAAGTTCTATCGGTAGGTGATAACCCAGCGTCAAGGTTTTCGGTAGTGTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTTAAATa                                                                                                        <  2:416320/251‑1 (MQ=255)
                                                                                                                     gTAGGCGTTCGT‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTCCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACGGGCGGTGATCAGGTTATCTCCAACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCa                                                                                                     <  1:64281/250‑1 (MQ=255)
                                                                                                                             tCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTa                                                                                             <  2:261771/250‑1 (MQ=255)
                                                                                                                             tCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTa                                                                                             <  2:309498/250‑1 (MQ=255)
                                                                                                                                g‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCGGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAAc                                                                                          <  1:359925/250‑1 (MQ=255)
                                                                                                                                     tGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGagag                                                                                                                                   <  2:129226/205‑1 (MQ=255)
                                                                                                                                     tGGTATGAAGTTCTATCGGTAGGTGATAAACCAGAGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGagag                                                                                                                                   >  1:129226/1‑205 (MQ=255)
                                                                                                                                               ttCTATCGGTAGGAGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGAATTCATCCGGAACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCACCGACTTTTAACGTTTGTAATAAAACGGATTGATTAATAAAGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGATTTTGTGCGAAGTCCAATTGCAAAAACCTAAACTCATGCAGCCACTg                                                                            >  1:28876/1‑250 (MQ=255)
                                                                                                                                                       gTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGCTTTTGCGCAACGTCAAATAGCAAAAACCTaaa                                                                                           >  1:325538/1‑227 (MQ=255)
                                                                                                                                                       gTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCACCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTaaa                                                                                           <  2:325538/227‑1 (MQ=255)
                                                                                                                                                                   ccAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGc                                                       <  1:170013/251‑1 (MQ=255)
                                                                                                                                                                           aaGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACtttt                                                                                                                                                       <  1:184850/147‑1 (MQ=255)
                                                                                                                                                                           aaGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGCGATAAGGTTTTCTCCCACTTCTAACGTTTGTAATCAACCGGATTGATGAATAATGGAGATCACtttt                                                                                                                                                       >  2:184850/1‑147 (MQ=255)
                                                                                                                                                                              gTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGttt                                            <  1:229871/251‑1 (MQ=255)
                                                                                                                                                                                         tttGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTgg                                 >  1:69974/1‑250 (MQ=255)
                                                                                                                                                                                                              ccAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGAGTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGATCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAAc            <  2:302226/251‑1 (MQ=255)
                                                                                                                                                                                                                 aTGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGAAGCTCCTCCGGTTTGCCGGTGGTGATTTGCAAGGGATGCCCACAACTgg         >  2:384275/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCaa  <  1:226075/251‑1 (MQ=255)
                                                                                                                                                                                                                        aTAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAAATGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGGTTGCCGATGGGGGTTTGCACGGGATGCCGGCAACTGGTTCCCa   >  1:155133/1‑250 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
AATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGT‑CGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTG‑CCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCAA  >  NC_000913/465921‑466385

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