Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,219,033	C→T	S17S (TCC→TCT)	ymgF →	inner membrane division septum protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,219,033	0	C	T	97.0%	220.5 / ‑3.8	66	S17S (TCC→TCT)	ymgF	inner membrane division septum protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (28/36);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCCATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGACGATAAATTTGCTTAAT  >  NC_000913/1218936‑1219126                                                                                                  |                                                                                              gtttgCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTAt                                                                                             <  1:1435995/100‑1 (MQ=255)  tttgCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTAtt                                                                                            >  1:1566699/1‑100 (MQ=255)   ttgCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTg                                                                                           <  1:300240/100‑1 (MQ=255)     gCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCa                                                                                         >  1:550201/1‑100 (MQ=255)     gCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCa                                                                                         >  1:1411908/1‑100 (MQ=255)      cTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCAt                                                                                        <  2:1231261/100‑1 (MQ=255)      cTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCCATTGCAt                                                                                        <  1:824422/100‑1 (MQ=255)       tataATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCAtt                                                                                       <  1:1134731/100‑1 (MQ=255)         taATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCCCTCTCTATATCATATTTTCTATTGCATTTa                                                                                     <  2:869493/100‑1 (MQ=255)         taATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTa                                                                                     >  1:69020/1‑100 (MQ=255)         taATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTa                                                                                     >  1:556204/1‑100 (MQ=255)          aaTTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTAt                                                                                    >  1:1011830/1‑100 (MQ=255)             tGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCt                                                                                 <  2:1542698/100‑1 (MQ=255)              gTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTg                                                                                >  2:893680/1‑100 (MQ=255)              gTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTg                                                                                >  1:1198691/1‑100 (MQ=255)               tttGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGa                                                                               >  2:981168/1‑100 (MQ=255)               tttGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGa                                                                               <  2:64530/100‑1 (MQ=255)               tttGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGa                                                                               >  2:380884/1‑100 (MQ=255)                 tGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATc                                                                             >  2:1467191/1‑100 (MQ=255)                 tGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATc                                                                             <  2:1466472/100‑1 (MQ=255)                 tGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATc                                                                             >  2:947682/1‑100 (MQ=255)                  gAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCa                                                                            >  2:1475290/1‑100 (MQ=255)                  gAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCa                                                                            <  1:981176/100‑1 (MQ=255)                  gAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCa                                                                            >  1:856086/1‑100 (MQ=255)                   aaaGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCAc                                                                           <  2:92954/100‑1 (MQ=255)                    aaGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCAcc                                                                          >  2:981686/1‑100 (MQ=255)                          tGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTg                                                                    >  2:384373/1‑100 (MQ=255)                           gACAGGGTTGCCCTTAGCACCATGAACCATAGTAATAACTTGGCTTCTTTCACTCTCTTTATCATATTTTCTATTGCATTTATGCTAATCACCCTCCTgg                                                                   >  2:796874/1‑100 (MQ=255)                           gACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTgg                                                                   <  1:1007994/100‑1 (MQ=255)                              aGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCa                                                                <  2:20092/100‑1 (MQ=255)                               ggATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCAt                                                               >  1:118683/1‑100 (MQ=255)                                 aTTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATcc                                                             <  2:821243/100‑1 (MQ=255)                                 aTTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATcc                                                             <  2:556201/100‑1 (MQ=255)                                    gCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCCATAGCATTTATGCTGATCACCCTCCTGGTCATCCTTa                                                          >  1:1488476/1‑100 (MQ=255)                                     ccATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTAt                                                         <  1:357793/100‑1 (MQ=255)                                      cATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTAtt                                                        >  2:250819/1‑100 (MQ=255)                                          agtagCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCaa                                                    <  1:1134405/100‑1 (MQ=255)                                             agCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAAc                                                 >  1:1150427/1‑100 (MQ=255)                                             agCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAAc                                                 <  1:1219785/100‑1 (MQ=255)                                               cATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAAccc                                               <  1:426243/100‑1 (MQ=255)                                               cATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAAccc                                               <  1:1211571/100‑1 (MQ=255)                                                aTGAACAATAGTAATAATCTGGATTATTTCACTCTCTATGTCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCa                                              <  1:1212079/100‑1 (MQ=255)                                                aTGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCa                                              <  1:571128/100‑1 (MQ=255)                                                 tGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAg                                             <  1:874611/100‑1 (MQ=255)                                                  gAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGt                                            >  2:1472745/1‑100 (MQ=255)                                                       aTAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCgg                                       >  2:431202/1‑100 (MQ=255)                                                        tAGTAATATTCTCGATTATTTCACTCTCCAGATCATGTTTTCTATTGCATTTATGCTGTTCCCCCGCCTGGTCATCCTTATTGCAAAACCCAGTACCggg                                      <  2:985179/100‑1 (MQ=255)                                                         aGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGc                                     <  2:69020/100‑1 (MQ=255)                                                          gTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCt                                    >  2:1114822/1‑100 (MQ=255)                                                           taataaTCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTg                                   >  1:1306171/1‑100 (MQ=255)                                                            aataatCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTgg                                  >  2:187859/1‑100 (MQ=255)                                                               aatCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGgag                               <  2:1299750/100‑1 (MQ=255)                                                                atCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGgaga                              <  2:1566689/100‑1 (MQ=255)                                                                 tCTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTCTGCTGATCACCCTCCTGGTCAGCCTTATTGCAAAACCCAGTACCGGGCTGGGAGaa                             <  2:1002374/100‑1 (MQ=255)                                                                  cTGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAg                            <  1:455536/100‑1 (MQ=255)                                                                   tGGATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGt                           <  2:1474122/100‑1 (MQ=255)                                                                     gATTATTTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGc                         <  2:850341/100‑1 (MQ=255)                                                                       ttattTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCtt                       <  1:1434804/100‑1 (MQ=255)                                                                        tattTCACTCTCTATATCATATTTTCTATTGCATTTCTGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTtg                      <  1:1104836/100‑1 (MQ=255)                                                                         attTCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTtgt                     <  2:1202539/100‑1 (MQ=255)                                                                            tCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGAc                  <  1:790922/100‑1 (MQ=255)                                                                            tCACTCTCTATATCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGAc                  <  2:163639/100‑1 (MQ=255)                                                                                      tatCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGATCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGACGATAAATTTg        >  1:1566389/1‑100 (MQ=255)                                                                                        tCATATTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGACGATAAATTTGCt      <  1:1171336/100‑1 (MQ=255)                                                                                         cATATTTTCTATTGCAAATATGCAGATCACCCTCCTGGTTATCCTTTTTTCTAAATCCCATACCGGGCTAATAGAAATGCTTTTGATGATAAAtat         >  2:1332858/1‑94 (MQ=255)                                                                                            atTTTCTATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGACGATAAATTTGCTTAAt  >  2:15069/1‑100 (MQ=255)                                                                                                  |                                                                                              GTTTGCTATAATTGTTTGAAAGTTTTGACAGGATTGCCATTAGTAGCATGAACAATAGTAATAATCTGGATTATTTCACTCTCTATATCATATTTTCCATTGCATTTATGCTGATCACCCTCCTGGTCATCCTTATTGCAAAACCCAGTACCGGGCTGGGAGAAGTGCTTGTGACGATAAATTTGCTTAAT  >  NC_000913/1218936‑1219126

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