Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 609,731 A→G D40G (GAT→GGT)  entF → enterobactin non‑ribosomal peptide synthetase EntF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273609,7310AG94.7% 51.8 / ‑3.9 19D40G (GAT→GGT) entFenterobactin non‑ribosomal peptide synthetase EntF
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (9/9);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

AGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGA  >  NZ_CP009273/609644‑609815
                                                                                       |                                                                                    
aGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAggttggtt                                                                                    >  2:37653/1‑90 (MQ=255)
    cGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGcc                                                                                <  1:65621/90‑1 (MQ=255)
           tGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTg                                                                         >  1:158491/1‑90 (MQ=255)
              aTGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGcc                                                                      <  2:107302/90‑1 (MQ=255)
                 gCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCACTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCcgc                                                                   >  2:335295/1‑90 (MQ=255)
                                       ccTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAgcgc                                             <  2:6308/90‑1 (MQ=255)
                                       ccTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAgcgc                                             <  1:610121/90‑1 (MQ=255)
                                         tCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCaa                                           <  2:252984/90‑1 (MQ=255)
                                          ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                          >  1:355036/1‑90 (MQ=255)
                                          ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                          >  2:469515/1‑90 (MQ=255)
                                          ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                          >  1:550868/1‑90 (MQ=255)
                                          ccGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAg                                          >  2:189306/1‑90 (MQ=255)
                                                     gTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTg                               >  1:87082/1‑90 (MQ=255)
                                                           cATTACGTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATg                         <  1:403515/90‑1 (MQ=255)
                                                                 gTTGAGTTAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGtttt                   <  2:422184/90‑1 (MQ=255)
                                                                        tAACCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAg            <  1:658140/90‑1 (MQ=255)
                                                                         aaCCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGAc                                                   >  1:465108/1‑50 (MQ=255)
                                                                         aaCCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGAc                                                   <  2:465108/50‑1 (MQ=255)
                                                                         aaCCGGAGAGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGa           >  2:477010/1‑90 (MQ=255)
                                                                              gagaGGTTGGTTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACg      >  1:59045/1‑90 (MQ=255)
                                                                                  ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa  <  1:507361/90‑1 (MQ=255)
                                                                                  ggttggttCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGa  <  1:291766/90‑1 (MQ=255)
                                                                                       |                                                                                    
AGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATACGCTGCGGATGCGTTTTACGGAAGATAACGGCGA  >  NZ_CP009273/609644‑609815

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