Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,888,711 A→G V143A (GTC→GCC)  ygcQ ← electron transfer flavoprotein subunit alpha/FixB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,888,7110AG98.0% 152.3 / ‑3.8 49V143A (GTC→GCC) ygcQelectron transfer flavoprotein subunit alpha/FixB family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (28/20);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

CAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGCG  >  NZ_CP009273/2888623‑2888792
                                                                                        |                                                                                 
caGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGAc                                                                                  <  2:697076/90‑1 (MQ=255)
 aGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCa                                                                                 <  2:1034383/90‑1 (MQ=255)
  gCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCaa                                                                                >  2:432671/1‑90 (MQ=255)
  gCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCaa                                                                                <  1:329356/90‑1 (MQ=255)
  gCACCCGTCTGGCCTCAGCTAATGGGCCGCGGGTAACATTTTTCAGACCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCaa                                                                                <  1:938859/90‑1 (MQ=255)
     cccGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGt                                                                             <  2:574512/90‑1 (MQ=255)
     cccGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGt                                                                             <  1:678854/90‑1 (MQ=255)
      ccGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGtt                                                                            >  2:35456/1‑90 (MQ=255)
       cGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTc                                                                           >  1:1223346/1‑90 (MQ=255)
       cGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTc                                                                           >  1:1560754/1‑90 (MQ=255)
        gTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCa                                                                          <  2:1318486/90‑1 (MQ=255)
        gTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCa                                                                          >  1:1445434/1‑90 (MQ=255)
        gTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCa                                                                          <  2:1087319/90‑1 (MQ=255)
         tCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAg                                                                         >  2:1494133/1‑90 (MQ=255)
            ggCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCaa                                                                                >  1:412949/1‑80 (MQ=255)
            ggCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCaa                                                                                <  2:412949/80‑1 (MQ=255)
             gCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAgttgt                                                                     <  2:1336014/90‑1 (MQ=255)
                tCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGc                                                                  <  1:636331/90‑1 (MQ=255)
                 cAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCa                                                                 >  1:1526578/1‑90 (MQ=255)
                 cAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCa                                                                 >  1:1524341/1‑90 (MQ=255)
                  aGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCAt                                                                <  1:905035/90‑1 (MQ=255)
                   gCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCAtt                                                               >  1:168153/1‑90 (MQ=255)
                           gTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAggg                                                       >  2:306465/1‑90 (MQ=255)
                           gTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGCTGTTGCATTCCAGAggg                                                       >  2:576195/1‑90 (MQ=255)
                             cgcgGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCa                                                     >  2:1507481/1‑90 (MQ=255)
                              gcgGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGttgttg                                                                   <  2:1123273/75‑1 (MQ=255)
                              gcgGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGttgttg                                                                   >  1:1123273/1‑75 (MQ=255)
                                    aaCATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGc                                              >  2:716683/1‑90 (MQ=255)
                                     aCATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCg                                             >  2:1031150/1‑90 (MQ=255)
                                             cAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCgc                                     <  1:239244/90‑1 (MQ=255)
                                               gATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCgctg                                   >  2:808075/1‑90 (MQ=255)
                                               gATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCgctg                                   >  1:1127357/1‑90 (MQ=255)
                                                  ccTCTGTGCTACCCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTc                                <  2:984322/90‑1 (MQ=255)
                                                  ccTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTc                                <  2:1381241/90‑1 (MQ=255)
                                                  ccTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTc                                >  2:1425846/1‑90 (MQ=255)
                                                    tctGTTCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTccc                              >  1:309669/1‑90 (MQ=255)
                                                    tctGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTccc                              >  2:573436/1‑90 (MQ=255)
                                                          cTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGa                        <  2:242259/90‑1 (MQ=255)
                                                            aaCCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGAcg                      >  2:1110673/1‑90 (MQ=255)
                                                            aaCCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGAcg                      >  2:324427/1‑90 (MQ=255)
                                                            aaCCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGAcg                      >  1:718067/1‑90 (MQ=255)
                                                            aaCCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGAcg                      >  2:68447/1‑90 (MQ=255)
                                                                 gccagTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCTCTGCTCCCGCCTGACGCGCCa                 >  2:913675/1‑90 (MQ=255)
                                                                 gccagTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCa                 >  1:1278812/1‑90 (MQ=255)
                                                                      tCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGaa            <  1:951002/90‑1 (MQ=255)
                                                                      tCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGaa            <  1:461875/90‑1 (MQ=255)
                                                                      tCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGaa            <  2:1526578/90‑1 (MQ=255)
                                                                               gcACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGc   >  2:1039328/1‑90 (MQ=255)
                                                                                cACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGCg  <  1:591946/90‑1 (MQ=255)
                                                                                        |                                                                                 
CAGCACCCGTCTGGCCTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGCG  >  NZ_CP009273/2888623‑2888792

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