Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,882,935 A→C Y357D (TAT→GAT)  fadD ← long‑chain‑fatty‑acid‑‑CoA ligase FadD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,882,9350AC95.0% 140.6 / ‑0.4 40Y357D (TAT→GAT) fadDlong‑chain‑fatty‑acid‑‑CoA ligase FadD
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (22/16);  total (23/17)
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.10e-01

GTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAACTCTTTATTGTTCAGCA  >  NZ_CP009273/1882803‑1883066
                                                                                                                                    |                                                                                                                                   
gtggtACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCtt                                                                                                                                <  2:355700/138‑1 (MQ=255)
         cATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGAt                                                                                                                       <  1:648584/138‑1 (MQ=255)
             atcatcatcatcCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCGGCAGATACTGt                                                                                                                  >  1:768738/1‑139 (MQ=255)
                   atcatcCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCGGCAGATACTGTCCGGt                                                                                                             <  2:768738/138‑1 (MQ=255)
                   atcatcCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGt                                                                                                             <  1:361299/138‑1 (MQ=255)
                               tttGGCTTCCGACGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCtt                                                                                                                                <  1:147077/107‑1 (MQ=255)
                               tttGGCTTCCGACGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCtt                                                                                                                                >  2:147077/1‑107 (MQ=255)
                                ttGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGt                                                                                                             <  2:299250/125‑1 (MQ=255)
                                ttGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGt                                                                                                             >  1:299250/1‑125 (MQ=255)
                                 tGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACg                                                                                              >  1:700058/1‑139 (MQ=255)
                                 tGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACg                                                                                              >  1:105266/1‑139 (MQ=255)
                                    cTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGctc                                                                                           <  1:125257/139‑1 (MQ=255)
                                                              aTGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCcag                                                                                                                            <  2:225255/80‑1 (MQ=255)
                                                              aTGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCcag                                                                                                                            >  1:225255/1‑80 (MQ=255)
                                                                    ccACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGc                                                           >  1:7257/1‑139 (MQ=255)
                                                                         aTGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCg                                                      >  1:597734/1‑139 (MQ=255)
                                                                          tGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGCCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCg                                                      >  1:92805/1‑138 (MQ=255)
                                                                                 cAATATCATATGGTTTAACGCTGCCCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCCCCGCCTGCGGAAAGATg                                              <  2:709501/139‑1 (MQ=255)
                                                                                  aaTATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATg                                              <  1:533006/138‑1 (MQ=255)
                                                                                    tatCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATg                                              >  1:657454/1‑136 (MQ=255)
                                                                                    tatCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATg                                              >  2:749355/1‑136 (MQ=255)
                                                                                    tatCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATg                                              >  2:270812/1‑136 (MQ=255)
                                                                                          atGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTg                                       >  2:361407/1‑137 (MQ=255)
                                                                                           tGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGga                                     <  2:7257/138‑1 (MQ=255)
                                                                                              gTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGt                                                                                                             <  1:581095/63‑1 (MQ=255)
                                                                                              gTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGt                                                                                                             >  2:581095/1‑63 (MQ=255)
                                                                                              gTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGaaa                                 >  1:96193/1‑139 (MQ=255)
                                                                                                tAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATc                               >  2:481699/1‑139 (MQ=255)
                                                                                                 aaCGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATc                               >  2:242461/1‑138 (MQ=255)
                                                                                                  aCGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCtt                                                                                                                                <  2:35921/40‑1 (MQ=255)
                                                                                                  aCGCTGACCAGCGGCGCACACTCGGTAAGGCCATCGCCtt                                                                                                                                >  1:35921/1‑40 (MQ=255)
                                                                                                   cGCTGACCAGCGGCGCACACTCGGTAAGCCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAg                            >  2:288200/1‑139 (MQ=255)
                                                                                                        aCCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAtcttct                       >  2:540069/1‑139 (MQ=255)
                                                                                                        aCCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAgctgct                       <  1:288200/139‑1 (MQ=255)
                                                                                                        aCCAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAgctgct                       >  1:419238/1‑139 (MQ=255)
                                                                                                         ccAGCGGCGCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGctgctg                      <  1:481699/139‑1 (MQ=255)
                                                                                                            gcggcgCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGaa                   <  2:419238/139‑1 (MQ=255)
                                                                                                             cggcgCACACTCGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAAc                  >  2:659665/1‑139 (MQ=255)
                                                                                                                        cGGTAAGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAACTCTTTATTGtt       >  1:180681/1‑139 (MQ=255)
                                                                                                                             aGGCCATCGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAACTCTTTATTGTTCAGCa  <  1:540069/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAACTCTTTATTGTTCAGCA  >  NZ_CP009273/1882803‑1883066

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: