Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,954,230	A→G	8.2%	intergenic (‑77/+41)	fatA_2 ← / ← fatB	Ferric‑anguibactin receptor FatA/Ferric‑anguibactin‑binding protein FatB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,954,230	0	A	G	8.2%	147.8 / 2.1	48	intergenic (‑77/+41)	fatA_2/fatB	Ferric‑anguibactin receptor FatA/Ferric‑anguibactin‑binding protein FatB
Reads supporting (aligned to +/- strand):  ref base A (31/13);  new base G (2/2);  total (33/15)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCGCAGGCAATGCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCGTGACGTACCAGGCCTGCGGGTCGA  >  ACB143/2954095‑2954359                                                                                                                                        |                                                                                                                                  ggCAGCGCAGGCAATGCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCa                                                                                                                                   >  6:309016/1‑136 (MQ=255) ggCAGCGCAGGCAATGCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCa                                                                                                                                   >  2:133632/1‑136 (MQ=255)           gCAATGCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACg                                                                                                                         >  8:220213/1‑136 (MQ=255)                gCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGCCACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAAGGAGAATCGCGTGCTACCGCAGTg                                                                                                                    >  4:109242/1‑136 (MQ=255)                gCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTGGTGTTTTTGCAGTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCGCGTGCTAACGCAGTg                                                                                                                    >  4:104538/1‑136 (MQ=255)                            gACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTg                                                                                                         >  3:143258/1‑135 (MQ=255)                            gACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTg                                                                                                         >  3:138485/1‑135 (MQ=255)                             aCGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGc                                                                                                        >  2:179690/1‑135 (MQ=255)                             aCGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGc                                                                                                        >  6:52102/1‑135 (MQ=255)                                  agTGTGTTCAAGGTACGGATGCACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAg                                                                                                  <  7:196057/136‑1 (MQ=255)                                  agTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAg                                                                                                  <  5:177962/136‑1 (MQ=255)                                   gtgtgtGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAg                                                                                                  >  3:128980/1‑135 (MQ=255)                                       gtgCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGt                                                                                             >  6:206137/1‑136 (MQ=255)                                         gCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGGATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGAAAAAAGGAAGTGc                                                                                           >  3:170261/1‑136 (MQ=255)                                           aaGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGcc                                                                                          >  8:255949/1‑135 (MQ=255)                                            aGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTa                                                                                        <  6:205279/136‑1 (MQ=255)                                                aCGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTcc                                                                                    >  1:56174/1‑136 (MQ=255)                                                     ttgACACGCGATGGCTCCAGAGGCACAGATGGCGGGGTTTTTTTAATCATTTAAAAAGATCAGAATTAAAGCAAATAAGAAACACGTCATAACGCAATGGCCGGCATTTACAAAAAACAAGTTCCTAGGCCTGaaa                                                                               >  2:114645/3‑134 (MQ=255)                                                     tGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGAGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGa                                                                                 >  2:234968/1‑134 (MQ=255)                                                            gcgATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCGAAAAGGAAGt                                                                                             <  3:214032/115‑1 (MQ=255)                                                            gcgATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCGAAAAGGAAGt                                                                                             >  4:214032/1‑115 (MQ=255)                                                            gcgATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCt                                                                         >  5:146679/1‑135 (MQ=255)                                                            gcgATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCt                                                                         >  3:186664/1‑135 (MQ=255)                                                               aTGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATc                                                                     >  8:275560/1‑136 (MQ=255)                                                                          gggACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCcgc                                                          <  3:146100/136‑1 (MQ=255)                                                                           ggACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCcgc                                                          >  6:35698/1‑135 (MQ=255)                                                                              cAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATc                                                      <  3:217524/136‑1 (MQ=255)                                                                                gATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCa                                                     >  5:266964/1‑135 (MQ=255)                                                                                 aTGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCa                                                     >  2:41303/1‑134 (MQ=255)                                                                                 aTGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCa                                                     >  6:73186/1‑134 (MQ=255)                                                                                   ggCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGc                                                 >  2:101007/1‑136 (MQ=255)                                                                                    gCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGcc                                                <  2:31390/136‑1 (MQ=255)                                                                                    gCGGGTTGTGATTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATGACGTGCTAACGCAGTGGCATGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGcc                                                <  4:276830/136‑1 (MQ=255)                                                                                      gggTTGTGTTTTTGCATTTAAAAAGCGAATAAGTAAGGCAAATGAGAATCACGTGCTACCGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTTATAGCCCCTACTCACGTCCTCCGCGATACAAGCCAg                                              >  8:259879/1‑136 (MQ=255)                                                                                        gTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGt                                             >  8:120426/1‑135 (MQ=255)                                                                                          tgtgTTATTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAg                                          >  1:228240/1‑136 (MQ=255)                                                                                              tttttGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCg                                      >  7:85705/1‑136 (MQ=255)                                                                                              tttttGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCg                                      >  1:27689/1‑136 (MQ=255)                                                                                              tttttGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCg                                      >  1:208939/1‑136 (MQ=255)                                                                                                          acacgCGAAGAAGTCAGGCAAATGCGCATCGCGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCg                          <  3:109242/132‑1 (MQ=255)                                                                                                          aaaaGCGAAGAAGTAAGGCAAATGAGAATCGCGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCg                          <  3:104538/136‑1 (MQ=255)                                                                                                          aaaaGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCg                          <  3:282274/136‑1 (MQ=255)                                                                                                          aaaaGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCg                          <  3:244726/136‑1 (MQ=255)                                                                                                          aaaaGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCg                          <  3:208751/136‑1 (MQ=255)                                                                                                          aaaaGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGGGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCg                          <  8:256223/136‑1 (MQ=255)                                                                                                                    aagTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCGTGACGTAc                  >  5:50294/1‑134 (MQ=255)                                                                                                                                 gagaATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCGTGACGTACCAGGCCTGCGGGTc    >  1:39546/1‑135 (MQ=255)                                                                                                                                  agaATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCGTGACGTACCAGGCCTGCGGGTcga  >  6:113206/1‑136 (MQ=255)                                                                                                                                       cACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAg                                                                                              <  2:319020/39‑1 (MQ=255)                                                                                                                                       cACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAg                                                                                              >  1:319020/1‑39 (MQ=255)                                                                                                                                        |                                                                                                                                  GGCAGCGCAGGCAATGCCAGGGCAGGGGACGAGAGTGTGTGCAAGGTACGGATGGACACGCGATGGCTCCAGAGGGACAGATGGCGGGTTGTGTTTTTGCATTTAAAAAGCGAAGAAGTAAGGCAAATGAGAATCACGTGCTAACGCAGTGGCAGGCATTTGCAAAAAGGAAGTGCCTAGTCCTGATAGCCCCTGATCACGTCATCCGCGATCAAAGCCAGTGAGGTCGGGCTGGCGGCCGTGACGTACCAGGCCTGCGGGTCGA  >  ACB143/2954095‑2954359

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