Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,564,985 T→C 11.1% F21L (TTT→CTT)  yihS → predicted glucosamine isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,564,9850TC11.1% 163.6 / 12.4 63F21L (TTT→CTT) yihSpredicted glucosamine isomerase
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (2/5);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

ATAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCA  >  AP009048/3564916‑3565051
                                                                     |                                                                  
aTAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCtt                                                                   >  1:4293772/1‑71 (MQ=255)
 tAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCtt                                                                   <  1:4429998/70‑1 (MQ=255)
     aGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAtt                                                              >  1:823562/1‑71 (MQ=255)
     aGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAtt                                                              <  1:1961512/71‑1 (MQ=255)
      ggAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt                                                             <  1:1289491/71‑1 (MQ=255)
         aTGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTgg                                                          >  1:2683154/1‑71 (MQ=255)
            aaaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTa                                                        <  1:712962/70‑1 (MQ=255)
            aaaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTa                                                        >  1:180208/1‑70 (MQ=255)
             aaTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa                                                      <  1:2079373/71‑1 (MQ=255)
              aTGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa                                                      >  1:3071793/1‑70 (MQ=255)
                  tttAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAAtt                                                  <  1:3839123/70‑1 (MQ=255)
                   ttAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTaa                                                       >  1:71471/1‑64 (MQ=255)
                    tAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTgg                                                          >  1:1789263/1‑60 (MQ=255)
                     aaCACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCg                                               <  1:3259557/70‑1 (MQ=255)
                      acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGCCCGCATCTTTg                                                                 >  1:1939196/1‑51 (MQ=255)
                      acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaaa                                                     >  1:3784561/1‑63 (MQ=255)
                      acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGt                                              <  1:1992138/70‑1 (MQ=255)
                      acacCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGt                                              <  1:2737687/70‑1 (MQ=255)
                         cccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt                                                             >  1:2540093/1‑52 (MQ=255)
                          ccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTc                                                 <  1:2476993/63‑1 (MQ=255)
                          ccTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGc                                         >  1:2970742/1‑71 (MQ=255)
                          ccTAAGCCACAACCTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAg                                            <  1:4409477/67‑1 (MQ=255)
                           cTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGcc                                        <  1:4180512/71‑1 (MQ=255)
                            tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt                                                             <  1:3723819/49‑1 (MQ=255)
                            tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt                                                             >  1:501572/1‑49 (MQ=255)
                            tAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGcc                                        >  1:2638092/1‑70 (MQ=255)
                              aGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAgaga                                         >  1:1660167/1‑64 (MQ=255)
                                  acaACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGt                                                         <  1:923471/47‑1 (MQ=255)
                                  acaACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTaaa                                                      <  1:1872682/50‑1 (MQ=255)
                                  acaACCGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGt                                 <  1:1425666/71‑1 (MQ=255)
                                   caACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtt                                <  1:1337454/71‑1 (MQ=255)
                                    aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAAt                                                   <  1:2737976/51‑1 (MQ=255)
                                    aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTcc                                                >  1:682267/1‑54 (MQ=255)
                                    aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGttt                               <  1:1679833/71‑1 (MQ=255)
                                    aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGttt                               <  1:4223309/71‑1 (MQ=255)
                                    aaCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGtt                                <  1:3578276/70‑1 (MQ=255)
                                     aCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGAttt                                                             >  1:1105966/1‑40 (MQ=255)
                                     aCCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCTTAGTGCCGACTGGtttt                              <  1:762792/71‑1 (MQ=255)
                                       cGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTgg                            <  1:606733/71‑1 (MQ=255)
                                       cGTTGGCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAg                                            <  1:1634492/55‑1 (MQ=255)
                                         ttGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTa                                             <  1:4499125/52‑1 (MQ=255)
                                           ggCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTgg                        >  1:3151819/1‑71 (MQ=255)
                                           ggCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTgg                        >  1:4535636/1‑71 (MQ=255)
                                           ggCTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTgg                        >  1:3454736/1‑71 (MQ=255)
                                             cTGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGt                                              <  1:1934764/47‑1 (MQ=255)
                                              tGGAACAGGAAACCGACCGCATCCTTGATTTTGGTAAAAATTCCGt                                              <  1:4593622/46‑1 (MQ=255)
                                                gAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAgg                   >  1:4248318/1‑71 (MQ=255)
                                                  aCAGGAAACCGTCCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAg                    <  1:4395735/68‑1 (MQ=255)
                                                  aCAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGc                           >  1:217913/1‑61 (MQ=255)
                                                   cAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTa                                             >  1:2423005/1‑42 (MQ=255)
                                                    aGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGc                                         >  1:653636/1‑45 (MQ=255)
                                                     ggAAACCGACAGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAg                    >  1:3008781/1‑65 (MQ=255)
                                                      gAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaa             <  1:1870827/71‑1 (MQ=255)
                                                       aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGa                                      >  1:2922118/1‑45 (MQ=255)
                                                       aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGa                                      >  1:917744/1‑45 (MQ=255)
                                                       aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaaa            <  1:4478676/71‑1 (MQ=255)
                                                       aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATaaa            <  1:894624/71‑1 (MQ=255)
                                                       aaaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATa              <  1:630018/69‑1 (MQ=255)
                                                        aaCCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAg           <  1:2633438/71‑1 (MQ=255)
                                                            gaccgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGc        >  1:2300993/1‑70 (MQ=255)
                                                             accgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCa       >  1:1724915/1‑70 (MQ=255)
                                                               cgCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAAt    >  1:4220883/1‑71 (MQ=255)
                                                                 cATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGttt                               >  1:449899/1‑42 (MQ=255)
                                                                 cATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCa  <  1:4207656/71‑1 (MQ=255)
                                                                   tCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCa  >  1:1303636/1‑69 (MQ=255)
                                                                     |                                                                  
ATAAAAGGAATGAAATGGTTTAACACCCTAAGCCACAACCGTTGGCTGGAACAGGAAACCGACCGCATCTTTGATTTTGGTAAAAATTCCGTAGTGCCGACTGGTTTTGGCTGGTTAGGCAATAAAGGGCAAATCA  >  AP009048/3564916‑3565051

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: