Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,628,787 A→T 5.4% L4* (TTA→TAA)  lplA ← lipoate‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,628,7870AT 5.4% 376.4 / 10.3 129L4* (TTA→TAA) lplAlipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base A (48/74);  new base T (5/2);  total (53/76)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 137 total aligned reads displayed.

CGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGG  >  AP009048/4628721‑4628857
                                                                  |                                                                      
cGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAt                                                                      <  1:2725363/69‑1 (MQ=255)
  aaaaaTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgtg                                                                   >  1:360174/1‑70 (MQ=255)
  aaaaaTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTgg                                                                  <  1:334901/71‑1 (MQ=255)
   aaaaTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                                 >  1:4500375/1‑71 (MQ=255)
   aaaaTACACTCTTCCACCGCCAGGTTAAAACACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTTAtgtcc                                                                  <  1:2894330/70‑3 (MQ=255)
    aaaTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                                 >  1:7236868/1‑70 (MQ=255)
     aaTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACa                                                               >  1:5806767/1‑71 (MQ=255)
      aTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACAt                                                              <  1:6782399/71‑1 (MQ=255)
      aTACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACAt                                                              <  1:160231/71‑1 (MQ=255)
         cacTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTTATGTCCACATTAc                                                           <  1:7579810/71‑1 (MQ=255)
         cacTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAt                                                                      >  1:1320778/1‑60 (MQ=255)
              ttCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTaa                                                                       >  1:2574199/1‑54 (MQ=255)
              ttCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAttt                                                      >  1:2877287/1‑71 (MQ=255)
              ttCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAttt                                                      >  1:2168693/1‑71 (MQ=255)
              ttCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAttt                                                      >  1:878291/1‑71 (MQ=255)
               tCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAt                                                        <  1:5354878/68‑1 (MQ=255)
                ccACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAttt                                                      <  1:5040021/69‑1 (MQ=255)
                ccACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAtt                                                       <  1:6193956/68‑1 (MQ=255)
                     gCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCt                                               >  1:1206215/1‑71 (MQ=255)
                       cAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTg                                             <  1:3189530/71‑1 (MQ=255)
                       cAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTg                                             <  1:7288346/71‑1 (MQ=255)
                        aGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGa                                                                 <  1:65077/50‑1 (MQ=255)
                         ggTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTg                                             >  1:6216253/1‑69 (MQ=255)
                         ggTTAAACCACGGGGCGTAAGAGTCAGAGATGAGCAGGCGTAAtg                                                                     <  1:4941653/45‑1 (MQ=255)
                          gTTAATCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtag                                          <  1:7241158/71‑1 (MQ=255)
                          gTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAtt                                                       >  1:862802/1‑58 (MQ=255)
                           ttAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:1180728/71‑1 (MQ=255)
                           ttAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                           >  1:5577105/1‑69 (MQ=255)
                            tAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAt                                                        <  1:6594733/55‑1 (MQ=255)
                            tAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:4599959/70‑1 (MQ=255)
                               aCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTTATGTCCACATTACg                                                          >  1:149776/1‑50 (MQ=25)
                               aCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTTATGTCCACATTACg                                                          >  1:135510/1‑50 (MQ=25)
                               aCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTTATGTCCACATTACg                                                          >  1:443140/1‑50 (MQ=25)
                               aCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTTATGTCCACATTACg                                                          >  1:116862/1‑50 (MQ=25)
                               aCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATaa                                                            >  1:1352668/1‑48 (MQ=255)
                               aCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                              <  1:7020758/62‑1 (MQ=255)
                                ccACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCt                                                   <  1:7765850/56‑1 (MQ=255)
                                ccACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:442348/66‑1 (MQ=255)
                                   cGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCttt                                                 >  1:5905610/1‑55 (MQ=255)
                                    gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCa                                <  1:4859283/71‑1 (MQ=255)
                                    gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCa                                <  1:6833754/71‑1 (MQ=255)
                                    gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCa                                <  1:1570848/71‑1 (MQ=255)
                                     ggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTaa                                       >  1:7007333/1‑63 (MQ=255)
                                     ggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                                   <  1:28435/67‑1 (MQ=255)
                                      gTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCAct                              >  1:1013613/1‑71 (MQ=255)
                                       tCGTAAGAGTCTGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                                            <  1:3055081/56‑1 (MQ=255)
                                       tCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                                            <  1:7017034/56‑1 (MQ=255)
                                         gTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAAt                                      <  1:4308026/60‑1 (MQ=255)
                                         gTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtt                           >  1:7268925/1‑71 (MQ=255)
                                          tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTc                                                     <  1:3632860/44‑1 (MQ=255)
                                          tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTc                                                     <  1:5738135/44‑1 (MQ=255)
                                          tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttt                          >  1:3711910/1‑71 (MQ=255)
                                          tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttt                          >  1:4054585/1‑71 (MQ=255)
                                          tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCATTTCTTGTAGTAATGGGTCACTCtt                           <  1:6140530/70‑1 (MQ=255)
                                            agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                           >  1:5248469/1‑52 (MQ=255)
                                            agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                           >  1:5239195/1‑52 (MQ=255)
                                            agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt                        >  1:7906479/1‑71 (MQ=255)
                                            agagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt                         <  1:5109593/70‑1 (MQ=255)
                                             gagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttttt                       <  1:3275154/71‑1 (MQ=255)
                                             gagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttt                        >  1:4733739/1‑70 (MQ=255)
                                               gTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtct                     >  1:2629302/1‑71 (MQ=255)
                                                tCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                                <  1:5475361/43‑1 (MQ=255)
                                                tCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                                <  1:6764099/43‑1 (MQ=255)
                                                tCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:6022441/50‑1 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                              >  1:5712674/1‑43 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                              >  1:330133/1‑43 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCtt                                              >  1:2381618/1‑43 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                           >  1:2083740/1‑46 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                           >  1:7693732/1‑46 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttt                         <  1:2617210/64‑1 (MQ=255)
                                                  agagATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGc                  >  1:1667951/1‑71 (MQ=255)
                                                   gagaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCt                 <  1:4186680/71‑1 (MQ=255)
                                                    agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtt                           <  1:471384/60‑1 (MQ=255)
                                                    agaTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCtt                >  1:1078311/1‑71 (MQ=255)
                                                      aTGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                           <  1:4731837/42‑1 (MQ=255)
                                                       tGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttct             >  1:3828998/1‑71 (MQ=255)
                                                       tGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttct             <  1:6464177/71‑1 (MQ=255)
                                                       tGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttct             <  1:2624570/71‑1 (MQ=255)
                                                       tGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttct             <  1:3309712/71‑1 (MQ=255)
                                                         aGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtttttt                       <  1:2439227/59‑1 (MQ=255)
                                                           cAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTtcttcct         <  1:6388380/71‑1 (MQ=255)
                                                            aGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTtcttcctc        <  1:3664986/71‑1 (MQ=255)
                                                             ggCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTtcttcctct       <  1:4899671/71‑1 (MQ=255)
                                                              tcGTTATGTCCACATTACGATTTCCTTTCTTGTAGTAATGGGTCAct                              >  1:4517679/2‑47 (MQ=21)
                                                              gCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:6167475/36‑1 (MQ=255)
                                                              gCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:3764139/36‑1 (MQ=255)
                                                              gCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:7566921/36‑1 (MQ=255)
                                                              gCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:3246970/36‑1 (MQ=255)
                                                              gCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                         <  1:2317342/36‑1 (MQ=255)
                                                              gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTtcttcctctt      >  1:3786515/1‑71 (MQ=255)
                                                               cGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTttct             <  1:232349/63‑1 (MQ=255)
                                                                gTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCtt                           <  1:1933367/48‑1 (MQ=255)
                                                                 tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                                   <  1:6829180/39‑1 (MQ=255)
                                                                 tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                                   <  1:6792920/39‑1 (MQ=255)
                                                                 tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                                   <  1:3840524/39‑1 (MQ=255)
                                                                 tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATggg                                   <  1:7541120/39‑1 (MQ=255)
                                                                 tAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActc                             >  1:2132945/1‑45 (MQ=255)
                                                                  aaTGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCgg  >  1:6304201/1‑71 (MQ=255)
                                                                  aaTGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCgg  >  1:1064420/1‑71 (MQ=255)
                                                                  aaTGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCg   <  1:7628867/70‑1 (MQ=255)
                                                                  |                                                                      
CGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGG  >  AP009048/4628721‑4628857

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: