Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,525,135	(A)6→5	coding (9/765 nt)	yfeN →	putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,525,130	0	A	.	95.2%	259.9 / 6.2	62	coding (4/765 nt)	yfeN	putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (31/28);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCG  >  NC_000913/2525031‑2525223                                                                                                    |                                                                                              ggCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGa                                                                                               >  2:855211/1‑100 (MQ=255)      tatTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaac                                                                                          <  1:1387840/100‑2 (MQ=255)      tatTAATGCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaac                                                                                          <  1:252965/100‑2 (MQ=255)       atTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaaca                                                                                         >  1:1194043/1‑98 (MQ=255)       atTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaaca                                                                                         >  1:929664/1‑98 (MQ=255)        ttAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaacat                                                                                        <  2:1050101/100‑4 (MQ=255)          aaTTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCt                                                                                     >  2:787650/1‑100 (MQ=255)            ttCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCttt                                                                                   >  1:1509492/1‑100 (MQ=255)              ccTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTa                                                                                 >  2:1576568/1‑100 (MQ=255)              ccTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTa                                                                                 <  2:452470/100‑1 (MQ=255)              ccTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTa                                                                                 <  2:1009802/100‑1 (MQ=255)               cTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTaa                                                                                >  1:531478/1‑100 (MQ=255)                tttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAc                                                                               >  2:454436/1‑100 (MQ=255)                 ttGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAAct                                                                              <  2:1518416/100‑1 (MQ=255)                  tGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActc                                                                             >  1:415495/1‑100 (MQ=255)                  tGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActc                                                                             >  1:203890/1‑100 (MQ=255)                   gCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActct                                                                            >  1:295761/1‑100 (MQ=255)                   gCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActct                                                                            >  2:1080104/1‑100 (MQ=255)                   gCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAActct                                                                            >  1:1520619/1‑100 (MQ=255)                     aTTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTg                                                                           <  1:623507/99‑1 (MQ=255)                      ttAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGac                                                                         <  1:356296/100‑1 (MQ=255)                       tAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGaca                                                                        >  2:452597/1‑100 (MQ=255)                       tAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGaca                                                                        >  2:174948/1‑100 (MQ=255)                         aCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTAT‑AAAAAACATCTTTTAACTCTGACACt                                                                      <  2:941972/100‑1 (MQ=255)                               tGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTcc                                                                  >  2:181267/1‑98 (MQ=255)                                gTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCctct                                                               <  1:862465/100‑1 (MQ=255)                                 tCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCta                                                              <  2:1509483/100‑1 (MQ=255)                                   gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtata                                                            >  1:277027/1‑100 (MQ=255)                                   gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtata                                                            <  1:12164/100‑1 (MQ=255)                                   gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCtata                                                            <  2:1003722/100‑1 (MQ=255)                                     ttttAGCAAACCACTTTATTAAGCAATCTGCTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTatat                                                           >  2:1558589/1‑99 (MQ=255)                                      tttAGCAAACCACTTTATTTAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATAtt                                                          >  1:612153/1‑99 (MQ=255)                                      tctAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTa                                                         >  1:492568/3‑100 (MQ=255)                                        tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                       >  2:282910/1‑100 (MQ=255)                                        tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                       <  1:1396958/100‑1 (MQ=255)                                          gCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCTTCTATATTAGCGa                                                     >  2:121652/1‑100 (MQ=255)                                           cAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCTTAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGAt                                                    <  2:1110662/100‑1 (MQ=255)                                            aaaCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATa                                                   <  1:774869/100‑1 (MQ=255)                                             aaCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATAc                                                  >  1:1006100/1‑100 (MQ=255)                                                  cTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgtt                                             >  1:479293/1‑100 (MQ=255)                                                  cTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgt                                              <  1:1547006/99‑1 (MQ=255)                                                   tttATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgttg                                            <  2:994878/100‑1 (MQ=255)                                                    ttattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgttgt                                           <  2:1260785/100‑1 (MQ=255)                                                    ttattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgttgt                                           <  1:1243598/100‑1 (MQ=255)                                                         aaGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                      <  1:257857/100‑1 (MQ=255)                                                         aaGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                      <  1:1047351/100‑1 (MQ=255)                                                         aaGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                      >  2:565993/1‑100 (MQ=255)                                                           gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCAt                                    >  2:1561016/1‑100 (MQ=255)                                                           gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCAt                                    >  1:380520/1‑100 (MQ=255)                                                            cAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATg                                   >  2:1557019/1‑100 (MQ=255)                                                            cAATCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATg                                   >  2:884246/1‑100 (MQ=255)                                                              aTCTACTAATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCa                                 >  1:1102776/1‑100 (MQ=255)                                                                ctactaATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGa                               <  2:430728/100‑1 (MQ=255)                                                                  actaATAACGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAt                             <  2:1161912/100‑1 (MQ=255)                                                                     aataaCGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTa                          >  2:923066/1‑100 (MQ=255)                                                                     aataaCGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTa                          <  2:1387897/100‑1 (MQ=255)                                                                        aaCGATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTaaa                        <  2:1032847/99‑1 (MQ=255)                                                                           gATCGGCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCg                    <  1:846879/100‑1 (MQ=255)                                                                               ggCATAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGttt                <  1:158172/100‑1 (MQ=255)                                                                                  aTAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGc             <  1:736080/100‑1 (MQ=255)                                                                                   tAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGc             >  1:1098942/1‑100 (MQ=255)                                                                                   tAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCa            <  1:174948/100‑1 (MQ=255)                                                                                   tAGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCGATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGATTTGCa            >  2:1229865/1‑100 (MQ=255)                                                                                    aGGGAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAg           <  2:1145774/100‑1 (MQ=255)                                                                                     gggAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGa          >  1:199463/1‑100 (MQ=255)                                                                                       gAATTATGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGata        >  2:1545376/1‑100 (MQ=255)                                                                                            aTGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCgg    <  2:402135/100‑1 (MQ=255)                                                                                             tGCTATG‑AAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCg  <  1:872690/100‑1 (MQ=255)                                                                                                    |                                                                                              GGCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCG  >  NC_000913/2525031‑2525223

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