Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,652 G→A 8.4% P610L (CCG→CTG)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,6520GA 8.4% 142.0 / 2.7 49P610L (CCG→CTG) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base A (2/2);  total (23/26)
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.10e-01

AATTTCTTCGTTTTCCATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATC  >  NC_000913/843504‑843800
                                                                                                                                                    |                                                                                                                                                    
aaTTTCTTCGTTTTCCATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACg                                                                                                                                                      <  6:50693/149‑1 (MQ=255)
     cttcGTTTTCCATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTGCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAat                                                                                                                                                  <  4:264763/148‑1 (MQ=255)
        cGTTTTCCATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAatgat                                                                                                                                               <  6:274692/148‑1 (MQ=255)
        cGTTTTCCATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAatgat                                                                                                                                               <  2:23531/148‑1 (MQ=255)
         gTTTTCCATTAATACCGCTACCGCAACATTCACTGCCTGATAGGCTTTATCAGCAACTTCATGGCGAACAACAACATAAATTCCCACCACCCAACCAAGCCCGCGGACAAAGTTAGCAAAGCTGGTTATCGAAGACCACGGAacgatct                                                                                                                                             >  8:65637/1‑147 (MQ=255)
              ccATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATa                                                                                                                                         <  8:14221/148‑1 (MQ=255)
                         gCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGa                                                                                                                                           >  2:325750/1‑135 (MQ=255)
                         gCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGa                                                                                                                                           <  1:325750/135‑1 (MQ=255)
                                cATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGt                                                                                                                      <  7:172176/149‑1 (MQ=255)
                                             cTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGc                                                                                                         <  3:126667/149‑1 (MQ=255)
                                             cTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCATACGCCCACGGAGc                                                                                                         <  6:138177/149‑1 (MQ=255)
                                                    gCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGa                                                                                                  <  2:166469/149‑1 (MQ=255)
                                                      tatATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATTCCGCGGACAAAGTAAGCAAAGGTGGTTATCCAAGACCACGGAATCATGTGATACGCCCCGGTATCCTGACGTACGCCCACGGAGCGAATCGACa                                                                                                >  7:138554/3‑149 (MQ=255)
                                                          tCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCg                                                                                            <  5:274903/149‑1 (MQ=255)
                                                             gCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTc                                                                                         <  6:284890/149‑1 (MQ=255)
                                                               aTCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCAC‑GAAGGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCAc                                                                                      <  8:320254/149‑1 (MQ=255)
                                                                tCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCAc                                                                                      >  8:317648/1‑149 (MQ=255)
                                                                    cATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACt                                                                                     <  5:325151/146‑1 (MQ=255)
                                                                       ggCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCg                                                                               >  3:290780/1‑149 (MQ=255)
                                                                           aTCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTc                                                                            >  8:36351/1‑148 (MQ=255)
                                                                                 atcatAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGcc                                                                     <  4:31612/149‑1 (MQ=255)
                                                                                    atAATTTGCCACCACCGAACCACTGCCGCGTACATAGTTAGCAAAGGTGGTTATCGACGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGTAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGa                                                                  <  8:138554/149‑1 (MQ=255)
                                                                                    atAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGa                                                                  <  1:264640/149‑1 (MQ=255)
                                                                                     tAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGa                                                                  >  4:163380/1‑148 (MQ=255)
                                                                                          tGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTc                                                             <  8:297510/148‑1 (MQ=255)
                                                                                            ccaccaCCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCAc                                                           >  6:50275/1‑148 (MQ=255)
                                                                                                   cGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATc                                                                                                                             >  1:320010/1‑75 (MQ=255)
                                                                                                   cGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATc                                                                                                                             <  2:320010/75‑1 (MQ=255)
                                                                                                   cGAACCAATGCCGCGGACAAAGGTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAAtctc                                                   <  7:147439/149‑1 (MQ=255)
                                                                                                    gAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGcc                                                                     <  8:318341/130‑1 (MQ=255)
                                                                                                    gAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGcc                                                                     >  7:318341/1‑130 (MQ=255)
                                                                                                     aaCCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCa                                                 <  7:317648/149‑1 (MQ=255)
                                                                                                      aCCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACAGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAg                                                >  3:134588/1‑149 (MQ=255)
                                                                                                      aCCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACAGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAg                                                >  3:139292/1‑149 (MQ=255)
                                                                                                         aaTGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAAGCTCCCgtgg                                             >  8:323275/1‑148 (MQ=255)
                                                                                                           tGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTc                                           <  4:290780/149‑1 (MQ=255)
                                                                                                                    cAAAGTTAGCAAAGGTGGTTATCGAAGACCACAGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGttt                                  <  4:139292/149‑1 (MQ=255)
                                                                                                                    cAAAGTTAGCAAAGGTGGTTATCGAAGACCACAGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGttt                                  <  4:134588/149‑1 (MQ=255)
                                                                                                                     aaaGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGtttt                                 >  7:139561/1‑149 (MQ=255)
                                                                                                                     aaaGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGtttt                                 >  7:144045/1‑149 (MQ=255)
                                                                                                                     aaaGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGtttt                                 >  7:307239/1‑149 (MQ=255)
                                                                                                                         ttAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCaaa                             >  7:229690/1‑149 (MQ=255)
                                                                                                                              aaaGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGaa                        >  2:266883/1‑149 (MQ=255)
                                                                                                                              aaaGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGaa                        >  2:183588/1‑149 (MQ=255)
                                                                                                                                     gtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                 >  6:256981/1‑149 (MQ=255)
                                                                                                                                           gAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGt           <  7:36351/149‑1 (MQ=255)
                                                                                                                                            aaGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGt           >  6:259720/1‑148 (MQ=255)
                                                                                                                                             aGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATc                                                                                                    <  6:325586/58‑1 (MQ=255)
                                                                                                                                             aGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATc                                                                                                    >  5:325586/1‑58 (MQ=255)
                                                                                                                                                 cACGGAATGATGTGCTACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataa      >  1:66923/1‑148 (MQ=255)
                                                                                                                                                    ggAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATc  <  8:229690/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATTTCTTCGTTTTCCATTAATTCCGCTACCGCATCTTTCAGTGCCTGATTGGCTTTATCAGCATCTTCATGGCGATCAACATCATAATTTGCCACCACCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATC  >  NC_000913/843504‑843800

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