Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,941,604 A→C G10G (GGT→GGG ruvC ← crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,941,6040AC91.7% 139.1 / 2.4 48G10G (GGT→GGGruvCcrossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base C (27/17);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941517‑1941686
                                                                                       |                                                                                  
atcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAccc                                                                                  >  2:208327/1‑90 (MQ=255)
 tcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 >  2:86207/1‑90 (MQ=255)
             cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  2:727277/90‑1 (MQ=255)
              ggATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCCCGCGCGACCCCGGGTCAATGccta                                                                    >  1:597491/1‑88 (MQ=255)
              ggATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCga                                                                    >  1:265697/1‑90 (MQ=255)
               gATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCgag                                                                   >  2:1567149/1‑90 (MQ=255)
                        gCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGc                                                          >  2:1378780/1‑90 (MQ=255)
                              cAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCAc                                                    <  2:847021/90‑1 (MQ=255)
                                ggTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCAcgc                                                  >  1:518320/1‑90 (MQ=255)
                                   aGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtc                                               <  2:1333163/90‑1 (MQ=255)
                                   aGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtc                                               >  2:584328/1‑90 (MQ=255)
                                     gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtctc                                             >  2:872261/1‑90 (MQ=255)
                                     gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtctc                                             >  1:245620/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATc                                                                            >  2:676906/1‑55 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATc                                                                            <  1:676906/55‑1 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     >  2:60985/1‑62 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  1:60985/62‑1 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:334875/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  2:1428128/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  1:808021/1‑90 (MQ=255)
                                         gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  1:463579/1‑90 (MQ=255)
                                             cctacctaCCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGc                                     <  1:1040150/90‑1 (MQ=255)
                                              ctacctacCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCt                                    >  1:1327313/1‑90 (MQ=255)
                                                     cctGGCGGATGACGCCGTAGCCGGTCACGCCAGACCCCGGATAAATGCCGAGAATAATAGCCATCACGCGTTTCCGTTTTGCTGTTTAGc                             <  1:1198035/90‑1 (MQ=255)
                                                      ctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCa                            >  2:560599/1‑90 (MQ=255)
                                                      ctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCa                            >  2:24048/1‑90 (MQ=255)
                                                       tGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtttt                                       <  1:614772/78‑1 (MQ=255)
                                                       tGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtttt                                       >  2:614772/1‑78 (MQ=255)
                                                       tGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  1:122212/90‑1 (MQ=255)
                                                       tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  2:1114538/90‑1 (MQ=255)
                                                       tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  2:923710/90‑1 (MQ=255)
                                                         gCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtttt                                       >  2:489761/1‑76 (MQ=255)
                                                         gCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtttt                                       <  1:489761/76‑1 (MQ=255)
                                                          cGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGcc                        >  1:1076659/1‑90 (MQ=255)
                                                          cGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGcc                        <  1:1345513/90‑1 (MQ=255)
                                                           ggATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCt                       >  1:351546/1‑90 (MQ=255)
                                                                 cGCCGTGGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCag                 <  2:862874/90‑1 (MQ=255)
                                                                 cGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCag                 <  2:326086/90‑1 (MQ=255)
                                                                  gCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCaga                >  2:53212/1‑90 (MQ=255)
                                                                     gTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCagagag             >  1:915000/1‑90 (MQ=255)
                                                                      tAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGt            <  2:536344/90‑1 (MQ=255)
                                                                      tAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGt            >  1:379429/1‑90 (MQ=255)
                                                                           ggTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTg       >  2:183324/1‑90 (MQ=255)
                                                                           ggTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTg       >  2:130311/1‑90 (MQ=255)
                                                                             tCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCg     >  1:462003/1‑90 (MQ=255)
                                                                               aCGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAc   >  2:1411420/1‑90 (MQ=255)
                                                                                cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  2:575507/90‑1 (MQ=255)
                                                                                cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  2:523064/90‑1 (MQ=255)
                                                                                cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  1:569663/90‑1 (MQ=255)
                                                                                cgcgcgAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  2:989211/90‑1 (MQ=255)
                                                                                       |                                                                                  
ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941517‑1941686

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