Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,271,272 T→C 7.6% L190P (CTC→CCC)  tyrB → tyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,271,2720TC 7.6% 220.6 / 9.1 79L190P (CTC→CCC) tyrBtyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base C (4/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGA  >  AP009048/4271202‑4271338
                                                                      |                                                                  
cTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtct                                                                    <  1:1126062/71‑1 (MQ=255)
   aaaaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtctc                                                                   <  1:3727634/69‑1 (MQ=255)
    aaaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCAc                                                                 >  1:1672128/1‑70 (MQ=255)
     aaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtct                                                                    <  1:1131848/66‑1 (MQ=255)
     aaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCAc                                                                 >  1:701761/1‑69 (MQ=255)
     aaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACt                                                                >  1:1857480/1‑70 (MQ=255)
     aaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTa                                                               <  1:129698/71‑1 (MQ=255)
      aCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTaa                                                              >  1:428713/1‑71 (MQ=255)
       cATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCa                                                                  >  1:4272304/1‑66 (MQ=255)
       cATTACCTGACCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtctc                                                                   >  1:625481/1‑65 (MQ=255)
        aTTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATg                                                            <  1:2769950/71‑1 (MQ=255)
          tACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGAt                                                          >  1:3152362/1‑71 (MQ=255)
           aCCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGAt                                                          >  1:27428/1‑70 (MQ=255)
            ccTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGAt                                                          >  1:3551249/1‑69 (MQ=255)
            ccTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATc                                                         <  1:3495985/70‑1 (MQ=255)
             cTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATc                                                         <  1:1557821/69‑1 (MQ=255)
             cTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAg                                                       >  1:3017628/1‑71 (MQ=255)
             cTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAg                                                       <  1:3287310/71‑1 (MQ=255)
              tGCCCGGAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCa                                                                  >  1:670247/1‑59 (MQ=255)
              tGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCa                                                                  >  1:676466/1‑59 (MQ=255)
               gCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACACCCCAACGGGTGCCGATCTCACTAATGATCAGTg                                                     >  1:2838835/1‑71 (MQ=255)
               gCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCa                                                        >  1:4402769/1‑68 (MQ=255)
               gCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCCCACTGATCATTAGt                                                      >  1:677179/1‑70 (MQ=255)
                cccGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTgg                                                    >  1:1082672/1‑71 (MQ=255)
                cccGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTg                                                     <  1:1230250/70‑1 (MQ=255)
                   gCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGAt                                                 >  1:1235772/1‑71 (MQ=255)
                    cAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATg                                                >  1:1201902/1‑71 (MQ=255)
                     aGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtct                                                                    >  1:2199556/1‑50 (MQ=255)
                     aGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCa                                                                  >  1:1750862/1‑52 (MQ=255)
                     aGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGAt                                                 >  1:1939055/1‑69 (MQ=255)
                     aGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGc                                               <  1:3758130/71‑1 (MQ=255)
                         ttGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGGCGATCTCACTAATGATCAGTGGGATGCgg                                             >  1:1346668/1‑69 (MQ=255)
                         ttGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCAc                                                                 >  1:1694052/1‑49 (MQ=255)
                         ttGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTa                                                               <  1:4577632/51‑1 (MQ=255)
                         ttGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGt                                            <  1:95554/70‑1 (MQ=255)
                         ttGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCCCACTGATCATTAGTgag                                                   <  1:2950483/63‑3 (MQ=37)
                         ttGTGTTGCTGCATCCATGTTGACACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGt                                            <  1:4371826/70‑1 (MQ=255)
                          tgtgTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTgg                                                    <  1:1988393/61‑1 (MQ=255)
                           gtgtTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGa                                          <  1:3510664/70‑1 (MQ=255)
                           gtgtTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGAt                                         >  1:4201746/1‑71 (MQ=255)
                           gtgtTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCCCACTGATCATTAGTGAGAt                                                 >  1:1786684/1‑63 (MQ=25)
                              ttGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTg                                       >  1:4452186/1‑70 (MQ=255)
                              ttGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGa                                      <  1:4943197/71‑1 (MQ=255)
                                gctgcATCCATGTTGCCACAACCCAACGGGTGCCGATCCCACTGATCATTAGTGAGAtcg                                               >  1:2747003/1‑58 (MQ=21)
                                 ctgcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGAt                                                          >  1:5028604/1‑48 (MQ=255)
                                 ctgcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGaaa                                    <  1:1614479/70‑1 (MQ=255)
                                    cATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAAt                                                             >  1:2264110/1‑42 (MQ=255)
                                    cATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAAtt                                  <  1:5160243/69‑1 (MQ=255)
                                     aTCCATGTTGCCACAACCCAACGGGTGCCGATCTCa                                                                  <  1:824804/36‑1 (MQ=255)
                                     aTCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAAt                                                             >  1:4020319/1‑41 (MQ=255)
                                      tCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCa                              <  1:4045594/71‑1 (MQ=255)
                                      tCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCAGTGATTGa                                      <  1:4687788/63‑1 (MQ=255)
                                       ccATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCa                              >  1:2405055/1‑70 (MQ=255)
                                        cATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGaaa                                    >  1:419802/1‑63 (MQ=255)
                                        cATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCaa                             <  1:3155511/70‑1 (MQ=255)
                                           gTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATtctc                               >  1:2170192/1‑65 (MQ=255)
                                            ttGCCACAACCCAACGTGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGa                                      <  1:4719678/57‑1 (MQ=255)
                                            ttGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATtct                                <  1:2707068/63‑1 (MQ=255)
                                             tGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGaaa                                    >  1:305705/1‑58 (MQ=255)
                                              gCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCcgc                      >  1:3567596/1‑71 (MQ=255)
                                               ccACATCCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCcgc                      <  1:1255203/70‑1 (MQ=255)
                                               ccACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGaaa                                    <  1:2856170/56‑1 (MQ=255)
                                               ccACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGaaa                                    <  1:3052372/56‑1 (MQ=255)
                                                cacaACCCAACGGGTGCCGATCCCACTGATCATTAGTgag                                                   <  1:765851/40‑3 (MQ=9)
                                                 acaACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAAtt                                  >  1:749214/1‑56 (MQ=255)
                                                 acaACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAg                   >  1:4516172/1‑71 (MQ=255)
                                                  caACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGc                  >  1:4140236/1‑71 (MQ=255)
                                                  caACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGc                  >  1:3238032/1‑71 (MQ=255)
                                                   aaCCCAACGGGTGCCGATCTCACTAATGATCAGTGGGAt                                                 >  1:2365293/1‑39 (MQ=255)
                                                   aaCCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATg                                                >  1:3196639/1‑40 (MQ=255)
                                                    aCCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCt                 <  1:3658522/70‑1 (MQ=255)
                                                     cccAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAAtt                                  >  1:3786668/1‑52 (MQ=255)
                                                     cccAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTa               <  1:4130075/71‑1 (MQ=255)
                                                      ccAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATtc                                 <  1:1082934/52‑1 (MQ=255)
                                                        aaCGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTattc            <  1:3572500/71‑1 (MQ=255)
                                                        aaCGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTattc            >  1:3241269/1‑71 (MQ=255)
                                                        aaCGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTattc            <  1:68833/71‑1 (MQ=255)
                                                        aaCGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTattc            <  1:2652641/71‑1 (MQ=255)
                                                           gggTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGa                                      <  1:338234/42‑1 (MQ=255)
                                                           gggTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTattcca          <  1:331009/70‑1 (MQ=255)
                                                               gCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTattccattcc      <  1:1161861/70‑1 (MQ=255)
                                                                ccGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTc    <  1:2286689/71‑1 (MQ=255)
                                                                 gcaTCCCACTGATCATTAGTGAGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTc    >  1:4093903/3‑70 (MQ=25)
                                                                  gATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGa  >  1:1480518/1‑71 (MQ=255)
                                                                     ctcACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCaaa                            >  1:4667695/1‑42 (MQ=255)
                                                                     ctcACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGccc                        >  1:2894187/1‑46 (MQ=255)
                                                                      |                                                                  
CTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGA  >  AP009048/4271202‑4271338

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: