Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,199,864 G→A pseudogene (239/246 nt)
pseudogene (8/1137 nt)
BW25113_RS26105 →
BW25113_RS05965 →
helix‑turn‑helix domain‑containing protein
terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,199,8640GA94.0% 170.7 / 1.0 51pseudogene (239/246 nt)
pseudogene (8/1137 nt)
BW25113_RS26105
BW25113_RS05965
helix‑turn‑helix domain‑containing protein
terminase large subunit
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (21/26);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAT  >  NZ_CP009273/1199729‑1199999
                                                                                                                                       |                                                                                                                                        
aCGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGtt                                                                                                                                        >  1:179289/1‑138 (MQ=255)
aCGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                                        <  2:524064/138‑1 (MQ=255)
aCGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                                        <  1:1251628/138‑1 (MQ=255)
aCGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                                        <  1:923329/138‑1 (MQ=255)
  gAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATtcat                                                                                                                                     <  2:471055/139‑1 (MQ=255)
         tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCa                                                                                                                              <  2:674834/139‑1 (MQ=255)
         tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCa                                                                                                                              <  2:1049808/139‑1 (MQ=255)
         tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCa                                                                                                                              <  1:439041/139‑1 (MQ=255)
           tgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACt                                                                                                                            <  2:861422/139‑1 (MQ=255)
                gCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGcc                                                                                                                       >  1:244348/1‑139 (MQ=255)
                   aTCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCATCCCACTGTGCCgtg                                                                                                                    <  2:52803/139‑1 (MQ=255)
                     cAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCgtggt                                                                                                                 >  1:1039697/1‑139 (MQ=255)
                     cAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCgtggt                                                                                                                 >  2:150756/1‑139 (MQ=255)
                                  tGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCAc                                                                                                    >  2:248918/1‑139 (MQ=255)
                                       gTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCa                                                                                               <  1:150756/139‑1 (MQ=255)
                                        tCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCAc                                                                                              >  2:64734/1‑139 (MQ=255)
                                                gcgcGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAg                                                                                       <  2:1221288/138‑1 (MQ=255)
                                                 cgcgGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                                     <  1:1243208/139‑1 (MQ=255)
                                                         tGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTa                                                                             <  2:67224/139‑1 (MQ=255)
                                                               ccGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTAcacca                                                                        <  1:1265037/138‑1 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCGGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCgtg                                                                                                                    <  2:848661/92‑1 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                                        >  2:211902/1‑72 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                                        <  1:211902/72‑1 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATtcat                                                                                                                                     <  1:186972/75‑1 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATtcat                                                                                                                                     >  2:186972/1‑75 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCgtg                                                                                                                    >  1:848661/1‑92 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTAcacca                                                                        <  1:716795/135‑1 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTAcacca                                                                        >  2:716795/1‑135 (MQ=255)
                                                                  tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGAt                                                                    >  1:404383/1‑139 (MQ=255)
                                                                      aaCGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCt                                                                 >  1:331252/1‑138 (MQ=255)
                                                                                     gCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGgcg                                                 <  1:248918/139‑1 (MQ=255)
                                                                                          aaGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACg                                            <  2:331252/139‑1 (MQ=255)
                                                                                              aGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                                     <  2:1257924/94‑1 (MQ=255)
                                                                                              aGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                                     >  1:1257924/1‑94 (MQ=255)
                                                                                                          accgcgAAGGGCGGGCCGCCGGGGTGTTATTTCTTCCCCCCACGGGCCCTGGTTGGAGAAATATCCCGAGCACGCCACAAATGCGTTTTCCCCCCGGATGCTTACCGGGAGGGGGCGCGACGCCAGCCACTGATGTgg                            <  1:64734/136‑1 (MQ=255)
                                                                                                            cgcgTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGc                          >  1:996771/1‑139 (MQ=255)
                                                                                                                tAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACg                                            <  1:906319/117‑1 (MQ=255)
                                                                                                                tAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACg                                            >  2:906319/1‑117 (MQ=255)
                                                                                                                 aaGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTa                     >  1:822784/1‑139 (MQ=255)
                                                                                                                    gcgcGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTacca                  <  2:1039697/139‑1 (MQ=255)
                                                                                                                       cgTCGCCAGGGTGATGATTCATCACCCCACTGTGc                                                                                                                        <  1:1045637/35‑1 (MQ=255)
                                                                                                                       cgTCGCCAGGGTGATGATTCATCACCCCACTGTGc                                                                                                                        >  2:1045637/1‑35 (MQ=255)
                                                                                                                       cgTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccg               <  2:822784/139‑1 (MQ=255)
                                                                                                                        gTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTAc                                                              >  1:952055/1‑91 (MQ=255)
                                                                                                                        gTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTAc                                                              <  2:952055/91‑1 (MQ=255)
                                                                                                                        gTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgc              >  2:1021486/1‑139 (MQ=255)
                                                                                                                         tCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgcc             >  2:35603/1‑139 (MQ=255)
                                                                                                                         tCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgcc             <  2:996771/139‑1 (MQ=255)
                                                                                                                                    aTGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAt  >  1:1240810/1‑139 (MQ=255)
                                                                                                                                    aTGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAt  >  2:496292/1‑139 (MQ=255)
                                                                                                                                    aTGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAt  >  2:931763/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                        
ACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAT  >  NZ_CP009273/1199729‑1199999

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