Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,583,583 T→C 6.1% N219S (AAC→AGC)  ydeN ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,583,5830TC 6.1% 272.4 / 6.9 98N219S (AAC→AGC) ydeNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/51);  new base C (3/3);  total (44/54)
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.99e-01

Only 100 of 109 total aligned reads displayed.

AACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCTTC  >  AP009048/1583514‑1583653
                                                                     |                                                                      
aaCAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGtt                                                                       <  1:1411161/71‑1 (MQ=255)
     tGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGtttt                                                                     >  1:5677916/1‑68 (MQ=255)
       gaGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGt                                                                >  1:1616914/2‑71 (MQ=255)
       aaGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGtt                                                                       >  1:4001354/1‑64 (MQ=255)
       aaGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAgg                                                                 <  1:3427935/70‑1 (MQ=255)
       aaGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGt                                                                >  1:4022164/1‑71 (MQ=255)
               ttGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttc                                                        >  1:5515126/1‑71 (MQ=255)
                 gTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGtt                                                               >  1:4290927/1‑62 (MQ=255)
                 gTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCa                                                           <  1:3536923/66‑1 (MQ=255)
                 gTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttctt                                                      >  1:1573705/1‑71 (MQ=255)
                  tAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttcttc                                                     <  1:5689765/71‑1 (MQ=255)
                  tAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttctt                                                      <  1:3316718/70‑1 (MQ=255)
                  tAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttctt                                                      <  1:142639/70‑1 (MQ=255)
                   aaTATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTg                                                                    <  1:1562263/55‑1 (MQ=255)
                   aaTATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTg                                                                    <  1:5149254/55‑1 (MQ=255)
                   aaTATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTg                                                                    <  1:1916568/55‑1 (MQ=255)
                   aaTATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttcttc                                                     <  1:4808822/70‑1 (MQ=255)
                     tatGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGtt                                                                       <  1:4185527/50‑1 (MQ=255)
                     tatGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCAttc                                                        <  1:1987182/65‑1 (MQ=255)
                       tGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCg                                                   >  1:4461089/1‑68 (MQ=255)
                          cGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaa                                             <  1:128196/71‑1 (MQ=255)
                              ccTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGa                                                                   >  1:2293436/1‑45 (MQ=255)
                              ccTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAtgt                                         >  1:1949718/1‑71 (MQ=255)
                              ccTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAtg                                          >  1:1964725/1‑70 (MQ=255)
                              ccTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAtg                                          <  1:5622255/70‑1 (MQ=255)
                              ccTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAt                                           >  1:3542388/1‑69 (MQ=255)
                               cTGCAGCGTGGAATCTCATAAAGTAATCAAAGCCACGGt                                                                        <  1:5839476/39‑1 (MQ=255)
                                  cAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCg                                                   >  1:1576407/1‑57 (MQ=255)
                                  cAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCg                                                   >  1:2609237/1‑57 (MQ=255)
                                  cAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCg                                                   >  1:4700640/1‑57 (MQ=255)
                                      gTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaa                                              <  1:7209337/58‑1 (MQ=255)
                                       tGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAg                                                                  <  1:270569/37‑1 (MQ=255)
                                       tGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAg                                                                  <  1:7674799/37‑1 (MQ=255)
                                       tGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCa                                                           <  1:4984526/44‑1 (MQ=255)
                                       tGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCa                                                           <  1:3831895/44‑1 (MQ=255)
                                       tGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGtt                                <  1:7576591/71‑1 (MQ=255)
                                        ggAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            >  1:5664345/1‑58 (MQ=255)
                                          aaTCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaa                                              <  1:6001086/54‑1 (MQ=255)
                                          aaTCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaa                                              <  1:3077900/54‑1 (MQ=255)
                                          aaTCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaa                                              <  1:6430803/54‑1 (MQ=255)
                                          aaTCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGt                              >  1:4842763/1‑70 (MQ=255)
                                           aTCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCa                            >  1:7938394/1‑71 (MQ=255)
                                            tCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCAt                           >  1:4670889/1‑71 (MQ=255)
                                              ccATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGa                         >  1:2721647/1‑71 (MQ=255)
                                               cATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCg                                                   <  1:617384/44‑1 (MQ=255)
                                               cATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGAt                        <  1:474212/71‑1 (MQ=255)
                                                aTAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGa                         >  1:2359569/1‑69 (MQ=255)
                                                 tAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            <  1:6061026/49‑1 (MQ=255)
                                                 tAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAg                      >  1:2127334/1‑71 (MQ=255)
                                                  aaaGTAATCAAAGCCACGGTTTTGAGGTTGCCAttctt                                                      >  1:2126481/1‑38 (MQ=255)
                                                  aaaGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGt                     >  1:6853791/1‑71 (MQ=255)
                                                   aaGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            >  1:3896786/1‑47 (MQ=255)
                                                   aaGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            >  1:1649981/1‑47 (MQ=255)
                                                   aaGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTGCCGCAGAAAATGTGGTGAAGTTGTCATGATAGt                     >  1:7198253/1‑70 (MQ=255)
                                                    aGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCa                   >  1:3979236/1‑71 (MQ=255)
                                                      taataaAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGa                         <  1:5529710/63‑1 (MQ=255)
                                                      tAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            <  1:3734190/44‑1 (MQ=255)
                                                      tAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            <  1:260593/44‑1 (MQ=255)
                                                       aaTCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCAcgc                >  1:5371538/1‑71 (MQ=255)
                                                        aTCAAAGCCACGGTTTTGAGGTTGGCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCAcgcg               >  1:4726796/1‑71 (MQ=255)
                                                        aTCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCAcgcg               >  1:2584644/1‑71 (MQ=255)
                                                        aTCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCAcgcg               >  1:7598995/1‑71 (MQ=255)
                                                        aTCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCAcgcg               >  1:5044664/1‑71 (MQ=255)
                                                         tCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaaa                                             <  1:6294090/40‑1 (MQ=255)
                                                         tCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAt                                           >  1:38420/1‑42 (MQ=255)
                                                         tCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAt                                           >  1:4291067/1‑42 (MQ=255)
                                                         tCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgt              <  1:6924354/71‑1 (MQ=255)
                                                         tCAAAACCGTGGCTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            <  1:6894251/41‑1 (MQ=21)
                                                         tCAAAACCGTGGCTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            <  1:7124721/41‑1 (MQ=21)
                                                         tCAAAACCGTGGCTTTGAGGTTGCCATTCTTCCGCAGaaaa                                            <  1:3195172/41‑1 (MQ=21)
                                                         tCAAAACCGTGGCTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGTTGAAGTTGTCATGAt                        >  1:4855475/1‑61 (MQ=37)
                                                         tCAAAACCGTGGCTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGAt                        >  1:4393201/1‑61 (MQ=38)
                                                         tCAAAACCGTGGCTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGAt                        >  1:1492066/1‑61 (MQ=38)
                                                          cAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaa                                              <  1:2874409/38‑1 (MQ=255)
                                                          cAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGaa                                              <  1:5160233/38‑1 (MQ=255)
                                                          cAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtt             >  1:1694283/1‑71 (MQ=255)
                                                           aaaGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgttt            >  1:832500/1‑71 (MQ=255)
                                                            aaGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtttg           <  1:627765/71‑1 (MQ=255)
                                                            aaGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgttt            <  1:5826531/70‑1 (MQ=255)
                                                             aGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATg                          <  1:3372258/55‑1 (MQ=255)
                                                             aGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTc                    >  1:6570556/1‑61 (MQ=255)
                                                             aGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTc                    >  1:1042227/1‑61 (MQ=255)
                                                             aGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtttgt          <  1:3268177/71‑1 (MQ=255)
                                                             aGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtttgt          <  1:1877072/71‑1 (MQ=255)
                                                              gCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCAt                           <  1:3033778/53‑1 (MQ=255)
                                                              gCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCAcg                 <  1:7092227/63‑1 (MQ=255)
                                                               ccACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAAt                                           >  1:3422606/1‑36 (MQ=255)
                                                                   ggTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtt             <  1:856997/62‑1 (MQ=255)
                                                                   ggTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtt             <  1:3132980/62‑1 (MQ=255)
                                                                   ggTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCt    <  1:1106472/71‑1 (MQ=255)
                                                                    gTTTTGAGGTTGCCATTCTTCCGCAGAAAATgtgg                                       >  1:7001836/1‑35 (MQ=255)
                                                                    gTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGt                                 <  1:5619864/41‑1 (MQ=255)
                                                                    gTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGt                                 <  1:2919757/41‑1 (MQ=255)
                                                                    gTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCgtt             <  1:6871051/61‑1 (MQ=255)
                                                                    gTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCtt   <  1:3510760/71‑1 (MQ=255)
                                                                    gTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCtt   <  1:1325784/71‑1 (MQ=255)
                                                                     ttttGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCTTc  <  1:6235748/71‑1 (MQ=255)
                                                                     ttttGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCTTc  >  1:1780860/1‑71 (MQ=255)
                                                                     ttttGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCTTc  <  1:1773827/71‑1 (MQ=255)
                                                                     ttttGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCTTc  <  1:1449097/71‑1 (MQ=255)
                                                                     |                                                                      
AACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGCAGAAAATGTGGTGAAGTTGTCATGATAGTCACGCGTTTGTTTATCTTC  >  AP009048/1583514‑1583653

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: