Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,340,291 T→C 6.9% N195S (AAC→AGC)  ribA ← GTP cyclohydrolase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,340,2910TC 6.9% 245.8 / 8.1 87N195S (AAC→AGC) ribAGTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base T (49/32);  new base C (3/3);  total (52/35)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

ACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTACAA  >  AP009048/1340223‑1340360
                                                                    |                                                                     
acTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGtt                                                                      >  1:2484317/1‑70 (MQ=255)
  tATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAATTGAATTACACAATGCAAGAGGGTTATTTGTTCa                                                                    >  1:3292158/1‑70 (MQ=255)
     ttgcttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACCCAATGCAAGAGGGTTATTTGTTCAGCa                                                                 >  1:3167364/1‑70 (MQ=255)
       gcttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAt                                                              >  1:3210961/1‑71 (MQ=255)
        cttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:7219618/67‑1 (MQ=255)
        cttgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATg                                                             <  1:6767150/71‑1 (MQ=255)
          tgcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGc                                                           >  1:6071155/1‑71 (MQ=255)
           gcCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATg                                                             >  1:6429116/1‑68 (MQ=255)
            ccGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:795270/63‑1 (MQ=255)
            ccGGTCATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGccc                                                         <  1:4050754/71‑1 (MQ=255)
                ttATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCaaa                                                               >  1:863737/1‑61 (MQ=255)
                ttATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCaaa                                                               >  1:5327763/1‑61 (MQ=255)
                ttATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAttt                                                     >  1:5527855/1‑71 (MQ=255)
                 tATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAttt                                                     >  1:4312907/1‑70 (MQ=255)
                 tATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAttt                                                     >  1:5609298/1‑70 (MQ=255)
                   ttttGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGc                                                           >  1:4536052/1‑62 (MQ=255)
                   ttttGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTGATTTGTTCAGCAAATGGCCCATTTtc                                                   <  1:187716/70‑1 (MQ=255)
                    tttGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTtc                                                   <  1:31619/69‑1 (MQ=255)
                      tGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCa                                                                    >  1:7404150/1‑50 (MQ=255)
                        cTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGc                                              <  1:6028568/70‑1 (MQ=255)
                         ttCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:7442397/50‑1 (MQ=255)
                           ccGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:68148/69‑1 (MQ=255)
                           ccGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:5031449/69‑1 (MQ=255)
                            cGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTgg                                         >  1:5322413/1‑71 (MQ=255)
                             ggCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAAt                                                              >  1:5816580/1‑49 (MQ=255)
                              gcaagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       >  1:5138689/1‑71 (MQ=255)
                              gcaagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       >  1:4222810/1‑71 (MQ=255)
                               caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCgg                                               >  1:3297265/1‑62 (MQ=255)
                               caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCgg                                               >  1:1208512/1‑62 (MQ=255)
                               caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTAt                                      >  1:1889331/1‑71 (MQ=255)
                                aagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATc                                     >  1:6346793/1‑71 (MQ=255)
                                 agcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTtc                                                   <  1:7173522/56‑1 (MQ=255)
                                 agcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATc                                     >  1:1166856/1‑70 (MQ=255)
                                   caaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCt                                             <  1:1412664/60‑1 (MQ=255)
                                    aaaTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       <  1:2029671/65‑1 (MQ=255)
                                    aaaTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTAta                                     >  1:6725543/1‑66 (MQ=255)
                                    aaaTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCgaga                                 <  1:6332372/71‑1 (MQ=255)
                                      aTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                >  1:5917228/1‑70 (MQ=255)
                                      aTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                >  1:603649/1‑70 (MQ=255)
                                       tGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                <  1:2659261/69‑1 (MQ=255)
                                        gAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCAt                                                       <  1:4263042/45‑1 (MQ=255)
                                          aTTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:4278258/1‑71 (MQ=255)
                                           ttACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGtt                          >  1:4354366/1‑71 (MQ=255)
                                           ttACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           <  1:4396341/70‑1 (MQ=255)
                                                caATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTa                     <  1:3663820/71‑1 (MQ=255)
                                                caATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTa                     >  1:6171478/1‑71 (MQ=255)
                                                  aTGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAtt                   <  1:6781157/71‑1 (MQ=255)
                                                  aTGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAtt                   >  1:7755092/1‑71 (MQ=255)
                                                  aTGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAtt                   <  1:6124488/71‑1 (MQ=255)
                                                   tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:7213777/1‑71 (MQ=255)
                                                   tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:6683032/1‑71 (MQ=255)
                                                   tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:3972281/1‑71 (MQ=255)
                                                   tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:2976223/1‑71 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:7473043/1‑58 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:303449/1‑58 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:3096978/1‑58 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:4861256/1‑58 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2889967/1‑60 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:7122809/1‑60 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2487549/1‑60 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2088155/1‑60 (MQ=255)
                                                     cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:6350047/1‑60 (MQ=255)
                                                       agagGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:5126215/41‑1 (MQ=255)
                                                        gagGGTTATTTGTTAAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGCGt              >  1:1246257/1‑70 (MQ=255)
                                                           ggTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTa                                       <  1:4523446/42‑1 (MQ=255)
                                                           ggTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACg          <  1:1260752/71‑1 (MQ=255)
                                                            gTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    <  1:738768/60‑1 (MQ=255)
                                                            gTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGa         <  1:6937821/71‑1 (MQ=255)
                                                             ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:5191618/55‑1 (MQ=38)
                                                             ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:5200994/55‑1 (MQ=38)
                                                             ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:4781571/55‑1 (MQ=38)
                                                               aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:2060981/1‑57 (MQ=255)
                                                               aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:3200027/1‑57 (MQ=255)
                                                               aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:4634282/1‑57 (MQ=255)
                                                                tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCt      <  1:5333106/70‑1 (MQ=255)
                                                                tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCt      <  1:6351604/70‑1 (MQ=255)
                                                                tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATT‑GGGTTACGACCTa     >  1:5921977/1‑70 (MQ=255)
                                                                 ttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGt              <  1:1995570/61‑1 (MQ=255)
                                                                 ttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGt              <  1:1726977/61‑1 (MQ=255)
                                                                 ttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTa     >  1:32734/1‑70 (MQ=255)
                                                                  tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTAt                                      >  1:1447419/1‑36 (MQ=255)
                                                                  tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTTGTTATTGGGGTTACGACCTa     >  1:6742634/1‑69 (MQ=255)
                                                                  tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAc    <  1:382965/70‑1 (MQ=255)
                                                                  tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAc    <  1:4194662/70‑1 (MQ=255)
                                                                   gTTCAGCAAATGGCCCATTTTCTCGGGTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAcaa  >  1:1038347/1‑71 (MQ=255)
                                                                   gTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:7315210/1‑53 (MQ=255)
                                                                   gTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAcaa  >  1:4249525/1‑71 (MQ=255)
                                                                   gTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGGTACGACCTAca   <  1:5256991/70‑1 (MQ=255)
                                                                    |                                                                     
ACTATTTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTACAA  >  AP009048/1340223‑1340360

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: ‑