Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,351,422 A→C pseudogene (140/198 nt) BW25113_RS11725 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,351,4220AC98.0% 177.2 / ‑3.7 50pseudogene (140/198 nt)BW25113_RS11725hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (25/24);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-01

TATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGAT  >  NZ_CP009273/2351300‑2351560
                                                                                                                          |                                                                                                                                          
tATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGc                                                                                                                             >  2:59490/1‑138 (MQ=255)
  ttGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCa                                                                                                                           <  1:1178115/138‑1 (MQ=255)
  ttGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGcc                                                                                                                            <  1:921807/137‑1 (MQ=255)
        aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGc                                                                                                                     >  1:1100305/1‑138 (MQ=255)
             cTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCg                                                                                                               <  1:946268/139‑1 (MQ=255)
             cTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCg                                                                                                               <  1:1158019/139‑1 (MQ=255)
                          gaCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGc                                                                                                  <  2:242734/139‑1 (MQ=255)
                            cATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAg                                                                                                <  2:1100305/139‑1 (MQ=255)
                            cATCAAATGGGCGGGCCGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAg                                                                                                <  2:683205/139‑1 (MQ=255)
                                  aTTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAAt                                                                                           >  1:96270/1‑138 (MQ=255)
                                     ggctggcAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAAc                                                                                        >  2:137141/1‑138 (MQ=255)
                                     ggctggcAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAAc                                                                                        <  2:396532/138‑1 (MQ=255)
                                      gctggcAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCa                                                                                      <  1:833514/139‑1 (MQ=255)
                                                 ggaaaaTTAGAAGTTTTGCATTAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAacac                                                                           <  2:521846/136‑1 (MQ=255)
                                                 ggTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAacac                                                                           <  1:49690/139‑1 (MQ=255)
                                                             ggTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGcc                                                               >  2:167117/1‑139 (MQ=255)
                                                               tttGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGa                                                             <  1:117506/139‑1 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAg                                                                                                >  2:612463/1‑102 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAg                                                                                                <  1:612463/102‑1 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCa                                                                                      >  1:755414/1‑112 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCa                                                                                      <  2:755414/112‑1 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATc                                                           >  2:925030/1‑139 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATc                                                           >  2:841676/1‑139 (MQ=255)
                                                                 tGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATc                                                           >  1:1127847/1‑139 (MQ=255)
                                                                  gCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATcc                                                          >  1:1097253/1‑139 (MQ=255)
                                                                     tGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCa                                                       <  1:841676/139‑1 (MQ=255)
                                                                        aCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATgg                                                    >  1:1083129/1‑139 (MQ=255)
                                                                        aCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATgg                                                    >  2:1205370/1‑139 (MQ=255)
                                                                         cAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATggg                                                   <  1:468767/139‑1 (MQ=255)
                                                                             ccATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCtt                                               >  2:432395/1‑139 (MQ=255)
                                                                                 tAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGa                                            >  2:41643/1‑138 (MQ=255)
                                                                                   aaaTTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGatat                                         <  2:627710/139‑1 (MQ=255)
                                                                                             gcgcATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTa                                <  2:867057/138‑1 (MQ=255)
                                                                                             gcgcATGCTGGAACAGGGCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTa                                <  1:137141/138‑1 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCa                                                                                      <  2:1067946/68‑1 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCa                                                                                      >  1:1067946/1‑68 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAacac                                                                           >  2:50089/1‑79 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAacac                                                                           <  1:50089/79‑1 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGAt                                                            <  1:498734/94‑1 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGAt                                                            >  2:498734/1‑94 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGcc                >  2:408035/1‑138 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGcc                >  2:272352/1‑138 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGcc                >  2:586262/1‑138 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGcc                >  1:21144/1‑138 (MQ=255)
                                                                                                             ggCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGcc                >  1:874976/1‑138 (MQ=255)
                                                                                                              gCTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCa              >  1:1149307/1‑139 (MQ=255)
                                                                                                               cTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAAt                                                                                           >  1:1133757/1‑61 (MQ=255)
                                                                                                               cTTTGATCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAAt                                                                                           <  2:1133757/61‑1 (MQ=255)
                                                                                                                     tCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGc       <  2:379575/139‑1 (MQ=255)
                                                                                                                     tCGTGCCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGc       <  1:272352/139‑1 (MQ=255)
                                                                                                                          aCCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGAt  >  2:256664/1‑139 (MQ=255)
                                                                                                                          |                                                                                                                                          
TATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGAT  >  NZ_CP009273/2351300‑2351560

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