Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,566,793 G→A 8.5% D889D (GAC→GAT)  pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,566,7930GA 8.5% 141.4 / 4.9 47D889D (GAC→GAT) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCAT  >  CP009273/1566655‑1566941
                                                                                                                                          |                                                                                                                                                    
gggCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAAt                                                                                                                                           <  1:429529/150‑1 (MQ=255)
   cTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAg                                                                                                                                        <  2:714417/150‑1 (MQ=255)
   cTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAg                                                                                                                                        <  2:372918/150‑1 (MQ=255)
         ttattTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTAtt                                                                                                                                  >  2:965188/1‑150 (MQ=255)
         ttattTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTAtt                                                                                                                                  <  1:535996/150‑1 (MQ=255)
              ttGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACa                                                                                                                              <  2:142906/149‑1 (MQ=255)
                  cAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgt                                                                                                                         <  2:746088/150‑1 (MQ=255)
                  cAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgt                                                                                                                         <  1:190602/150‑1 (MQ=255)
                    ataataCTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAgtgt                                                                                                                       >  1:700516/1‑150 (MQ=255)
                            ccGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAAt                                                                                                               >  1:991548/1‑150 (MQ=255)
                              ggTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAAttg                                                                                                             <  2:249804/150‑1 (MQ=255)
                                 ataAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAAttgtt                                                                                                           >  2:251452/1‑149 (MQ=255)
                                          tACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAacgctacg                                                                                                 >  1:93318/1‑150 (MQ=255)
                                          tACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAacgctac                                                                                                  <  2:480899/149‑1 (MQ=255)
                                           aCCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAacgctacgc                                                                                                >  2:630554/1‑150 (MQ=255)
                                                       aGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCACCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAAGAGGGGTCTTTAATTGCTGAACGCTACGCTGTTGGatacc                                                                                     >  1:215765/1‑147 (MQ=255)
                                                            ttGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGc                                                                                <  1:920610/149‑1 (MQ=255)
                                                                ttAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTgcgc                                                                           >  1:636931/1‑150 (MQ=255)
                                                                  aaCGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCt                                                                          <  2:700516/149‑1 (MQ=255)
                                                                     ggCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAAc                                                                      >  2:282296/1‑150 (MQ=255)
                                                                            ttAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTg                                                                <  1:311706/149‑1 (MQ=255)
                                                                              aaCATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTg                                                             >  1:8663/1‑150 (MQ=255)
                                                                                   tctcTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTatt                                                        >  1:855611/1‑150 (MQ=255)
                                                                                    ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTatt                                                        >  2:802845/1‑149 (MQ=255)
                                                                                     tctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCATCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattca                                                      >  2:683772/1‑150 (MQ=255)
                                                                                         ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcat                                                     <  2:441015/147‑1 (MQ=255)
                                                                                         ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcat                                                     <  2:944126/147‑1 (MQ=255)
                                                                                            tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGtt                                               <  2:353932/150‑1 (MQ=255)
                                                                                               tttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                                            >  2:204930/1‑150 (MQ=255)
                                                                                               tttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                                            >  1:811778/1‑150 (MQ=255)
                                                                                               tttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                                            <  2:127352/150‑1 (MQ=255)
                                                                                                      caacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCAc                                      <  1:204930/149‑1 (MQ=255)
                                                                                                       aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACt                                     <  1:965188/149‑1 (MQ=255)
                                                                                                         caaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGAt                                  >  1:447214/1‑150 (MQ=255)
                                                                                                               cTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCtt                             <  1:906111/149‑1 (MQ=255)
                                                                                                               cTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCtt                             >  2:872958/1‑149 (MQ=255)
                                                                                                               cTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCtt                             >  2:316372/1‑149 (MQ=255)
                                                                                                                 ccTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAg                          >  2:462493/1‑150 (MQ=255)
                                                                                                                  cTGCTCAGTCCATGCTGCAGGATCATCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTg                                    >  1:399881/1‑139 (MQ=255)
                                                                                                                   tGCTCAGTCCATGCTGCAGGATCATCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa                                   <  2:399881/139‑1 (MQ=255)
                                                                                                                     cTCAGTCCATGCTGCAGGATCATCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAg                      <  1:683772/150‑1 (MQ=255)
                                                                                                                            cATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTa               <  2:465625/150‑1 (MQ=255)
                                                                                                                              tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAAc             >  1:308638/1‑150 (MQ=255)
                                                                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAAc             >  1:960254/1‑149 (MQ=255)
                                                                                                                                    aGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCAc       <  1:261760/150‑1 (MQ=255)
                                                                                                                                     ggATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACt      <  1:471824/150‑1 (MQ=255)
                                                                                                                                     ggATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACt      >  1:732246/1‑150 (MQ=255)
                                                                                                                                          gTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCAt  <  1:885059/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GGGCTGTTATTATTTTGGCAATAATACTCCGGTATAAGTATTTACCGGATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCAT  >  CP009273/1566655‑1566941

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