Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,340,295 G→A 7.0% L194L (CTG→TTG) ‡ ribA ← GTP cyclohydrolase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,340,2950GA 7.0% 268.1 / 6.8 85L194L (CTG→TTG) ‡ribAGTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base G (47/32);  new base A (3/3);  total (50/35)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTACAA  >  AP009048/1340228‑1340360
                                                                   |                                                                 
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)
                   cTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGc                                              <  1:6028568/70‑1 (MQ=255)
                    ttCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCa                                                                 <  1:7442397/50‑1 (MQ=255)
                      ccGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:5031449/69‑1 (MQ=255)
                      ccGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCtt                                            <  1:68148/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:1208512/1‑62 (MQ=255)
                          caagcaaATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCgg                                               >  1:3297265/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:603649/1‑70 (MQ=255)
                                 aTGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGAt                                >  1:5917228/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:2976223/1‑71 (MQ=255)
                                              tGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTg                  >  1:3972281/1‑71 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:4861256/1‑58 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:7473043/1‑58 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:3096978/1‑58 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAAt                             >  1:303449/1‑58 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2088155/1‑60 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2487549/1‑60 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:6350047/1‑60 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:7122809/1‑60 (MQ=255)
                                                cAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGt                           >  1:2889967/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:4781571/55‑1 (MQ=38)
                                                        ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:5200994/55‑1 (MQ=38)
                                                        ccaTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGtt                      <  1:5191618/55‑1 (MQ=38)
                                                          aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:3200027/1‑57 (MQ=255)
                                                          aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:4634282/1‑57 (MQ=255)
                                                          aTTTGCTGAACAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTAt                    >  1:2060981/1‑57 (MQ=255)
                                                           tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCt      <  1:6351604/70‑1 (MQ=255)
                                                           tttGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCt      <  1:5333106/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:4194662/70‑1 (MQ=255)
                                                             tGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTAc    <  1:382965/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)
                                                                   |                                                                 
TTGCTTGCCGGTTATTTTGCTTCCGGCAAGCAAATGAATTACACAATGCAAGAGGGTTATTTGTTCAGCAAATGGCCCATTTTCTCGGCTTTGGTATCGAGATAATGTTCGTTATTGGGGTTACGACCTACAA  >  AP009048/1340228‑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: ‑