Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,573,190	Δ1 bp	coding (1487/1542 nt)	yjiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,189	0	G	.	98.9%	394.4 / ‑2.5	90	coding (1486/1542 nt)	yjiU	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (43/45);  minor base C (1/0);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGC  >  NC_000913/4573114‑4573260                                                                            |                                                                        aaaGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAg                                                                         >  1:1714234/1‑76 (MQ=255)  aagaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAga                                                                        >  2:506908/1‑75 (MQ=255)   agaagCATGGGCCCAACTTGCCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgag                                                                       <  2:499200/76‑3 (MQ=255)   agaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgag                                                                       <  1:386438/76‑3 (MQ=255)    gaagCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAgagc                                                                      <  1:3796046/76‑4 (MQ=255)      atcATGGGCCCAACTTACCGTACCGCTTGCAGTACTGAAACCCGATGGCAACCTTTTTTATTTACCGCAA‑GAGCat                                                                   <  2:3238826/74‑1 (MQ=255)      agCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCat                                                                   >  2:4124180/1‑76 (MQ=255)       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCata                                                                  >  2:1836188/1‑76 (MQ=255)       gCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCata                                                                  <  2:3594715/76‑1 (MQ=255)        cATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCatat                                                                 >  2:621602/1‑76 (MQ=255)          tGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGa                                                               <  2:977758/76‑1 (MQ=255)           gggCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGAc                                                              >  1:2877239/1‑76 (MQ=255)               ccAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACtttt                                                          >  2:3695480/1‑76 (MQ=255)                 aaCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTc                                                        <  1:2514682/76‑1 (MQ=255)                 aaCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTc                                                        <  2:2946010/76‑1 (MQ=255)                  aCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCt                                                       <  1:3602253/76‑1 (MQ=255)                  aCTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCt                                                       >  2:3870575/1‑76 (MQ=255)                   cTTACCGGCCCTCTTGCAGTACTGAAACCCGCTTGCACCCGTTTTTATTTCCCG‑CACGAGCATATTACTTTTTcct                                                      >  2:578310/1‑74 (MQ=255)                   cTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCtt                                                      >  2:549364/1‑76 (MQ=255)                    ttACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCttt                                                     <  1:1333479/76‑1 (MQ=255)                     tACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTa                                                    >  1:1161377/1‑76 (MQ=255)                      tccGGACCGCTTGCAGTGCTGCAACCCGATGGCAGCCGTGTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTaa                                                   <  2:2444502/75‑1 (MQ=255)                      actggACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTaa                                                   >  2:3068487/4‑76 (MQ=255)                      aCCTGACCGCTTGCAGTACTGAAACACGATGGCAACCGTTTTTATTTAACGAAA‑CAGCATATGACTTTTTAGTTaa                                                   <  2:2444890/76‑1 (MQ=255)                      aCCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTaa                                                   <  2:3949545/76‑1 (MQ=255)                      aCCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTaa                                                   >  2:3591475/1‑76 (MQ=255)                       ccGGACCGCTTGGAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAAt                                                  <  2:4510930/76‑1 (MQ=255)                       ccGGACCGCTTGCAGTACTGAAATCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAAt                                                  >  2:2537963/1‑76 (MQ=255)                       ccGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGActtcttctt                                                      <  2:1351374/72‑1 (MQ=255)                       ccGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGActtcttctt                                                      >  1:1351375/1‑72 (MQ=255)                       ccGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAAt                                                  <  1:2817120/76‑1 (MQ=255)                       ccGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAAt                                                  >  1:426248/1‑76 (MQ=255)                        tggACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATc                                                 >  2:4079609/2‑76 (MQ=255)                         ggACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCa                                                <  1:2639901/76‑1 (MQ=255)                          gACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCAc                                               <  2:2877239/76‑1 (MQ=255)                          gACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCAc                                               <  1:3346460/76‑1 (MQ=255)                          gACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCAc                                               <  1:1945238/76‑1 (MQ=255)                           aCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTCACCCAAA‑GAGCATATCACTTTTTCTTTAATCatg                                              >  2:72387/1‑74 (MQ=255)                            ccGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACgg                                             <  1:720911/76‑1 (MQ=255)                              gCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGAt                                           >  2:3723802/1‑76 (MQ=255)                              gCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGAt                                           >  1:742961/1‑76 (MQ=255)                               cTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATg                                          <  1:231342/76‑1 (MQ=255)                               cTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATg                                          >  2:702958/1‑76 (MQ=255)                                ttGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGc                                         <  1:709448/76‑1 (MQ=255)                                 tGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGcc                                        >  1:1044350/1‑76 (MQ=255)                                 tGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGcc                                        <  1:4142118/76‑1 (MQ=255)                                  gCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCt                                       >  2:1581292/1‑76 (MQ=255)                                   cAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTg                                      <  2:2187384/76‑1 (MQ=255)                                     gTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGa                                    <  2:3306294/76‑1 (MQ=255)                                       aCTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                                  >  1:2386570/1‑76 (MQ=255)                                       aCTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                                  >  1:2829235/1‑76 (MQ=255)                                       aCTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                                  >  2:2521848/1‑76 (MQ=255)                                        cTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAttt                                 <  2:3813902/76‑1 (MQ=255)                                        cTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAttt                                 <  2:1851363/76‑1 (MQ=255)                                        cTGAAACCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAttt                                 >  2:3038443/1‑76 (MQ=255)                                             aCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCg                            <  1:2357784/76‑1 (MQ=255)                                             aCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCg                            >  1:2870019/1‑76 (MQ=255)                                             aCCCGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCg                            <  2:3626922/76‑1 (MQ=255)                                              cccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCgg                           <  1:917002/76‑1 (MQ=255)                                              cccGATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCgg                           >  2:4807086/1‑76 (MQ=255)                                               ccGGTGTCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGt                          <  2:3406624/76‑1 (MQ=255)                                                 gATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                        >  1:3248513/1‑76 (MQ=255)                                                 gATGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTaa                        <  1:4211027/76‑1 (MQ=255)                                                  aTGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAc                       >  2:2689189/1‑76 (MQ=255)                                                   tGGCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAAcc                      <  2:3174264/76‑1 (MQ=255)                                                    ggCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCa                     <  1:2935151/76‑1 (MQ=255)                                                     gCAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCaa                    >  1:4582150/1‑76 (MQ=255)                                                      cAACCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAAt                   <  2:3563863/76‑1 (MQ=255)                                                      cAACCGTTTTTATTTACC‑AAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAAt                   <  1:54487/76‑1 (MQ=255)                                                       aaCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATa                  <  2:2386569/76‑1 (MQ=255)                                                       aaCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATa                  >  1:2917151/1‑76 (MQ=255)                                                       aaCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATa                  >  2:2169540/1‑76 (MQ=255)                                                       aaCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATa                  >  2:3812539/1‑76 (MQ=255)                                                        aCCGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaa                 >  1:452569/1‑76 (MQ=255)                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaaa               >  1:869144/1‑76 (MQ=255)                                                          cGTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaaaa               >  1:3976524/1‑76 (MQ=255)                                                           gTTTTTATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAAt              <  1:3068487/76‑1 (MQ=255)                                                            tttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATc             <  1:3154584/76‑1 (MQ=255)                                                            tttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATc             <  1:1798282/76‑1 (MQ=255)                                                             ttttATTTACCGACA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATcc            <  1:2521848/76‑1 (MQ=255)                                                             ttttATTTACCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATcc            >  2:1808089/1‑76 (MQ=255)                                                              tttatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATCAAATCCt           >  2:2031142/1‑76 (MQ=255)                                                              tttatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCt           <  2:2459724/76‑1 (MQ=255)                                                                tatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTgg         <  2:1703495/76‑1 (MQ=255)                                                                tatttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTgg         >  2:4544259/1‑76 (MQ=255)                                                                 atttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGa        >  2:2228624/1‑76 (MQ=255)                                                                  tttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGAt       <  1:3368957/76‑1 (MQ=255)                                                                  tttaCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGAt       <  1:1444888/76‑1 (MQ=255)                                                                    taCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc     >  2:114845/1‑76 (MQ=255)                                                                    taCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc     <  1:3151768/76‑1 (MQ=255)                                                                    taCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATAc     <  1:3849568/76‑1 (MQ=255)                                                                     aCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACg    <  1:330919/76‑1 (MQ=255)                                                                     aCCGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACg    >  1:4588498/1‑76 (MQ=255)                                                                      ccGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACgg   >  2:1625594/1‑76 (MQ=255)                                                                       cGAAA‑GAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGc  >  1:2451516/1‑76 (MQ=255)                                                                            |                                                                        AAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAAAATCCTGGATACGGC  >  NC_000913/4573114‑4573260

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