Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	750,896	Δ1 bp	coding (1274/2448 nt)	ybgQ ←	putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	750,893	0	C	.	93.2%	308.3 / 16.7	73	coding (1277/2448 nt)	ybgQ	putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base . (39/29);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCG  >  NC_000913/750796‑750991                                                                                                  |                                                                                                   gCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGccc                                                                                                  <  1:1417478/100‑1 (MQ=255)   tgtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgc                                                                                                >  1:733234/1‑98 (MQ=255)    gtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgct                                                                                               >  2:1233738/1‑97 (MQ=255)     tgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGcccgctc                                                                                              >  2:1556243/1‑96 (MQ=255)             aCGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgct                                                                                     <  2:610953/100‑1 (MQ=255)              cGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgctg                                                                                    <  1:971902/100‑1 (MQ=255)              cGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCctgctg                                                                                    <  2:126981/100‑1 (MQ=255)               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTCATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTCCTGa                                                                                   >  1:1734033/1‑100 (MQ=255)               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCctgctg                                                                                    <  2:361609/100‑1 (MQ=255)               gATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGa                                                                                   >  2:788215/1‑100 (MQ=255)                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAg                                                                                  <  1:625765/100‑1 (MQ=255)                aTCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAg                                                                                  >  1:682102/1‑100 (MQ=255)                 tCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGt                                                                                 <  1:487359/100‑1 (MQ=255)                  cTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTa                                                                                >  2:775454/1‑100 (MQ=255)                    gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTAt                                                                               >  1:1151646/1‑100 (MQ=255)                    gagaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATc                                                                              <  1:1781460/100‑1 (MQ=255)                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaa                                                                            <  1:147750/100‑1 (MQ=255)                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaa                                                                            >  2:739308/1‑100 (MQ=255)                      gaAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaa                                                                            >  1:555420/1‑100 (MQ=255)                        aGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCaaaa                                                                          >  2:821129/1‑100 (MQ=255)                         gCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAAt                                                                         <  1:31980/100‑1 (MQ=255)                          cGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATg                                                                        >  1:1249440/1‑100 (MQ=255)                           gATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATgg                                                                       >  1:1467005/1‑100 (MQ=255)                            aTAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGc                                                                      <  1:1476399/100‑1 (MQ=255)                                gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACt                                                                  >  1:1805844/1‑100 (MQ=255)                                gcggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTCTAATTAAAACGGTCGCTCAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCt                                                                     >  1:918270/1‑97 (MQ=255)                                  ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                                <  2:504175/100‑1 (MQ=255)                                  ggcgAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTgg                                                                <  2:115436/100‑1 (MQ=255)                                    caagcgaAATCGTGCTGTTGGTGGCATCCACATGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGcc                                                              >  2:1314458/3‑100 (MQ=255)                                     gagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGcc                                                              <  1:1494526/100‑1 (MQ=255)                                     gagcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGccc                                                             >  2:448751/1‑100 (MQ=255)                                      agcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCa                                                            <  2:1094525/100‑1 (MQ=255)                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                           >  2:565261/1‑100 (MQ=255)                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                           >  2:500340/1‑100 (MQ=255)                                       gcgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAg                                                           >  1:1730216/1‑100 (MQ=255)                                        cgaAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAgg                                                          <  2:250444/100‑1 (MQ=255)                                          aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                        >  2:602214/1‑100 (MQ=255)                                          aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                        <  2:345526/100‑1 (MQ=255)                                          aaaTCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                        >  1:199961/1‑100 (MQ=255)                                            tcttccgatctggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGAc                                                      <  1:1751985/90‑1 (MQ=255)                                             tCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACa                                                     <  2:1503570/100‑1 (MQ=255)                                              cGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACAt                                                    >  2:1864025/1‑100 (MQ=255)                                               gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                                   <  2:879298/100‑1 (MQ=255)                                               gTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATc                                                   <  2:1557520/100‑1 (MQ=255)                                                    gTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAgg                                              >  1:1675459/1‑100 (MQ=255)                                                      tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                            >  2:995521/1‑100 (MQ=255)                                                      tggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGAc                                            <  2:1675469/100‑1 (MQ=255)                                                       ggtggCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACagatcggaaga                                           >  2:1751996/1‑90 (MQ=255)                                                                cgcTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                                  >  1:485395/3‑100 (MQ=255)                                                                cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                                  >  1:1496921/1‑100 (MQ=255)                                                                cACTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCaa                                  >  1:989144/1‑100 (MQ=255)                                                                  cTTGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGc                                >  2:1722403/1‑100 (MQ=255)                                                                   ttGTTTGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGcc                               <  1:1686312/100‑1 (MQ=255)                                                                      tttGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcaca                            <  1:1306460/100‑1 (MQ=255)                                                                      tttGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcaca                            >  2:507021/1‑100 (MQ=255)                                                                       ttGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCcaca                            >  2:933000/1‑100 (MQ=255)                                                                       ttGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAg                           <  1:950992/100‑1 (MQ=255)                                                                        tGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGc                          >  1:1227160/1‑100 (MQ=255)                                                                        tGCTGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGc                          <  2:592392/100‑1 (MQ=255)                                                                           tGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCata                       >  2:980052/1‑100 (MQ=255)                                                                           tGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCata                       >  2:439209/1‑100 (MQ=255)                                                                           tGTAATTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCata                       >  2:460841/1‑100 (MQ=255)                                                                              aaTTAAAACGGTAGCTTAAG‑CCTGCTCGTCCGGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt                    <  1:1314451/100‑1 (MQ=255)                                                                              aaTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAttt                    >  1:292119/1‑100 (MQ=255)                                                                               aTTAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt                   >  1:1940865/1‑100 (MQ=255)                                                                                ttAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTg                  <  1:274176/100‑1 (MQ=255)                                                                                ttAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTg                  >  1:377135/1‑100 (MQ=255)                                                                                 tAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc                 <  2:550456/100‑1 (MQ=255)                                                                                 tAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc                 >  1:378554/1‑100 (MQ=255)                                                                                 tAAAACGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGc                 >  2:1527505/1‑100 (MQ=255)                                                                                  aaaaCGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc                <  2:1467013/100‑1 (MQ=255)                                                                                   aaaCGGTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGccc               <  2:1029110/100‑1 (MQ=255)                                                                                       ggTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATa           >  1:1023130/1‑100 (MQ=255)                                                                                       ggTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATa           >  2:1750995/1‑100 (MQ=255)                                                                                        gTAGCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAc          <  2:1190016/100‑1 (MQ=255)                                                                                           gCTTAAG‑CCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCa       >  2:1056823/1‑100 (MQ=255)                                                                                                 gCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCg  >  2:1115094/1‑100 (MQ=255)                                                                                                  |                                                                                                   GCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCG  >  NC_000913/750796‑750991

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