Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,053,471 A→G 28.4% F9S (TTT→TCT)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,053,4710AG28.4% 66.6 / 34.7 47F9S (TTT→TCT) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base A (11/22);  new base G (8/5);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

TGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGT  >  W3110S.gb/3053404‑3053533
                                                                   |                                                              
tGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAAct                                                             <  1:1711723/71‑1 (MQ=255)
  cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAActct                                                           <  1:3800460/71‑1 (MQ=255)
  cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAActct                                                           >  1:1592691/1‑71 (MQ=255)
      aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGcc                                                       >  1:4637555/1‑71 (MQ=255)
      aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGcc                                                       >  1:3456131/1‑71 (MQ=255)
                gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCtt                                                          <  1:1236565/58‑1 (MQ=255)
                gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCtt                                                          <  1:3203678/58‑1 (MQ=255)
                gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATAtct                                              >  1:4201751/1‑70 (MQ=255)
                   gcgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGaaa                                                               >  1:4250872/1‑50 (MQ=255)
                   gcgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATctcact                                          >  1:2106505/1‑71 (MQ=255)
                     gcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCggg                                                    >  1:103290/1‑59 (MQ=255)
                      cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCAt                                       <  1:4077075/71‑1 (MQ=255)
                      cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCAt                                       <  1:434373/71‑1 (MQ=255)
                       tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCAt                                       >  1:644661/1‑70 (MQ=255)
                        gCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATaa                                     >  1:3434758/1‑71 (MQ=255)
                              ttGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCgg                                                     <  1:733157/49‑1 (MQ=255)
                                gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGa                                                   >  1:3471002/1‑49 (MQ=255)
                                   tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCgg                                                     <  1:2597382/44‑1 (MQ=255)
                                   tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCAt                                       <  1:2919962/58‑1 (MQ=255)
                                   tttCCTCCACCATGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCAt                                       <  1:209920/58‑1 (MQ=255)
                                    ttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCgg                                                     <  1:308273/43‑1 (MQ=255)
                                    ttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTa                          >  1:3444310/1‑70 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACg                        >  1:3731593/1‑71 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAaca                                   >  1:2615433/1‑60 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAaca                                   >  1:3058595/1‑60 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAaca                                   >  1:859230/1‑60 (MQ=255)
                                       cTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtt                      <  1:2722581/71‑1 (MQ=255)
                                        tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCt                            <  1:649309/64‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATctcactca                                        <  1:2263086/47‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt                    <  1:4435024/67‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt                    <  1:140624/67‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGt                <  1:1386768/71‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGt                <  1:4187512/71‑1 (MQ=255)
                                              aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATctcactc                                         >  1:3878927/1‑45 (MQ=255)
                                              aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtt               >  1:3878919/1‑71 (MQ=255)
                                              aGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCtt                           >  1:4574794/1‑59 (MQ=255)
                                              aGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt                    <  1:730476/66‑1 (MQ=255)
                                                 gCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGttttt            <  1:1384933/71‑1 (MQ=255)
                                                  ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGttttt            <  1:4395774/70‑1 (MQ=255)
                                                   cTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt             >  1:2411232/1‑68 (MQ=255)
                                                    tGACGGCGACGCTGAGACTCTTGCCGGGATATctcactca                                        <  1:3648436/40‑1 (MQ=255)
                                                    tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt                    <  1:1199105/60‑1 (MQ=255)
                                                    tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt                    <  1:3710309/60‑1 (MQ=255)
                                                      aCGGCGACGCTGATACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGt                <  1:2381641/62‑1 (MQ=255)
                                                      aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttttt                  <  1:1000211/60‑1 (MQ=255)
                                                           gACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTCTTTGTTTTTAGTGTagagg  <  1:3138385/71‑2 (MQ=255)
                                                            aCGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttttt                  <  1:4412524/54‑1 (MQ=255)
                                                               cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt                    >  1:3982930/1‑49 (MQ=255)
                                                                   |                                                              
TGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGT  >  W3110S.gb/3053404‑3053533

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: