Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 5,593 A→T 12.7% intergenic (+63/+90) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported5,5930AT12.7% 182.5 / 17.7 71intergenic (+63/+90)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/25);  new base T (5/4);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

ACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCA  >  Exported/5523‑5663
                                                                      |                                                                      
aCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                        <  1:82975/71‑1 (MQ=255)
      aaaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                                  >  1:4630571/1‑71 (MQ=255)
          gACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGt                                                               >  1:3674851/1‑70 (MQ=255)
           aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGtt                                                              <  1:3212405/70‑1 (MQ=255)
           aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTa                                                             <  1:3452523/71‑1 (MQ=255)
            cAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                                            <  1:4412898/71‑1 (MQ=255)
            cAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                                                            >  1:4630314/1‑71 (MQ=255)
                 gCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                                                       >  1:1866947/1‑71 (MQ=255)
                     ggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                                                       <  1:910727/67‑1 (MQ=255)
                     ggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                                                       <  1:53032/67‑1 (MQ=255)
                         cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg                                               >  1:5540228‑M1/4‑53 (MQ=255)
                         aCAATCAGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatatat                                               >  1:740861/1‑71 (MQ=255)
                          cAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATAtt                                              <  1:2741661/71‑1 (MQ=255)
                            aTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTg                                             >  1:5468280/1‑70 (MQ=255)
                            aTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGa                                            <  1:4792640/71‑1 (MQ=255)
                               cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAgg                                                                      >  1:340860/1‑42 (MQ=255)
                               cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATc                                         >  1:4198018/1‑71 (MQ=255)
                               cGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggat                                         <  1:3323539‑M1/71‑25 (MQ=255)
                                ggCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggat                                         >  1:3154301‑M1/1‑46 (MQ=255)
                                 gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                       >  1:1033395/1‑71 (MQ=255)
                                  cATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                       >  1:3429588/1‑70 (MQ=255)
                                  cATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataag                                      <  1:577759‑M1/71‑28 (MQ=255)
                                    ttCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt                                        <  1:4489936/67‑1 (MQ=255)
                                    ttCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagc                                     >  1:264260‑M1/1‑42 (MQ=255)
                                    ttCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccgg                                           >  1:5316084‑M1/1‑42 (MQ=255)
                                     tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                                   <  1:5617861/39‑1 (MQ=255)
                                     tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCAt                                    >  1:1566045/1‑70 (MQ=255)
                                       aacgcctGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaag                                 >  1:2528465‑M1/3‑39 (MQ=255)
                                       aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt                                 >  1:626255/1‑71 (MQ=255)
                                        acgccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttg                                                   <  1:1843152‑M1/51‑15 (MQ=255)
                                        acgccTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaag                                 <  1:1266928‑M1/69‑33 (MQ=255)
                                          gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgc                              <  1:270646‑M1/71‑36 (MQ=255)
                                          gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcg                               >  1:5269286‑M1/1‑36 (MQ=255)
                                          gcCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgc                              <  1:4307580‑M1/71‑36 (MQ=255)
                                           ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgca                             <  1:1028253‑M1/71‑37 (MQ=255)
                                            cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgca                             <  1:3871861‑M1/70‑37 (MQ=255)
                                             tGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatc                           >  1:4757524‑M1/1‑33 (MQ=255)
                                             tGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccgg                                           <  1:2817629‑M1/55‑23 (MQ=255)
                                              gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg                          >  1:2140600/1‑71 (MQ=255)
                                              gATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgca                             >  1:3802728‑M1/1‑32 (MQ=255)
                                               aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                         >  1:2157516/1‑71 (MQ=255)
                                                tGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtag                                                >  1:1451296‑M1/1‑30 (MQ=255)
                                                tGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtag                                                >  1:1294376‑M1/1‑30 (MQ=255)
                                                tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                         >  1:3454739/1‑70 (MQ=255)
                                                 gCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatccgg                        >  1:1124834‑M1/1‑29 (MQ=255)
                                                  cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                      >  1:4075143/1‑71 (MQ=255)
                                                  cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                      >  1:4258947/1‑71 (MQ=255)
                                                  cGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                      <  1:5461251/71‑1 (MQ=255)
                                                  cGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggca                      >  1:5385333‑M1/1‑28 (MQ=255)
                                                   gACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatccggcat                     <  1:3170691‑M1/71‑45 (MQ=255)
                                                   gACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagcaaagcgcatccggcat                     <  1:3128084‑M1/71‑45 (MQ=255)
                                                   gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTg                                             <  1:697496/47‑1 (MQ=255)
                                                   gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAg                      >  1:3075199/1‑70 (MQ=255)
                                                     cGCTGGCGCGTCTTATCAGTCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                   <  1:2522751/71‑1 (MQ=255)
                                                      gCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacg                    >  1:781071‑M1/1‑24 (MQ=255)
                                                       cTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacg                    <  1:4337157‑M1/68‑46 (MQ=255)
                                                       cTGGCGCGTCTTATCAGGCCTACacgctgcgattttgta                                                 >  1:297465‑M1/1‑23 (MQ=255)
                                                       cTGGCGCGTCTTATCAGGCCTACacgctgcgattttgta                                                 >  1:3077778‑M1/1‑23 (MQ=255)
                                                        tGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                       <  1:4933170/48‑1 (MQ=255)
                                                          gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCa                                     >  1:4682794/1‑48 (MQ=255)
                                                          gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg               <  1:2819361/70‑1 (MQ=255)
                                                             cgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt           <  1:2253435/71‑1 (MQ=255)
                                                               tCTTATCTGGCCTACacgctgcgattttgtaggccggataagtaaagcgcatccgg                        <  1:239818‑M1/56‑42 (MQ=255)
                                                               tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa                  <  1:1069708/62‑1 (MQ=255)
                                                                 ttATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacgaa                  >  1:3674019‑M1/1‑13 (MQ=255)
                                                                 ttATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacgaa                  >  1:3667330‑M1/1‑13 (MQ=255)
                                                                 ttATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacgaa                  >  1:275120‑M1/1‑13 (MQ=255)
                                                                 ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTg                             >  1:467451/1‑49 (MQ=255)
                                                                  tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGt                            >  1:5188918/1‑49 (MQ=255)
                                                                  tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc      >  1:3893917/1‑71 (MQ=255)
                                                                   aTCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccgg                        >  1:182445‑M1/1‑11 (MQ=255)
                                                                   aTCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccgg                        >  1:3957820‑M1/1‑11 (MQ=255)
                                                                      aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc   <  1:2775898/70‑1 (MQ=255)
                                                                      aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc   <  1:2229800/70‑1 (MQ=255)
                                                                      aGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCa  >  1:4362773/1‑71 (MQ=255)
                                                                      |                                                                      
ACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCA  >  Exported/5523‑5663

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: