Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,224,919 G→A 7.5% intergenic (‑423/‑123) ymgG ← / → ymgH hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,224,9190GA 7.5% 167.9 / 3.4 53intergenic (‑423/‑123)ymgG/ymgHhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base A (2/2);  total (25/28)
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.73e-01

GATATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAAAA  >  AP009048/1224851‑1224986
                                                                    |                                                                   
gaTATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                   >  1:2303994/1‑71 (MQ=255)
   atTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGa                                                                >  1:4581760/1‑71 (MQ=255)
         caGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGta                                                           >  1:345771/1‑70 (MQ=255)
           gCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtatta                                                        <  1:2428986/71‑1 (MQ=255)
                gATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAg                                                    >  1:3291733/1‑70 (MQ=255)
                 aTTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                  <  1:2685975/71‑1 (MQ=255)
                  ttACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACa                                                              >  1:3117156/1‑58 (MQ=255)
                    aCATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGta                                                           >  1:1834425/1‑59 (MQ=255)
                     cATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTAGGACAGTATTATCAGAAAag                                               <  1:337039/70‑1 (MQ=255)
                          gTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGa                                          <  1:1250380/70‑1 (MQ=255)
                          gTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGa                                          <  1:863122/70‑1 (MQ=255)
                            cATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAttt                                       >  1:3566041/1‑71 (MQ=255)
                              ttttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGtatt                                                         <  1:3177254/51‑1 (MQ=255)
                                ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTAc                                                                  >  1:388365/1‑40 (MQ=255)
                                ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                  <  1:2573700/56‑1 (MQ=255)
                                ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAtttttt                                    <  1:2176622/70‑1 (MQ=255)
                                ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTa                                   >  1:1144001/1‑71 (MQ=255)
                                ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTa                                   >  1:3113394/1‑71 (MQ=255)
                                  attttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATc                                 >  1:3159854/1‑71 (MQ=255)
                                   ttttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAttttt                                     <  1:644116/66‑1 (MQ=255)
                                      tGCCGACGGCCTCATTGTCGAAAGATAAGAGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAAc                             >  1:3185368/1‑71 (MQ=255)
                                        ccGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACt                            <  1:3407922/70‑1 (MQ=255)
                                         cGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCaa                              >  1:2300986/1‑67 (MQ=255)
                                           aCGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaaaa                                                <  1:1355052/47‑1 (MQ=255)
                                           aCGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTacac                        >  1:367901/1‑71 (MQ=255)
                                           aCGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTacac                        >  1:1376402/1‑71 (MQ=255)
                                            cGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACt                       >  1:4009834/1‑71 (MQ=255)
                                              gCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAtttt                                      <  1:2125909/54‑1 (MQ=255)
                                               ccTCATTGTCGAAAGATAAGCGTACGACAGtattat                                                       <  1:2525986/36‑1 (MQ=255)
                                               ccTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAttttt                                     >  1:2600325/1‑54 (MQ=255)
                                                cTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                  <  1:1162030/40‑1 (MQ=255)
                                                 tCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAtttttt                                    >  1:817655/1‑53 (MQ=255)
                                                  cATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                  <  1:2792172/38‑1 (MQ=255)
                                                  cATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                  <  1:220278/38‑1 (MQ=255)
                                                  cATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAgcg                 >  1:2414267/1‑71 (MQ=255)
                                                   aTTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAgcgc                <  1:335325/71‑1 (MQ=255)
                                                   aTTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAgc                  <  1:1225221/69‑1 (MQ=255)
                                                     tGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCAc              >  1:4176830/1‑71 (MQ=255)
                                                      gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCaa                              >  1:1929358/1‑54 (MQ=255)
                                                      gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCAc              <  1:2598176/70‑1 (MQ=255)
                                                      gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACt             >  1:1181154/1‑71 (MQ=255)
                                                      gTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACt             >  1:370740/1‑71 (MQ=255)
                                                       tCGAAAGATAAGCGTACGACAGTATTATCAGaaaa                                                <  1:20477/35‑1 (MQ=255)
                                                       tCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACtt                      <  1:185676/61‑1 (MQ=255)
                                                       tCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACtt                      <  1:1644721/61‑1 (MQ=255)
                                                        cGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAtttttt                                    >  1:2508966/1‑46 (MQ=255)
                                                         gAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCg          <  1:2496130/71‑1 (MQ=255)
                                                          aaaGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACtt                      <  1:2677999/58‑1 (MQ=255)
                                                             gATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTa      <  1:4397477/71‑1 (MQ=255)
                                                                aaGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaa    >  1:3588280/1‑70 (MQ=255)
                                                                aaGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaa     <  1:4072723/69‑1 (MQ=255)
                                                                 aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaaa  <  1:1574449/71‑1 (MQ=255)
                                                                 aGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTaaaaa  <  1:2668209/71‑1 (MQ=255)
                                                                    gTACGACAGTATTATCAGAAAAGAGTGATTATTTATCCAACTACACTTCAGCGCACt             <  1:3236337/57‑1 (MQ=255)
                                                                    |                                                                   
GATATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAAAA  >  AP009048/1224851‑1224986

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: