Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,858,499 G→C W131C (TGG→TGC yeaD → D‑hexose‑6‑phosphate mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,858,4990GC98.2% 184.7 / ‑5.0 54W131C (TGG→TGCyeaDD‑hexose‑6‑phosphate mutarotase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (22/31);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGACGCTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGGCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCGATCTTGAATCACATGGCGAATTTGAAACC  >  NZ_CP009273/1858414‑1858583
                                                                                     |                                                                                     
ggACGCTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGTCGAAAAAATTCTGCCCGc                                                                                   <  2:403626/90‑1 (MQ=255)
ggACGCTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGc                                                                                   <  2:99201/90‑1 (MQ=255)
  aCGCTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCAc                                                                                 <  2:580141/90‑1 (MQ=255)
     cTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGAc                                                                              >  1:1385978/1‑90 (MQ=255)
     cTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGAc                                                                              >  1:1875412/1‑90 (MQ=255)
          aTCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTAc                                                                         <  2:1766488/90‑1 (MQ=255)
          aTCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTAc                                                                         <  2:2293263/90‑1 (MQ=255)
           tCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACg                                                                        <  2:2453559/90‑1 (MQ=255)
            cacaTCATGAAGATGCTGATGGCGTAGCACTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGc                                                                       <  1:954794/90‑1 (MQ=255)
                tcatGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGtt                                                                   <  2:881970/90‑1 (MQ=255)
                tcatGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGtt                                                                   <  1:2041148/90‑1 (MQ=255)
                tcatGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGtt                                                                   <  1:28042/90‑1 (MQ=255)
                 catGAAGATGCTTATAGCGTTGCGCTGACTTTTGAATTGACGCAAAGAGAAGAGCCGAAAAAATTCTGCCCGCACGACTTTACGCTGTTa                                                                  >  1:2293263/1‑90 (MQ=255)
                 catGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCAc                                                                                 <  1:51840/75‑1 (MQ=255)
                 catGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCAc                                                                                 >  2:51840/1‑75 (MQ=255)
                 catGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCAc                                                                                 >  2:907703/1‑75 (MQ=255)
                 catGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCAc                                                                                 <  1:907703/75‑1 (MQ=255)
                   tGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAgc                                                                >  1:345299/1‑90 (MQ=255)
                   tGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAgc                                                                >  1:73688/1‑90 (MQ=255)
                        aTGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCAtt                                                           >  1:1190101/1‑90 (MQ=255)
                        aTGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCAtt                                                           >  2:181242/1‑90 (MQ=255)
                          gCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTc                                                        >  2:140671/1‑90 (MQ=255)
                          gCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATATc                                                        >  1:1020439/1‑90 (MQ=255)
                             gATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGGACGACTTTACGCTGTTAGCGCATTTCcgc                                                     >  1:676754/1‑90 (MQ=255)
                              aTGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCcgcg                                                    <  2:2190882/90‑1 (MQ=255)
                              aTGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCcgcg                                                    >  2:721082/1‑90 (MQ=255)
                              aTGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCcgcg                                                    <  2:1190101/90‑1 (MQ=255)
                              aTGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCcgcg                                                    <  1:781210/90‑1 (MQ=255)
                               tGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGt                                                   <  1:870038/90‑1 (MQ=255)
                                 gCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACg                                                                        >  2:307589/1‑68 (MQ=255)
                                 gCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACg                                                                        <  1:307589/68‑1 (MQ=255)
                                    tAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCGGGGCACACGACTATACGATGTTAGCGCAGTTGCGCGCggaaa                                              >  1:2301047/1‑87 (MQ=255)
                                           gACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCCGCGTGGGTAAAACtt                                        >  2:943721/1‑90 (MQ=255)
                                           gACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTtg                                       >  2:917879/1‑90 (MQ=255)
                                              ttttGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGGACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGa                                    <  2:676754/90‑1 (MQ=255)
                                               tttGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGaa                                   >  1:2162462/1‑90 (MQ=255)
                                                  gAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATc                                <  2:2162462/90‑1 (MQ=255)
                                                  gAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATc                                <  2:388344/90‑1 (MQ=255)
                                                  gAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATc                                <  1:2448663/90‑1 (MQ=255)
                                                  gAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATc                                <  1:1439653/90‑1 (MQ=255)
                                                  gAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATc                                <  2:1533320/90‑1 (MQ=255)
                                                  gAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATc                                <  1:721082/90‑1 (MQ=255)
                                                   aaTTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCg                               <  1:189789/90‑1 (MQ=255)
                                                    aTTGACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCGa                              <  2:2279903/90‑1 (MQ=255)
                                                       gactcAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCGATCt                           >  1:962378/5‑90 (MQ=255)
                                                       gACGCAAAGCGAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCGATCt                           >  1:903669/1‑90 (MQ=255)
                                                                 gAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGtt                                                                   >  2:881986/1‑41 (MQ=255)
                                                                 gAAGAGACGAAAAAATTCTGCCCGCACGACTTTACGCTGtt                                                                   <  1:881986/41‑1 (MQ=255)
                                                                         gAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTaaaa                                           >  1:1616556/1‑56 (MQ=255)
                                                                         gAAAAAATTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTaaaa                                           <  2:1616556/56‑1 (MQ=255)
                                                                              aaTTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGt                                                   >  1:1848485/1‑43 (MQ=255)
                                                                              aaTTCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGt                                                   <  2:1848485/43‑1 (MQ=255)
                                                                                ttCTGCCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCGATCTTGAATCACATGGCGAATTTGAAacc  <  2:1523053/90‑1 (MQ=255)
                                                                                ttCTGCCCGCACGACTTTACGCTGTTAGCGCATATCCGCGTGGGTAAAACTTGTGAAATCGATCTTGAATCACATGGCGAATTTGAAacc  <  2:1020439/90‑1 (MQ=255)
                                                                                     |                                                                                     
GGACGCTGAAATCACATCATGAAGATGCTGATGGCGTAGCGCTGACTTTTGAATTGACGCAAAGCGAAGAGACGAAAAAATTCTGGCCGCACGACTTTACGCTGTTAGCGCATTTCCGCGTGGGTAAAACTTGTGAAATCGATCTTGAATCACATGGCGAATTTGAAACC  >  NZ_CP009273/1858414‑1858583

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