Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 58,662 G→A 24.6% V63V (GTG→GTA yabP → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904858,6620GA24.6% 107.9 / 33.5 57V63V (GTG→GTAyabPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (6/8);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GGAGACGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATGAGAGAGAA  >  AP009048/58593‑58731
                                                                     |                                                                     
ggAGACGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTg                                                                       <  1:2688425/70‑1 (MQ=255)
  agaCGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAAt                                                                    <  1:3576439/71‑1 (MQ=255)
    aCGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATaa                                                                  >  1:4286560/1‑71 (MQ=255)
       ttACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCa                                                               <  1:1181276/71‑1 (MQ=255)
       ttACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCa                                                               <  1:130695/71‑1 (MQ=255)
             cATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTgg                                                         <  1:110174/71‑1 (MQ=255)
              aTCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGc                                                        <  1:4228433/71‑1 (MQ=255)
              aTCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGc                                                        <  1:4233536/71‑1 (MQ=255)
              aTCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGc                                                        <  1:1312154/71‑1 (MQ=255)
                cTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCAGTTCAATTTCGTGAATAACCAGACTGGCaa                                                      >  1:4263959/1‑71 (MQ=255)
                 tGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGc                                                        >  1:2953097/1‑68 (MQ=255)
                     cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                                  <  1:4203722/70‑1 (MQ=255)
                     cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                                  <  1:3736416/70‑1 (MQ=255)
                     cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                                  <  1:3537273/70‑1 (MQ=255)
                        gTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATATCGTGAATAACCAGACTGGCAACATCAACg                                              >  1:2879971/1‑71 (MQ=255)
                         tCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                                  <  1:1593826/66‑1 (MQ=255)
                         tCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                                  <  1:3207613/66‑1 (MQ=255)
                         tCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAAcc                                                                <  1:914976/52‑1 (MQ=255)
                           aGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGcc                                            <  1:1260782/70‑1 (MQ=255)
                              tAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAg                                        <  1:3632422/71‑1 (MQ=255)
                              tAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAg                                        <  1:3417533/71‑1 (MQ=255)
                               aaaCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAg                                        >  1:1273233/1‑70 (MQ=255)
                               aaaCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAg                                        >  1:2888874/1‑70 (MQ=255)
                                 aCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACa                                                    >  1:2410985/1‑56 (MQ=255)
                                 aCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGAc                                                            <  1:2166498/48‑1 (MQ=255)
                                 aCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGAc                                                            <  1:3669643/48‑1 (MQ=255)
                                   aGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTgg                                                         >  1:714636/1‑49 (MQ=255)
                                   aGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATAc                                   >  1:1025582/1‑71 (MQ=255)
                                    gACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGAc                                                            >  1:2120597/1‑45 (MQ=255)
                                      cTGAGGCTTATAAACTTCTTTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACttt                                <  1:158327/71‑1 (MQ=255)
                                       tGAGGCTTATAAACTTCTGTTCAATTTCGTGAATaa                                                                  >  1:1723243/1‑36 (MQ=255)
                                       tGAGGCTTATAAACTTCTGTTCAATTCCGTGAATaa                                                                  <  1:757896/36‑1 (MQ=255)
                                         aGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGc                                                        <  1:2025423/44‑1 (MQ=255)
                                          ggCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGc                                             >  1:2001665/1‑54 (MQ=255)
                                           gCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTgg                           >  1:286485/1‑71 (MQ=255)
                                              tataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGc                        >  1:2350032/1‑71 (MQ=255)
                                               ataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCaa                                                >  1:422135/1‑46 (MQ=255)
                                               ataAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACTCCTCAGAATACt                                  >  1:235809/1‑60 (MQ=255)
                                                taAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCaa                                                      >  1:2402484/1‑39 (MQ=255)
                                                taAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCa                                         >  1:1337340/1‑52 (MQ=255)
                                                taAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTg                            >  1:3315072/1‑65 (MQ=255)
                                                 aaaCTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCaa                                                >  1:217949/1‑44 (MQ=255)
                                                 aaaCTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTg                            <  1:416193/64‑1 (MQ=255)
                                                     ttCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGc                                             >  1:3010307/1‑43 (MQ=255)
                                                     ttCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCaa                                                >  1:1439605/1‑40 (MQ=255)
                                                      tCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTa                               >  1:857182/1‑56 (MQ=255)
                                                       cTGTTGAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAAc                 >  1:1960610/1‑69 (MQ=255)
                                                       cTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAAc                 <  1:4417807/69‑1 (MQ=255)
                                                         gTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACtt                                 >  1:2458242/1‑51 (MQ=255)
                                                          ttCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACga               <  1:1608828/68‑1 (MQ=255)
                                                            cAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAAt           >  1:1877988/1‑70 (MQ=255)
                                                             aaTTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATga         <  1:3832651/71‑1 (MQ=255)
                                                             aaTTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATga         <  1:1638989/71‑1 (MQ=255)
                                                              aTTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATgag        <  1:1235608/71‑1 (MQ=255)
                                                               tttCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTg                            <  1:4600835/50‑1 (MQ=255)
                                                                ttCGTGAATAACCAGACTGGCAACATCAACGCCTc                                          >  1:1212484/1‑35 (MQ=255)
                                                                ttCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTAt                     >  1:2182764/1‑56 (MQ=255)
                                                                    tGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATGAGAGAGaa  <  1:4348569/71‑1 (MQ=255)
                                                                     gAATAACCAGACTGGCAACATCAACGCCTCAGAATACt                                  >  1:2457440/1‑38 (MQ=255)
                                                                     |                                                                     
GGAGACGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATGAGAGAGAA  >  AP009048/58593‑58731

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