Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 231,008 T→G 6.4% intergenic (‑127/‑114) yafC ← / → yafD predicted DNA‑binding transcriptional regulator/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048231,0080TG 6.4% 163.5 / 6.6 63intergenic (‑127/‑114)yafC/yafDpredicted DNA‑binding transcriptional regulator/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/21);  new base G (2/2);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTC  >  AP009048/230938‑231077
                                                                      |                                                                     
tCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACt                                                                       >  1:2385274/1‑71 (MQ=255)
  aGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTg                                                                     >  1:3048742/1‑71 (MQ=255)
     aTGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCtt                                                                  <  1:2488257/71‑1 (MQ=255)
     aTGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCt                                                                   >  1:4674490/1‑70 (MQ=255)
       gCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTa                                                                 >  1:1523056/1‑70 (MQ=255)
        cGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTAc                                                                <  1:640310/70‑1 (MQ=255)
          tGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACtt                                                                      <  1:1426492/62‑1 (MQ=255)
          tGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATa                                                             <  1:4993042/71‑1 (MQ=255)
            atatGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATaa                                                            <  1:2826497/70‑1 (MQ=255)
            atatGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATaa                                                            >  1:3750984/1‑70 (MQ=255)
              atGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACaa                                                         >  1:59144/1‑71 (MQ=255)
              atGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACa                                                          >  1:5156349/1‑70 (MQ=255)
                   aTGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATa                                                             >  1:4679280/1‑62 (MQ=255)
                    tGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACa                                                   <  1:3596644/71‑1 (MQ=255)
                     gCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACaa                                                  >  1:991693/1‑71 (MQ=255)
                     gCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACaa                                                  >  1:4438344/1‑71 (MQ=255)
                      cGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACa                                                               <  1:418708/57‑1 (MQ=39)
                        aTTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTg                                               >  1:3016972/1‑71 (MQ=255)
                         tttAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGc                                              >  1:4958473/1‑71 (MQ=255)
                         tttAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGc                                              >  1:3030507/1‑71 (MQ=255)
                          ttAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACAt                                                              >  1:2242889/1‑54 (MQ=255)
                          ttAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCt                                             >  1:5062722/1‑71 (MQ=255)
                             aTGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACac                                                              <  1:3915449/51‑2 (MQ=38)
                                 tctcCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTg                                        >  1:1613826/1‑69 (MQ=255)
                                  ctcCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAtt                                                >  1:1596254/1‑60 (MQ=255)
                                  ctcCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAtt                                                >  1:3200645/1‑60 (MQ=255)
                                    ccATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACaa                                                  >  1:264200/1‑56 (MQ=255)
                                    ccATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGtt                                         >  1:1375991/1‑65 (MQ=255)
                                     cATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTa                                                                 <  1:1968415/40‑1 (MQ=255)
                                     cATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAagag                                  >  1:863238/1‑71 (MQ=255)
                                     cATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAag                                    <  1:3486183/69‑1 (MQ=255)
                                      aTAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACa                                                          >  1:2301146/1‑46 (MQ=255)
                                       tAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACAt                                                              <  1:4065742/41‑1 (MQ=255)
                                       tAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGa                                       <  1:4119541/64‑1 (MQ=255)
                                        aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACAt                                                              <  1:1963799/40‑1 (MQ=255)
                                        aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAtt                                                >  1:1229223/1‑54 (MQ=255)
                                        aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAtt                                                >  1:4241346/1‑54 (MQ=255)
                                        aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGc                                              >  1:4888372/1‑56 (MQ=255)
                                        aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGt                                 >  1:4978691/1‑69 (MQ=255)
                                         aTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAg                              >  1:3995355/1‑71 (MQ=255)
                                              cAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTa                           >  1:5074523/1‑69 (MQ=255)
                                                 aaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATc                      >  1:5048291/1‑71 (MQ=255)
                                                  aTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATc                      <  1:1659075/70‑1 (MQ=255)
                                                   ttCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGc                                              >  1:3642205/1‑45 (MQ=255)
                                                   ttCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATccc                    >  1:4830033/1‑71 (MQ=255)
                                                     cTGACCGGTTTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCta                  >  1:2986235/1‑71 (MQ=255)
                                                      tGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCta                  >  1:1591323/1‑70 (MQ=255)
                                                        aCCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATaa               <  1:484766/71‑1 (MQ=255)
                                                         ccGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCtata                <  1:2806203/69‑1 (MQ=255)
                                                         ccGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATaa               >  1:5076693/1‑70 (MQ=255)
                                                          cGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACa             <  1:1711017/71‑1 (MQ=255)
                                                            gtgtAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAg            >  1:3504920/1‑70 (MQ=255)
                                                            gtgtAAGCACTTGCTTACATAACAATATACAATAGCTCGTTGAAaga                                   >  1:4219753/1‑47 (MQ=255)
                                                            atgtAAGCAAGTGCTTACACAACAATATACAATTGCTCGTTGAAaga                                   >  1:1894433/2‑47 (MQ=25)
                                                            atgtAAGCAAGTGCTTACACAACAATATACAATTGCTCGTTGAAaga                                   >  1:1522256/2‑47 (MQ=25)
                                                              gtAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTaa                          <  1:1085154/54‑1 (MQ=255)
                                                              gtAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACa             >  1:1221296/1‑67 (MQ=255)
                                                               tAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAaga                                   <  1:1750761/44‑1 (MQ=255)
                                                               tAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTaaa                         <  1:1358451/54‑1 (MQ=255)
                                                               tAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAg         <  1:1101318/70‑1 (MQ=255)
                                                                aaGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCACTATAACAGTAGa        >  1:4706063/1‑70 (MQ=255)
                                                                 aGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTg                                <  1:3265520/45‑1 (MQ=255)
                                                                  gCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGaa       <  1:5128778/69‑1 (MQ=255)
                                                                  gCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAAcc     >  1:493722/1‑71 (MQ=255)
                                                                     cTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTc  <  1:2377671/71‑1 (MQ=255)
                                                                      |                                                                     
TCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTC  >  AP009048/230938‑231077

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