Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,495,173 T→A A40A (GCA→GCT)  yfdY ← DUF2545 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,495,1730TA94.9% 261.9 / 4.5 78A40A (GCA→GCT) yfdYDUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (38/36);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

AGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCG  >  NC_000913/2495076‑2495269
                                                                                                 |                                                                                                
aGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGc                                                                                                >  2:1614428/1‑100 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                         >  1:1054121/1‑100 (MQ=255)
       catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                         >  1:728549/1‑100 (MQ=255)
        atcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                         >  2:1332179/1‑99 (MQ=255)
        atcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                        <  1:1283137/100‑1 (MQ=255)
        atcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                        <  2:597839/100‑1 (MQ=255)
          caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATc                                                                                      >  1:1789142/1‑100 (MQ=255)
          caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                      >  1:1049512/1‑100 (MQ=255)
          caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                      >  1:423807/1‑100 (MQ=255)
           aaGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCa                                                                                     >  1:308374/1‑100 (MQ=255)
            aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                    <  2:735975/100‑1 (MQ=255)
             gCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAcc                                                                                   >  1:1042174/1‑100 (MQ=255)
              ccACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCa                                                                                  >  2:1570708/1‑100 (MQ=255)
              ccACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCa                                                                                  <  2:728550/100‑1 (MQ=255)
                 caAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATCACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGc                                                                               <  1:1239699/100‑1 (MQ=255)
                  aaGCCGGGGCAAGATACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGcc                                                                              <  1:1324140/100‑1 (MQ=255)
                   aGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCg                                                                             >  1:767829/1‑100 (MQ=255)
                    gCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCGGCGAGGATCACCATGCCga                                                                            >  1:387258/1‑100 (MQ=255)
                      cGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaag                                                                          <  1:317752/100‑1 (MQ=255)
                        gggCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaagaa                                                                        <  1:1385967/100‑1 (MQ=255)
                            aaGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAagaaagaa                                                                    >  1:1847209/1‑100 (MQ=255)
                            aaGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaag                                                                      <  1:1178978/98‑1 (MQ=255)
                            aaGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACACGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAagaaagaa                                                                    <  1:666371/100‑1 (MQ=255)
                             aGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAaagaag                                                                   <  1:1707856/100‑1 (MQ=255)
                             aGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAaagaag                                                                   >  1:690315/1‑100 (MQ=255)
                             aGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAaagaag                                                                   >  1:542286/1‑100 (MQ=255)
                              gAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAagaaga                                                                  >  1:333205/1‑100 (MQ=255)
                               aaaCATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAagag                                                                 >  1:1745274/1‑100 (MQ=255)
                               aaaCATAAAGATCCCCGTCACTAGTTCATTACCGCCAGTGAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAagag                                                                 <  1:872785/100‑1 (MQ=255)
                                  cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAcc                                                              >  2:1808716/1‑100 (MQ=255)
                                  cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAcc                                                              <  1:1693009/100‑1 (MQ=255)
                                     aaaTATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccgcc                                                           >  1:728542/1‑100 (MQ=255)
                                       atatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGa                                                         >  1:313510/1‑100 (MQ=255)
                                         atCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATa                                                       <  1:1852512/100‑1 (MQ=255)
                                           ccccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATAcc                                                     <  2:1297692/100‑1 (MQ=255)
                                            cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    >  2:1284516/1‑100 (MQ=255)
                                             ccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCAc                                                   <  2:1553503/100‑1 (MQ=255)
                                               gTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    <  1:682664/97‑1 (MQ=255)
                                                 cACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTc                                               <  2:1509452/100‑1 (MQ=255)
                                                   cTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGCCCGCCGATACCACGTTCAg                                             >  2:1330608/1‑100 (MQ=255)
                                                     aGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTa                                           >  1:1398590/1‑100 (MQ=255)
                                                     aGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTa                                           >  2:1093842/1‑100 (MQ=255)
                                                       ttCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcc                                         >  1:997594/1‑100 (MQ=255)
                                                       ttCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcc                                         <  2:310550/100‑1 (MQ=255)
                                                        tCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACTGCCGATACCACGTTCAGTACCt                                        >  1:463275/1‑100 (MQ=255)
                                                        tCATGACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcc                                         <  2:247648/99‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  1:1807799/100‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAcc                                         <  2:1520299/98‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  1:1300963/100‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  2:1867671/100‑1 (MQ=255)
                                                         cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  2:1343168/100‑1 (MQ=255)
                                                           ttACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGAc                                     <  2:1139515/100‑1 (MQ=255)
                                                              ccGCCAGTTAACCACATGGTGAAGTAATAAATGATAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCAGACCGa                                  >  1:230522/1‑100 (MQ=255)
                                                              ccGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                  >  2:695131/1‑100 (MQ=255)
                                                                 ccAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGt                               <  1:834901/100‑1 (MQ=255)
                                                                 ccAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGt                               <  2:690316/100‑1 (MQ=255)
                                                                 ccAGTTAACCACATGATGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGt                               <  2:1915990/100‑1 (MQ=255)
                                                                  cAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTa                              >  1:1626484/1‑100 (MQ=255)
                                                                  cAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTa                              >  1:1619015/1‑100 (MQ=255)
                                                                   aGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAg                             >  2:337239/1‑100 (MQ=255)
                                                                   aGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAATAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAg                             <  1:1936131/100‑1 (MQ=255)
                                                                    gTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGc                            <  2:153239/100‑1 (MQ=255)
                                                                    gTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGc                            <  2:1232806/100‑1 (MQ=255)
                                                                       aaCCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGc                         >  2:1199317/1‑100 (MQ=255)
                                                                          cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                      >  1:1259569/1‑100 (MQ=255)
                                                                           acaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTAc                     >  1:293406/1‑100 (MQ=255)
                                                                              tGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCa                  <  2:1042175/100‑1 (MQ=255)
                                                                                  gAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACa              <  1:976375/100‑1 (MQ=255)
                                                                                    aGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAAt            >  2:969458/1‑100 (MQ=255)
                                                                                     gTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATa           >  2:1767104/1‑100 (MQ=255)
                                                                                      taataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATac          >  2:507175/1‑100 (MQ=255)
                                                                                       aataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATaca         >  2:362930/1‑100 (MQ=255)
                                                                                         taaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATAcaca       >  2:1288175/1‑100 (MQ=255)
                                                                                         taaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATAcaca       <  2:1247737/100‑1 (MQ=255)
                                                                                          aaaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAg      <  2:1515980/100‑1 (MQ=255)
                                                                                           aaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAgg     >  2:1787135/1‑100 (MQ=255)
                                                                                           aaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAgg     >  1:128583/1‑100 (MQ=255)
                                                                                            aTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAggg    <  2:1467935/100‑1 (MQ=255)
                                                                                              gagaGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCg  >  2:1644501/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCG  >  NC_000913/2495076‑2495269

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: ‑