Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,075,220	(T)5→6	coding (256/357 nt)	yggP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,075,215	1	.	T	98.7%	238.8 / ‑3.2	76	coding (261/357 nt)	yggP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (43/32);  total (44/32)
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.87e-01

CAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑TTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGA  >  NC_000913/3075140‑3075285                                                                             |                                                                        cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGcc                                                                          <  1:3155846/76‑1 (MQ=255)  aaaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑t                                                                        <  2:1765529/76‑1 (MQ=255)   aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCCCCCCCGGTTGCCCTACACAACCCGGGTC‑tt                                                                       >  2:291122/1‑76 (MQ=255)   aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATCCAAAACCCTGCCC‑tt                                                                       >  2:253193/1‑76 (MQ=255)   aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGccat                                                                        >  1:1517873/1‑74 (MQ=255)   aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tt                                                                       <  1:1330380/76‑1 (MQ=255)   aaaTATTCGTAGCTGTAGCCCGCAGAATATCCGCTTGGTACCCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tt                                                                       <  1:1089483/76‑1 (MQ=255)    aaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt                                                                      <  1:2149317/76‑1 (MQ=255)    aaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt                                                                      >  2:4601612/1‑76 (MQ=255)     atatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑t                                                                    <  1:3074023/76‑1 (MQ=255)     atatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑t                                                                    <  2:2207349/76‑1 (MQ=255)       atTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttt                                                                  <  2:3204085/76‑2 (MQ=255)       atTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttt                                                                  <  2:1742667/76‑2 (MQ=255)       atTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttt                                                                  <  1:807400/76‑2 (MQ=255)       atTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAAAGCTGGCC‑ttt‑ttt                                                                  >  1:3073227/1‑75 (MQ=255)       agtCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttt                                                                  <  2:2425257/74‑2 (MQ=255)        ttCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttta                                                                 <  1:1722563/76‑3 (MQ=255)        ttCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttta                                                                 <  1:2338172/76‑3 (MQ=255)        ttCGTAGCGGTAGCCCGCTGAATAGCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑ttta                                                                 <  2:2709743/76‑3 (MQ=255)         tCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑tttat                                                                >  1:1098292/1‑73 (MQ=255)         tCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑tttat                                                                <  1:2378336/76‑4 (MQ=255)         tCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt‑tttat                                                                >  1:1247404/1‑73 (MQ=255)                 gTAGCCCGCTGAATATACGCTTGGTAACCCCATCGCCGCTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCa                                                         >  2:2605624/1‑76 (MQ=255)                  tAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAg                                                        <  1:4669785/76‑1 (MQ=255)                  tAGCCCGCTGAATATCCGCTTGGGAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAg                                                        >  2:3642430/1‑76 (MQ=255)                   aGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGt                                                       >  2:847490/1‑76 (MQ=255)                    gCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCGATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTg                                                      <  1:2358482/76‑1 (MQ=255)                    gCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATCCAAAACGCTGGCCTTTT‑TTATTTTTGCCAGTg                                                      >  2:363832/1‑76 (MQ=255)                    gCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTg                                                      <  2:2604078/76‑1 (MQ=255)                    gCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTg                                                      >  1:2415375/1‑76 (MQ=255)                      ccGCTGAATATCCGCTTGGTACCCCCATCGCCGGTTGCAATACACAACGCTGGCCTTTT‑TTATATTTGCCATTgag                                                    <  2:4426220/76‑1 (MQ=255)                       cGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTCACCTTTT‑TTATATTTGTCTGTgagc                                                   >  1:1884433/1‑75 (MQ=255)                        gCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGAt                                                  >  2:2881853/1‑76 (MQ=255)                        gCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGAt                                                  >  1:4500801/1‑76 (MQ=255)                          tGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGAttt                                                >  1:2448885/1‑76 (MQ=255)                           gAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGAtttt                                               >  2:948792/1‑76 (MQ=255)                            aaTATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGAttttt                                              <  1:3753772/76‑1 (MQ=255)                             atatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTa                                             <  1:847490/76‑1 (MQ=255)                              tcgccGCTGGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTAc                                            <  2:702883/73‑1 (MQ=255)                               atCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTAcc                                           >  1:306368/1‑76 (MQ=255)                                tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTAccc                                          <  1:391165/76‑1 (MQ=255)                                tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTAccc                                          <  2:708602/76‑1 (MQ=255)                                tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTAcca                                          <  2:1179199/76‑2 (MQ=255)                                    cTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACtt                                      <  1:2000941/76‑1 (MQ=255)                                     ttGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTc                                     >  2:3248826/1‑76 (MQ=255)                                     ttGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTc                                     <  2:549870/76‑1 (MQ=255)                                     ttGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTc                                     >  1:4615051/1‑76 (MQ=255)                                     ctGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTc                                     <  2:3377847/75‑1 (MQ=255)                                      tGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCg                                    <  2:3471021/76‑1 (MQ=255)                                        gTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGAc                                  <  1:1123627/76‑1 (MQ=255)                                         tAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                 <  1:1240620/76‑1 (MQ=255)                                         tAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                 >  2:3942218/1‑76 (MQ=255)                                           aCCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACaan                               <  1:3182969/76‑2 (MQ=255)                                           aCCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAAt                               <  1:1757172/76‑1 (MQ=255)                                            ccccATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                              >  1:2551037/1‑76 (MQ=255)                                            ccccATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                              >  2:4836804/1‑76 (MQ=255)                                             cccATCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATaa                             >  2:3667186/1‑76 (MQ=255)                                              ccATCGCCGGATGCAATACAAAACGATGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACCTCGACAATAAc                            >  1:1837038/1‑76 (MQ=255)                                               cATCGCCGGTTGCAATACAAAACCCCGGCC‑TTTCTTATATTTGCCAGTGAAATTTTTACCCACCTCGACAATAAcc                           >  2:3355985/1‑76 (MQ=255)                                                aTCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAccc                          >  1:1753864/1‑76 (MQ=255)                                                aTCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAccc                          >  2:2707253/1‑76 (MQ=255)                                                aTCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAccc                          >  2:2684379/1‑76 (MQ=255)                                                 tCGCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                         >  2:3711726/1‑76 (MQ=255)                                                   gCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCgg                       >  1:659202/1‑76 (MQ=255)                                                   gCCGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCgg                       >  1:2422540/1‑76 (MQ=255)                                                     cGGTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGca                     >  1:2210966/1‑76 (MQ=255)                                                      ggTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGcac                    >  2:929893/1‑76 (MQ=255)                                                      ggTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGcac                    >  2:2186980/1‑76 (MQ=255)                                                       gTTGCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGcaca                   >  2:251010/1‑76 (MQ=255)                                                          gCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCGGGGGGTTTTTTCCCCCTTTCGCCATTACCCCCGGCCCTTTc                >  2:3967181/1‑76 (MQ=255)                                                          gCAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTc                <  1:1034271/76‑1 (MQ=255)                                                           cAATACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCa               <  2:2415374/76‑1 (MQ=255)                                                            aaTACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCAt              <  1:1871434/76‑1 (MQ=255)                                                            aaTACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCAt              >  2:3136946/1‑76 (MQ=255)                                                            aaTACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCAt              <  1:2578232/76‑1 (MQ=255)                                                             aTACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTCCCCACTTCGACAATAACCCCGGCACATTCATg             <  1:36214/76‑1 (MQ=255)                                                             aTACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATg             <  1:185463/76‑1 (MQ=255)                                                             aTACAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATg             >  2:3252318/1‑76 (MQ=255)                                                               aCAAAACGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGcc           >  2:910794/1‑76 (MQ=255)                                                                 aaaaCGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCg         <  2:2613039/76‑1 (MQ=255)                                                                   aaCGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGt       >  2:1724264/1‑76 (MQ=255)                                                                    aCGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTa      >  1:3900596/1‑76 (MQ=255)                                                                    aCGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTa      >  2:2093748/1‑76 (MQ=255)                                                                     cGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCTGCACATTCATGCCCTGTaa     >  1:1169582/1‑76 (MQ=255)                                                                     cGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTaa     <  1:3822486/76‑1 (MQ=255)                                                                     cGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTaa     >  1:814855/1‑76 (MQ=255)                                                                     cGCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTaa     >  1:1165251/1‑76 (MQ=255)                                                                      ggtGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt    <  2:96682/74‑1 (MQ=255)                                                                      gCTGGCCTTTT‑TTATATTTGCCAGTGCGATTTTTACCCCCTTCTACACTAACCCCGGCACATTCATGCCCGGTAAt    <  2:4408007/76‑1 (MQ=255)                                                                      gCTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt    >  1:389996/1‑76 (MQ=255)                                                                       cTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg   <  2:419654/76‑1 (MQ=255)                                                                       cTGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg   <  1:4225469/76‑1 (MQ=255)                                                                        tGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAt                  >  1:1475650/1‑60 (MQ=255)                                                                        tGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAt                  <  2:1475649/60‑1 (MQ=255)                                                                        tGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa  <  2:313627/76‑1 (MQ=255)                                                                        tGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa  >  1:1103672/1‑76 (MQ=255)                                                                        tGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa  >  1:3867964/1‑76 (MQ=255)                                                                        tGGCCTTTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGCCACATTCATGCCCGGTAATGa  >  2:1580788/1‑76 (MQ=255)                                                                             |                                                                        CAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑TTT‑TTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGA  >  NC_000913/3075140‑3075285

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