Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	957,203	G→A	V148I (GTA→ATA)	ycaP →	UPF0702 family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	957,203	0	G	A	100.0%	256.7 / NA	75	V148I (GTA→ATA)	ycaP	UPF0702 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/37);  total (38/37)

TGGCGAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGGTACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAGTGA  >  NC_000913/957106‑957297                                                                                                   |                                                                                               tGGCGAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATa                                                                                              <  1:554933/100‑1 (MQ=255)  ggCGAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATAc                                                                                             <  1:1326539/100‑1 (MQ=255)  ggCGAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATAc                                                                                             >  1:979594/1‑100 (MQ=255)    cGNGCTNGTCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCNTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATAcn                                                                                            <  2:416929/99‑2 (MQ=255)    cGAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGt                                                                                           >  2:1132468/1‑100 (MQ=255)     gAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTc                                                                                          <  1:710391/100‑1 (MQ=255)     gAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGAGTGGCGTGGAGCAGCTGGGGCAGATACGTc                                                                                          >  2:685023/1‑100 (MQ=255)        cTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTgg                                                                                       <  1:1388709/100‑1 (MQ=255)        cTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTgg                                                                                       >  1:1401603/1‑100 (MQ=255)        cTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTgg                                                                                       >  1:1450297/1‑100 (MQ=255)         tGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTggn                                                                                      <  2:1689672/100‑2 (MQ=255)               ggTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATtct                                                                                <  1:554220/100‑1 (MQ=255)                gTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATtctc                                                                               >  2:144685/1‑100 (MQ=255)                 tCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGTAACAGCTGGGGCAGATACGTCTGGCGATTCTCg                                                                              >  1:674479/1‑100 (MQ=255)                   gAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGaa                                                                            >  1:1086808/1‑100 (MQ=255)                   gAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGaa                                                                            >  1:1675049/1‑100 (MQ=255)                   gAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGaa                                                                            >  1:259468/1‑100 (MQ=255)                   gAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCATTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGaa                                                                            >  1:965221/1‑100 (MQ=255)                    aaaCTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGaaa                                                                           <  2:843490/100‑1 (MQ=255)                         cAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAAc                                                                      >  1:710916/1‑100 (MQ=255)                         cAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAAc                                                                      <  2:1329839/100‑1 (MQ=255)                          aaTAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACg                                                                     >  2:842384/1‑100 (MQ=255)                              aCTCCAACATGACTGGAGATTG‑GCTCATGATGGATCGACGATTGCGTGGCGTGGTGCAGCTGTGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGc                                                                  >  1:1517724/1‑99 (MQ=255)                              aCTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCa                                                                 >  2:186991/1‑100 (MQ=255)                               cTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCaa                                                                >  2:94296/1‑100 (MQ=255)                                 ccAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAAt                                                              >  1:472482/1‑100 (MQ=255)                                 ccAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAAt                                                              <  2:369318/100‑1 (MQ=255)                                    aCATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAg                                                           <  1:951862/100‑1 (MQ=255)                                      aTGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAgtg                                                         <  2:1086818/100‑1 (MQ=255)                                      aTGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAgtg                                                         <  2:1002110/100‑1 (MQ=255)                                              aaTTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTAttt                                                  >  1:652938/1‑100 (MQ=255)                                                 ttGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCtt                                               <  2:1486326/100‑1 (MQ=255)                                                  tGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCttt                                              >  2:1498364/1‑100 (MQ=255)                                                  tGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCttt                                              >  2:708731/1‑100 (MQ=255)                                                  tGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCGGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCttt                                              <  1:690148/100‑1 (MQ=255)                                                   gAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTg                                             <  2:237961/100‑1 (MQ=255)                                                    ngTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGa                                            >  2:788296/2‑100 (MQ=255)                                                    aGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGCAACCAACGGGCAAATCAGTGTCTATTTCTTTGa                                            <  1:1069618/100‑1 (MQ=255)                                                    aGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGa                                            <  2:1401614/100‑1 (MQ=255)                                                    aGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGa                                            <  1:268470/100‑1 (MQ=255)                                                     gTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGaa                                           >  1:790602/1‑100 (MQ=255)                                                      ttCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAg                                          <  2:579775/100‑1 (MQ=255)                                                       tCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGa                                         >  2:181713/1‑100 (MQ=255)                                                          ttATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGa                                      >  1:306501/1‑100 (MQ=255)                                                            aTGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACa                                    >  1:1578374/1‑100 (MQ=255)                                                            aTGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACa                                    >  1:340501/1‑100 (MQ=255)                                                               gAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAgg                                 >  1:38420/1‑100 (MQ=255)                                                               gAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAgg                                 <  2:1675062/100‑1 (MQ=255)                                                                aGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGt                                <  2:652942/100‑1 (MQ=255)                                                                 gCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTg                               <  1:910054/100‑1 (MQ=255)                                                                 gCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTg                               <  2:830688/100‑1 (MQ=255)                                                                  cTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGa                              >  1:484785/1‑100 (MQ=255)                                                                    aCGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGaaa                            >  1:1782637/1‑100 (MQ=255)                                                                    aCGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGaaa                            <  1:94296/100‑1 (MQ=255)                                                                     cGATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGTGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAAc                           <  2:472484/100‑1 (MQ=255)                                                                      gATTGCGTGGCGTGGAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTATTTAAAGATGACAAGATGAAAcc                          >  2:250351/1‑100 (MQ=255)                                                                         tGCGTGGCGTGGATCAGCTGTGGTAGATACGTCTGGTGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCggg                       <  2:194420/100‑1 (MQ=255)                                                                          gcgtggcgtgGAGCAGCAGGGGCAGATACGACTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGt                      >  1:308181/1‑100 (MQ=255)                                                                          gcgtggcgcgGTGCATCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTAATTCTTTGAAGATGACAAGGTGAAACCGGGt                      >  1:308927/1‑100 (MQ=255)                                                                           cgtggcgtggAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGtt                     >  1:1686594/1‑100 (MQ=255)                                                                           cgtggcgtggAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGtt                     <  1:907399/100‑1 (MQ=255)                                                                            gtggcgtggAGCAGATGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGttt                    <  2:1450308/100‑1 (MQ=255)                                                                               gcgtggAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACt                 <  2:1152151/100‑1 (MQ=255)                                                                                  tggAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTAt              >  1:155994/1‑100 (MQ=255)                                                                                   ggAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTTAAGATGACAAGATGAAACCGGGTTTACATa               >  2:557124/1‑98 (MQ=255)                                                                                    gAGCAGCTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTAttt            <  2:1029127/100‑1 (MQ=255)                                                                                     agcagcTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTAtttt           <  2:416573/100‑1 (MQ=255)                                                                                         gcTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTAccc       >  2:199762/1‑100 (MQ=255)                                                                                         gcTGGGGCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTAccc       <  1:697394/100‑1 (MQ=255)                                                                                            gtggCAGATACGTCTGGTTATTCTTGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAGt    <  2:710920/98‑1 (MQ=255)                                                                                            ggggCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAATGTGAAACCGTGTTTACTTATTTTACCCAGt    <  1:1099553/100‑1 (MQ=255)                                                                                            ggggCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAgg    <  1:700585/100‑2 (MQ=255)                                                                                            ggggCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAGt    >  2:1523037/1‑100 (MQ=255)                                                                                            ggggCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAGt    <  2:181224/100‑1 (MQ=255)                                                                                              ggCAGATACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAGTGa  >  2:1748091/1‑100 (MQ=255)                                                                                                   |                                                                                               TGGCGAGCTGGCCTGGTCGAAACTCAATAACTCCAACATGAC‑GGAATTTGAGTTCTTTATGGAGCTACGATTGCGTGGCGTGGAGCAGCTGGGGCAGGTACGTCTGGCGATTCTCGAAACCAACGGGCAAATCAGTGTCTATTTCTTTGAAGATGACAAGGTGAAACCGGGTTTACTTATTTTACCCAGTGA  >  NC_000913/957106‑957297

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