Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 587,625 A→T A119A (GCT→GCA ybcH ← PRK09936 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913587,6250AT96.1% 172.7 / ‑1.8 51A119A (GCT→GCAybcHPRK09936 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (21/28);  total (22/29)
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.75e-01

CGCCAGTTCAGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTC  >  NC_000913/587530‑587721
                                                                                               |                                                                                                
cgcCAGTTCAGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGctgcca                                                                                              >  1:135787/1‑100 (MQ=255)
   cAGTTCAGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCcagca                                                                                           <  2:71217/100‑1 (MQ=255)
    aGTTCAGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCcagcag                                                                                          <  1:568227/100‑1 (MQ=255)
        cAGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGa                                                                                      >  1:225513/1‑100 (MQ=255)
         aGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGAt                                                                                     <  1:165959/100‑1 (MQ=255)
           gtcgtcAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTa                                                                                   >  2:573423/1‑100 (MQ=255)
            tcgtcAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTaa                                                                                  <  1:729726/100‑1 (MQ=255)
             cgtcAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAg                                                                                 >  2:466363/1‑100 (MQ=255)
             cgtcAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAg                                                                                 >  2:466364/1‑100 (MQ=255)
                cAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATa                                                                              >  1:176193/1‑100 (MQ=255)
                 aaTTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAg                                                                             <  2:14330/100‑1 (MQ=255)
                  aTTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGc                                                                            <  2:703755/100‑1 (MQ=255)
                   tttCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCt                                                                           <  1:44642/100‑1 (MQ=255)
                     tCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCttt                                                                         >  2:468110/1‑100 (MQ=255)
                      ccGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTc                                                                        >  2:57857/1‑100 (MQ=255)
                         cgcTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAg                                                                     >  2:658820/1‑100 (MQ=255)
                         cgcTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAg                                                                     >  1:1092088/1‑100 (MQ=255)
                           cTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAgcg                                                                   >  2:371875/1‑100 (MQ=255)
                                 tACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCgg                                                             >  1:773163/1‑100 (MQ=255)
                                  aCCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGga                                                            >  1:693936/1‑100 (MQ=255)
                                  aCCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGga                                                            <  1:708064/100‑1 (MQ=255)
                                    cAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGgacg                                                          <  1:57858/100‑1 (MQ=255)
                                      gCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGgacgac                                                        >  1:610217/1‑100 (MQ=255)
                                          tCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGtt                                                    >  1:397816/1‑100 (MQ=255)
                                           ccGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGtt                                                    >  1:354824/1‑99 (MQ=255)
                                             ggCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCt                                                 <  1:841618/100‑1 (MQ=255)
                                               cGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTgg                                               <  2:971292/100‑1 (MQ=255)
                                               cGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTgg                                               <  1:762481/100‑1 (MQ=255)
                                                 ttATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTg                                             <  2:919312/100‑1 (MQ=255)
                                                 ttATGACAGGCGCGGCGCTCCATCCTCTGGATTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGAGTGGTTCCGGTg                                             <  2:240152/100‑1 (MQ=255)
                                                   aTGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCa                                           >  2:86083/1‑100 (MQ=255)
                                                    tGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCAt                                          >  1:906291/1‑100 (MQ=255)
                                                       cAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATaaa                                       <  2:178337/100‑1 (MQ=255)
                                                        aGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATaaaa                                      <  2:955214/100‑1 (MQ=255)
                                                             cggcgcTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTc                                 >  1:349156/1‑100 (MQ=255)
                                                              ggcgcTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTcc                                <  2:1091311/100‑1 (MQ=255)
                                                              ggcgcTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTcc                                <  2:322213/100‑1 (MQ=255)
                                                                cgcTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCgg                              <  2:113260/100‑1 (MQ=255)
                                                                 gcTCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCg                               <  2:135785/98‑1 (MQ=255)
                                                                   tCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATc                           <  1:23802/100‑1 (MQ=255)
                                                                   tCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATc                           <  1:675722/100‑1 (MQ=255)
                                                                   tCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATc                           >  1:768676/1‑100 (MQ=255)
                                                                   tCCATAATCTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATc                           >  1:768684/1‑100 (MQ=255)
                                                                     aaTAATCTCACTTGCTGGAGAACGGCTGCCAGCACCCGATTAAAAAAGCGTTCCAGCGTCGCGAACGTCTATTACAAATAAAAAAAAAATTCCGGATCg                          <  2:100707/98‑1 (MQ=255)
                                                                          tCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGTGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTc                    <  2:225511/100‑1 (MQ=255)
                                                                           cTGGCTTGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTc                    <  2:821358/99‑1 (MQ=255)
                                                                                 tGCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGa             <  2:454555/100‑1 (MQ=255)
                                                                                  gCTGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGAc            <  2:779282/100‑1 (MQ=255)
                                                                                    tGGAGATCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGTTGCATTAAAAATTCTGGATCGGCGTTCAGCCCGACaa          >  2:880892/1‑100 (MQ=255)
                                                                                       agaTCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATaa       <  2:176191/100‑1 (MQ=255)
                                                                                         aTCGGCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGc     <  1:705757/100‑1 (MQ=255)
                                                                                            ggCTGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTc  >  1:934037/1‑100 (MQ=255)
                                                                                               |                                                                                                
CGCCAGTTCAGGTCGTCAATTTCCGCGCTGATGTACCAGCCATCCGGCGTTATGCCAGGCGCGGCGCTCCATAATCTGGCTTGCTGGAGATCGGCAGCCAGCAGGCGATTAAGATAGCTTTCCAGCGCTGCGGACGACTGTTTCTGGTGCATAAAAAATTCCGGATCGGCGTTCAGCCCGACAATAAGCTTC  >  NC_000913/587530‑587721

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