Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 302,171 A→G 19.4% intergenic (‑374/‑44) yagT ← / → yagU predicted xanthine dehydrogenase, 2Fe‑2S subunit/conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb302,1710AG19.4% 144.7 / 30.6 67intergenic (‑374/‑44)yagT/yagUpredicted xanthine dehydrogenase, 2Fe‑2S subunit/conserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/22);  new base G (8/5);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTGGCGATTAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCAC  >  W3110S.gb/302102‑302239
                                                                     |                                                                    
ttGGCGATTAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATAc                                                                     <  1:4178313/71‑1 (MQ=255)
ttGGCGATTAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATAc                                                                     <  1:2269677/71‑1 (MQ=255)
      aTTAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATGCTCCTg                                                                <  1:555160/70‑1 (MQ=255)
       ttAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGaa                                                              <  1:1152946/71‑1 (MQ=255)
       ttAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGa                                                               <  1:3330739/70‑1 (MQ=255)
        tAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATGCTCCTg                                                                >  1:1097823/1‑68 (MQ=255)
        tAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTg                                                                >  1:1279211/1‑68 (MQ=255)
           gAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAgg                                                          <  1:3619574/71‑1 (MQ=255)
            aaaCTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGgt                                                         >  1:1995241/1‑71 (MQ=255)
             aaCTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATGCTCCTGAAGAGgtg                                                        >  1:2047627/1‑71 (MQ=255)
             aaCTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTg                                                                >  1:2953657/1‑63 (MQ=255)
             aaCTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGgtg                                                        <  1:2840865/71‑1 (MQ=255)
              aCTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTg                                                                <  1:388745/62‑1 (MQ=255)
                  tCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGgtg                                                        >  1:2215768/1‑66 (MQ=255)
                  tCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGtata                                                    <  1:2956546/70‑1 (MQ=255)
                  tCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGtat                                                     <  1:2101959/69‑1 (MQ=255)
                   cTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGtatat                                                   >  1:3566395/1‑69 (MQ=255)
                    ttGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACa                                                 >  1:2859086/1‑71 (MQ=255)
                    ttGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACa                                                 >  1:713539/1‑71 (MQ=255)
                     tGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACa                                                 >  1:4085053/1‑70 (MQ=255)
                     tGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACa                                                 >  1:3203857/1‑70 (MQ=255)
                      gTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATaa                                                   <  1:1714573/67‑1 (MQ=255)
                      gTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACAt                                                <  1:2951640/70‑1 (MQ=255)
                         ataCAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTc                                                                   >  1:2515091/1‑48 (MQ=255)
                         ataCAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATg                                            >  1:3374592/1‑71 (MQ=255)
                         ataCAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATg                                            >  1:128512/1‑71 (MQ=255)
                          taCAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATgg                                           <  1:1560661/71‑1 (MQ=255)
                           aCAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAaga                                                            <  1:1764691/53‑1 (MQ=255)
                           aCAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAaga                                                            <  1:1350164/53‑1 (MQ=255)
                            cAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTa                                         >  1:1686699/1‑71 (MQ=255)
                             aaaaCAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTAc                                        <  1:3973567/71‑1 (MQ=255)
                              aaaCAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAaga                                                            >  1:2850606/1‑50 (MQ=255)
                                aCAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAaga                                                            >  1:3971154/1‑48 (MQ=255)
                                 cAATACAGTTCTTTACATTTGCCTTGTTTTATGAATGCTCCTg                                                                <  1:64235/43‑1 (MQ=255)
                                 cAATACAGTTCTTTACATTTGCCTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACATAATGGTACAAg                                     >  1:1983735/1‑70 (MQ=255)
                                   aTACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAg                                  >  1:2528910/1‑71 (MQ=255)
                                   aTACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAg                                  >  1:2512976/1‑71 (MQ=255)
                                       aGTTCTTTACATTTGCCTTGTTTTGTGAATACTCCTGAAGAGGTGTATAACa                                                 <  1:534334/52‑1 (MQ=255)
                                        gTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAggg                                >  1:4086441/1‑68 (MQ=255)
                                        gTTCTTTACATTTCCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAATCAggg                                >  1:4038419/1‑68 (MQ=255)
                                         ttCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAggg                                >  1:3455867/1‑67 (MQ=255)
                                         ttCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTCGa                            <  1:1992200/71‑1 (MQ=255)
                                           cTTTACATTTGCCTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACa                                                 <  1:1497852/48‑1 (MQ=255)
                                           cTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGata                          >  1:2331303/1‑71 (MQ=255)
                                           cTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGat                           <  1:2422721/70‑1 (MQ=255)
                                               aCATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATaa                                                   >  1:866477/1‑42 (MQ=255)
                                                 aTTTGCCTTGTTTTATGAATGCTCCTGAAGACGTGTATAACATa                                               >  1:2664875/1‑44 (MQ=255)
                                                    tGCCTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGaa                      >  1:3950567/1‑66 (MQ=255)
                                                    tGCCTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGa                       >  1:2777451/1‑65 (MQ=255)
                                                    tGCCTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGa                       >  1:2772103/1‑65 (MQ=255)
                                                    tGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATaa                                              >  1:2170685/1‑42 (MQ=255)
                                                     gCCTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATAtt                >  1:957026/1‑71 (MQ=255)
                                                       cTTGTTTTATGAATGCTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTg              <  1:1404839/71‑1 (MQ=255)
                                                       cTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTg              <  1:2637007/71‑1 (MQ=255)
                                                          gTTTTATGAATACTCCTGAAGAGGTGTATAACATaa                                              >  1:3004387/1‑36 (MQ=255)
                                                          gTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGaa            >  1:3539625/1‑70 (MQ=255)
                                                          gTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGaa            >  1:4502182/1‑70 (MQ=255)
                                                           ttttATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGAAATGAATATATTTGAACa          <  1:4494516/71‑1 (MQ=255)
                                                             ttATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGt                               >  1:2597106/1‑48 (MQ=255)
                                                             ttATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACaaa        >  1:32301/1‑71 (MQ=255)
                                                              tATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAAc       >  1:786505/1‑71 (MQ=255)
                                                               aTGAATACTCCAGAAGAGGTGTATAACATAATGGTACAAGCAg                                  >  1:1031496/1‑43 (MQ=255)
                                                                  aaTGCTCCTGAAGAGGTGTATAACATAATGGTACAAGCAggg                                <  1:4338702/42‑1 (MQ=255)
                                                                  aaTACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGaa                      <  1:2727637/52‑1 (MQ=255)
                                                                  aaTACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGaa                      <  1:163335/52‑1 (MQ=255)
                                                                    tACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAata                    >  1:944318/1‑52 (MQ=255)
                                                                    tACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCAc  >  1:2610147/1‑70 (MQ=255)
                                                                     |                                                                    
TTGGCGATTAGGAAACTATCTTGTTATACAAAACAATACAGTTCTTTACATTTGCCTTGTTTTATGAATACTCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCAC  >  W3110S.gb/302102‑302239

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