Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,908,423 C→T Q42Q (CAG→CAA)  yebO ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,908,4230CT97.9% 170.8 / NA 47Q42Q (CAG→CAA) yebOputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/21);  minor base A (0/1);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATCTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTAGAGG  >  NC_000913/1908277‑1908569
                                                                                                                                                  |                                                                                                                                                  
aGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGtt                                                                                                                                                >  2:84924/1‑151 (MQ=255)
     aataatATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTa                                                                                                                                           <  1:265903/151‑1 (MQ=255)
      ataataTCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGtt                                                                                                                                            >  1:109563/1‑149 (MQ=255)
         atatCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTcc                                                                                                                                       >  1:152887/1‑151 (MQ=255)
         atatCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTcc                                                                                                                                       >  1:681274/1‑151 (MQ=255)
               tcttCTTCAACCGATTTTTGACTCGCAACGGTCTTTTTATCTGCCTTCTCTGGTTCGTTTGCTTCGCACAAACGCCGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGaa                                                                                                                                  <  2:122364/147‑1 (MQ=255)
                 gTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTg                                                                                                                               >  1:201325/1‑151 (MQ=255)
                  tattCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTgg                                                                                                                              >  1:221742/3‑151 (MQ=255)
                   cTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGc                                                                                                                             >  1:697410/1‑151 (MQ=255)
                   cTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGc                                                                                                                             >  2:438186/1‑151 (MQ=255)
                        aCCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGAt                                                                                                                        <  2:220362/151‑1 (MQ=255)
                          cGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGATTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTg                                                                                                                      >  1:325997/1‑151 (MQ=255)
                           gATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTg                                                                                                                      >  2:273312/1‑150 (MQ=255)
                                      ctcAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAAcc                                                                                                            <  2:268800/149‑1 (MQ=255)
                                         aaCGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACcaca                                                                                                         >  1:741045/1‑149 (MQ=255)
                                           cGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACcaca                                                                                                         >  2:356568/1‑147 (MQ=255)
                                                         gCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGAGCCAGCACCGCCTCAAACAGTGCAATATGGTCGTTAGGCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATa                                                                                           <  1:451690/147‑1 (MQ=255)
                                                          cTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATaaga                                                                                        >  2:227450/1‑149 (MQ=255)
                                                              ctctGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAacc                                                                                  >  1:623048/1‑151 (MQ=255)
                                                                     tCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCg                                                                           >  2:143601/1‑151 (MQ=255)
                                                                      cGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGa                                                                          <  2:96250/151‑1 (MQ=255)
                                                                                  acaAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGa                                                              <  2:325549/151‑1 (MQ=255)
                                                                                    aaaCGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACg                                                            <  1:438135/151‑1 (MQ=255)
                                                                                         acgTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATg                                                       <  2:550566/151‑1 (MQ=255)
                                                                                         acgTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATg                                                       <  2:463194/151‑1 (MQ=255)
                                                                                           gTAACAGTGCATTTTGACGTTCTTGCTGATCCAGCAACGCCTCAAGCAGTTCACCTTGTTCGTTAGTCCGTGAACTGACACGATTGATAAAAAACCACAAGTTGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAAt                                                        >  1:476709/1‑148 (MQ=255)
                                                                                             aaCAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTgca                                                   <  1:756178/151‑1 (MQ=255)
                                                                                                aGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTtcatc                                                 >  2:372310/1‑150 (MQ=255)
                                                                                                   gCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACg                                              <  1:395961/150‑1 (MQ=255)
                                                                                                    cATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGcc                                            <  1:549444/151‑1 (MQ=255)
                                                                                                      ttttGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTg                                          >  1:711643/1‑151 (MQ=255)
                                                                                                       tttGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTg                                          >  1:48505/1‑150 (MQ=255)
                                                                                                        ttGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGa                                         <  2:92331/150‑1 (MQ=255)
                                                                                                         tGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAATACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAAt                                       >  1:619425/1‑151 (MQ=255)
                                                                                                                tttGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACaa                                 <  2:90126/150‑1 (MQ=255)
                                                                                                                  tGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACtt                              >  2:443617/1‑151 (MQ=255)
                                                                                                                  tGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACt                               <  2:269035/150‑1 (MQ=255)
                                                                                                                  tGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACt                               <  1:532572/150‑1 (MQ=255)
                                                                                                                   gCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACtt                              <  2:255367/150‑1 (MQ=255)
                                                                                                                        tCCAGCGACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCTCCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCg                            <  1:59254/147‑1 (MQ=255)
                                                                                                                        tCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTc                         >  2:62168/1‑150 (MQ=255)
                                                                                                                          cAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCAtt                      <  2:200926/151‑1 (MQ=255)
                                                                                                                                    tCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACctc            >  2:583243/1‑151 (MQ=255)
                                                                                                                                     cAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACctct           <  1:628799/151‑1 (MQ=255)
                                                                                                                                         cAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGt       >  1:276035/1‑151 (MQ=255)
                                                                                                                                          aGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTa      <  1:475637/151‑1 (MQ=255)
                                                                                                                                              cAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTAGAgg  >  2:498936/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATCTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTAGAGG  >  NC_000913/1908277‑1908569

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