Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,759,232 (G)6→7 intergenic (‑97/‑45) yfjM ← / → rnlA DUF1524 domain‑containing protein/type II toxin‑antitoxin system toxin mRNA endoribonuclease RnlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,759,2261.G94.2% 184.9 / 4.4 52intergenic (‑91/‑51)yfjM/rnlADUF1524 domain‑containing protein/type II toxin‑antitoxin system toxin mRNA endoribonuclease RnlA
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (30/19);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

TTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTT  >  NZ_CP009273/2759137‑2759302
                                                                                          |                                                                            
ttGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCa                                                                               >  2:1884316/1‑90 (MQ=255)
 tgtgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                             <  1:96904/90‑1 (MQ=255)
   tgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                           <  1:2386280/90‑1 (MQ=255)
   tgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                           <  2:2238106/90‑1 (MQ=255)
   tgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                           <  2:2494577/90‑1 (MQ=255)
     gACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg                                                                         <  1:2257767/90‑1 (MQ=255)
     gACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg                                                                         <  2:688160/90‑1 (MQ=255)
       cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggggg                                                                       >  2:3018689/1‑89 (MQ=255)
       cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggggg                                                                       >  2:2637111/1‑89 (MQ=255)
       cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggggg                                                                       >  2:125367/1‑89 (MQ=255)
        tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggt                                                                      <  2:68901/90‑3 (MQ=255)
          cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggtat                                                                    <  1:1576739/90‑5 (MQ=255)
          cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggtat                                                                    <  2:2472326/90‑5 (MQ=255)
          cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggtat                                                                    <  1:1366920/90‑5 (MQ=255)
          cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggtat                                                                    <  1:1126566/90‑5 (MQ=255)
          cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggtat                                                                    <  2:1597845/90‑5 (MQ=255)
          cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggtat                                                                    <  1:1020787/90‑5 (MQ=255)
                   tAGACGATCCTGTCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTaga                                                            >  2:440855/1‑90 (MQ=255)
                   tAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTaga                                                            >  1:1944300/1‑90 (MQ=255)
                      aCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGtt                                                         <  2:987919/90‑1 (MQ=255)
                      aCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGtt                                                         <  2:1646095/90‑1 (MQ=255)
                      aCGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGtt                                                         <  1:1350820/90‑1 (MQ=255)
                       cGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGttt                                                        <  1:2637111/90‑1 (MQ=255)
                        gATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTc                                                       <  2:1124640/90‑1 (MQ=255)
                          tCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTccc                                                     <  2:2526653/90‑1 (MQ=255)
                          tCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTccc                                                     <  1:1884316/90‑1 (MQ=255)
                          tCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTccc                                                     <  2:944043/90‑1 (MQ=255)
                          tCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg                                                                         >  2:62313/1‑69 (MQ=255)
                          tCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggggg                                                                         <  1:62313/69‑1 (MQ=255)
                            cTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCa                                                   >  1:1062650/1‑90 (MQ=255)
                              gCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCata                                                 >  2:3017097/1‑90 (MQ=255)
                               cccTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCatat                                                <  2:906868/90‑1 (MQ=255)
                                  tataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGtt                                             <  2:3025801/90‑1 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  2:2274104/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  2:1552157/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  2:1748878/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  2:1075299/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  1:1862696/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  1:437581/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGttt                                            >  1:326000/1‑90 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggt                                                                      >  2:1876649/1‑61 (MQ=255)
                                   ataGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggggt                                                                      <  1:1876649/63‑3 (MQ=255)
                                       ttGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt                                        <  2:1471555/90‑1 (MQ=255)
                                        tGGAGTGCCCCCTTTTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTag                                                             <  2:888786/68‑1 (MQ=255)
                                        tGGAGTGCCCCCTTTTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTag                                                             >  1:888786/1‑68 (MQ=255)
                                        tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                       >  2:1497446/1‑90 (MQ=255)
                                        tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                       >  1:2882616/1‑90 (MQ=255)
                                        tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGCATTATTGAAGAGTTTCCCCATATGTTTCTATg                                       >  1:2532564/1‑90 (MQ=255)
                                        tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                           >  2:657246/1‑53 (MQ=255)
                                        tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                           <  1:657246/53‑1 (MQ=255)
                                        tGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATgg                                      <  2:708646/90‑1 (MQ=255)
                                           aGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGa                                    >  1:1001714/1‑90 (MQ=255)
                                             tGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtattat                                                                  >  2:2653784/1‑57 (MQ=255)
                                             tGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGCGGGGtattat                                                                  <  1:2653784/57‑1 (MQ=255)
                                              gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGttt                                                        <  1:101113/67‑1 (MQ=255)
                                              gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGttt                                                        >  2:101113/1‑67 (MQ=255)
                                              gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATcc                                 >  2:2020296/1‑90 (MQ=255)
                                              gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATcc                                 >  1:1439662/1‑90 (MQ=255)
                                              gCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATcc                                 >  2:60155/1‑90 (MQ=255)
                                                ccccTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAg                               >  2:2827888/1‑90 (MQ=255)
                                                   cTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaa                            >  2:2676862/1‑90 (MQ=255)
                                                     aTTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt                          <  2:1889627/90‑1 (MQ=255)
                                                       tCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGa                        <  1:3017097/90‑1 (MQ=255)
                                                         aGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACa                      >  1:1361541/1‑90 (MQ=255)
                                                             ggCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                  >  2:1022973/1‑90 (MQ=255)
                                                             ggCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                  >  2:675359/1‑90 (MQ=255)
                                                             ggCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                  >  2:778606/1‑90 (MQ=255)
                                                             ggCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                  >  2:1506079/1‑90 (MQ=255)
                                                             ggCTAAATTTAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                  >  1:1347341/1‑90 (MQ=255)
                                                                    tttAGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTa           >  2:2213467/1‑90 (MQ=255)
                                                                       aGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGtt                                                         <  2:459306/41‑1 (MQ=255)
                                                                       aGCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGtt                                                         >  1:459306/1‑41 (MQ=255)
                                                                        gCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa       <  1:2213467/90‑1 (MQ=255)
                                                                        gCAAATCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaaa       <  1:3018689/90‑1 (MQ=255)
                                                                             tCAATACACTTCAGGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACtt  <  1:2311476/90‑1 (MQ=255)
                                                                                          |                                                                            
TTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTT  >  NZ_CP009273/2759137‑2759302

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