Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,325,850 C→T R369C (CGC→TGC)  yciQ → enhancer of membrane protein expression; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,325,8500CT96.9% 213.3 / ‑4.1 64R369C (CGC→TGC) yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (29/33);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGC  >  NC_000913/1325754‑1325943
                                                                                                |                                                                                             
ggCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCa                                                                                            <  1:471395/100‑1 (MQ=255)
ggCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCa                                                                                            >  1:377333/1‑100 (MQ=255)
 gcgcGACCAGAAATAATGCGGTATTACTGCCCGCAGCGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGTGGAGATGCaa                                                                                           >  2:681712/1‑100 (MQ=255)
   gcgACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACa                                                                                         >  2:850828/1‑100 (MQ=255)
   gcgACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCCTACGGGGAGATGCAACa                                                                                         >  1:283230/1‑100 (MQ=255)
    cgACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACag                                                                                        <  2:1239044/100‑1 (MQ=255)
      atcAGAAATAATGCGGTATTACTGCTCGGAGAGCAGATATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACagag                                                                                      <  2:135267/98‑1 (MQ=255)
        cAGAAATAATGCGGTATTACTGCGCGCAGAGCAGTGATTACTGGAAGCGTTGTTTCGTAAAGGCGATAGGGTCGTTCTTACGTGGAGATGCAACAGAATc                                                                                    >  1:340997/1‑100 (MQ=255)
         aGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCt                                                                                   >  2:1109473/1‑100 (MQ=255)
          gAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGATATGCAACAGAGTCtt                                                                                  >  1:1546132/1‑100 (MQ=255)
          gAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCtt                                                                                  <  2:611305/100‑1 (MQ=255)
          gAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCtt                                                                                  >  1:1046159/1‑100 (MQ=255)
           aaaTAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTa                                                                                 <  2:516217/100‑1 (MQ=255)
             ataatGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACg                                                                               >  1:31828/1‑100 (MQ=255)
                 tGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAgg                                                                           >  2:656246/1‑100 (MQ=255)
                  gCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAggg                                                                          >  2:584260/1‑100 (MQ=255)
                  gCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAggg                                                                          <  2:1190364/100‑1 (MQ=255)
                   cGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGc                                                                         <  2:1121268/100‑1 (MQ=255)
                      tATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCAtt                                                                      >  2:1520269/1‑100 (MQ=255)
                         tACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCAttttt                                                                   <  2:81710/100‑1 (MQ=255)
                         tACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCAttttt                                                                   <  2:805514/100‑1 (MQ=255)
                                 gCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGTAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATg                                                            <  2:669656/99‑1 (MQ=255)
                                  cAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCa                                                          <  1:823119/100‑1 (MQ=255)
                                   agagCAGTTATTACTGGAAGGGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAg                                                         <  2:1408898/100‑1 (MQ=255)
                                      gCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGaaa                                                      >  1:259255/1‑100 (MQ=255)
                                        aGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAAtt                                                    >  1:507626/1‑100 (MQ=255)
                                           tattaCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTa                                                 >  1:1155262/1‑100 (MQ=255)
                                             ttaCTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATc                                               <  1:1458211/100‑1 (MQ=255)
                                                 tGGTAGCGTTATTCCGTAAAGGCGATAAGGTCGTTCTTTTGGTGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGcc                                           <  2:1281730/100‑1 (MQ=255)
                                                 tGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGcc                                           <  2:1282041/100‑1 (MQ=255)
                                                  ggAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcg                                          <  2:239666/100‑1 (MQ=255)
                                                  gcaaGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcg                                          >  1:1064137/3‑100 (MQ=255)
                                                   gAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCcgc                                         <  1:1569334/100‑1 (MQ=255)
                                                       cGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGCGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTa                                      <  1:1098468/99‑1 (MQ=255)
                                                        gTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTaaa                                    >  1:983678/1‑100 (MQ=255)
                                                        gTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTaaa                                    >  1:1081035/1‑100 (MQ=255)
                                                           aTTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAg                                 <  1:494455/100‑1 (MQ=255)
                                                             ttCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTc                               <  1:850834/100‑1 (MQ=255)
                                                                ttAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTc                            <  1:1207771/99‑1 (MQ=255)
                                                                 tAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTca                           >  2:1207764/1‑99 (MQ=255)
                                                                 tAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCg                           <  2:507623/100‑1 (MQ=255)
                                                                 tAAAGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCg                           <  1:210598/100‑1 (MQ=255)
                                                                  aaaGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGt                          <  1:537426/100‑1 (MQ=255)
                                                                   aaGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCTGCGTAAAAAGTCTTCGtt                         <  2:470272/100‑1 (MQ=255)
                                                                   aaGGCGATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGtt                         <  2:470393/100‑1 (MQ=255)
                                                                     ggCGCTAAGGTCGTTCTTCAGGGGCGACGAAATAGAGTATTACGCAGGGCATTTTTACGGATGCAGAAATTTGATCTGCCGCGGAAAAAGGCTTCTtttt                       >  2:1469707/1‑100 (MQ=255)
                                                                        gATAAGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTAcc                    <  1:1253101/100‑1 (MQ=255)
                                                                            aGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTc                >  2:1554413/1‑100 (MQ=255)
                                                                            aGGTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTc                <  2:1573875/100‑1 (MQ=255)
                                                                             ggTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCt               <  1:265868/100‑1 (MQ=255)
                                                                             ggTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCt               >  2:265867/1‑100 (MQ=255)
                                                                             ggTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCt               >  1:772389/1‑100 (MQ=255)
                                                                             gctCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCt               >  1:1067366/3‑100 (MQ=255)
                                                                              gTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTg              >  1:572065/1‑100 (MQ=255)
                                                                              gTCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTg              >  2:82685/1‑100 (MQ=255)
                                                                               tCGTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGa             <  1:238618/100‑1 (MQ=255)
                                                                                 gTTCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATa           >  1:488427/1‑100 (MQ=255)
                                                                                  ttCTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATAc          <  1:656249/100‑1 (MQ=255)
                                                                                    cTTACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACtt        >  2:313669/1‑100 (MQ=255)
                                                                                    cTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACtt        >  2:232411/1‑100 (MQ=255)
                                                                                     ttACGGGGAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACttt       >  1:335574/1‑100 (MQ=255)
                                                                                         ggggAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTg   <  1:1247737/100‑1 (MQ=255)
                                                                                          gggAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGc  >  1:937619/1‑100 (MQ=255)
                                                                                          gggAGATGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGc  <  1:32147/100‑1 (MQ=255)
                                                                                                |                                                                                             
GGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGGTCTGATACTTTTTTGC  >  NC_000913/1325754‑1325943

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