Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 1,070,248 +GA intergenic (+48/+94) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,070,2461.G91.7% 69.2 / 0.0 24intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (11/11);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
*AE0154511,070,2462.A91.3% 61.3 / ‑0.5 24intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (11/10);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACC  >  AE015451/1069960‑1070536
                                                                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                                                                    
gAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGC‑GAAGGAAGGCTGGCCAGGCAAAGGGCGGAT‑GGGTTTGGTTGGTAGGGGTTATGCGTTTAGTTTTCCGGCGCC‑‑tgt‑‑gggttcgagcgccc                                                                                                                                                                                                                                                                                        >  1:89629/1‑288 (MQ=255)
       aGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCgcgggg                                                                                                                                                                                                                                                                                 >  1:697066/1‑300 (MQ=255)
                            cGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATATCACAGGCgg                                                                                                                                                                                                                                                            >  2:134178/1‑300 (MQ=255)
                                                           ccttccacaaatcccgtcgtcACCCGCTGTCCAACTGTCTGTGCCTTGCTTTCATCTTGCTGATCCTGGTAGTTATGTACCTGACGCCGCCGATTCGCATCTAGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                                                                                                                                                                                                             <  2:241326/281‑1 (MQ=255)
                                                           ccTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCTCTCGATCTCACAGGCGATGAAAGCCTACAGCCATGACTCTGGCGGccc                                                                                                                                                                                                                             >  1:587012/1‑301 (MQ=255)
                                                                             ggCACCCGCTGACCAACTACCGGTGACTGGCGTTCATCGTGCTGATGCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                                                                                                                                                                                           <  2:232890/301‑1 (MQ=255)
                                                                              gctcccggtggccAACTACTTGCGCCTGGCCTTCAGCATGCTGCGCCTGGTAGTGATGTACCTGACGACGCCGATTCGTAGTGCGGTAGTGCGGACCCCGGCATTGATCGCTATGCTGTGGGTGGCTTTCAAGCGGAA‑GAAGGCGCGGCAGTCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGGGCGCCGCGCGGGCGGCGCTGGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGGCCTGACGCGTTCCCCTTCCAt                                                                                                                                                                                                          <  2:119842/290‑1 (MQ=255)
                                                                                                        tGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAGCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCCCAATGCTGCCTCCGCCAAgc                                                                                                                                                                                >  1:420569/1‑301 (MQ=255)
                                                                                                               cATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGcaca                                                                                                                                                                         >  1:466973/1‑301 (MQ=255)
                                                                                                                  cGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGTGCGGCGCTCGAGCTCACCGGCGATGAAAGCCTAAAAGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAGTCCCCCCAATGCTCCCTCCGTCCAGCACCCACAggg                                                                                                                                                                      >  1:654785/1‑299 (MQ=255)
                                                                                                                         gtcctggtagcgatgggactggggccgccGATTCGCATCTCGGTACTGCTGATCCCGGCATGGATCGCTGTGGTGTGGGTGGCTATCAAGCTGAA‑GAAGGCTCGCCAGGCGAAGGGCTGATCGGGTTT‑GTTCGTAAGAGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                                                                                                                                                               <  2:369867/278‑1 (MQ=255)
                                                                                                                                                     cgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTAc                                                                                                                                   >  1:84591/1‑301 (MQ=255)
                                                                                                                                                       gttcgcagatcGGTAAGGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGGGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACAt                                                                                                                                 <  1:147672/292‑1 (MQ=255)
                                                                                                                                                          cGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGa                                                                                                                              >  1:491454/1‑301 (MQ=255)
                                                                                                                                                            gaTCTCGGTAATGCTGCTCACGGCATGGGTCGCTGTGCTGTGGGGGGCTTTCAATCTGAA‑GAAGGCTCGCCATGCCGAGGGCTGATCGGGTTT‑GTTCGCAAGGGGTATACCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCGCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAAc                                                                                                                            <  2:194119/300‑1 (MQ=255)
                                                                                                                                                                tcGGTAATGCTGAACCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGGCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCa                                                                                                                        <  1:294222/301‑1 (MQ=255)
                                                                                                                                                                         cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCCCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCGGCCGGCCCTGACGCGTTCCCCTTCCAGGCCCCTCAATGCTGCCCCCGCCAAGCACGCACTGTGGGGTACAGCACCACCGTAGCGAGCACTGGCAGTACACAGAACAGCATCAATCGcc                                                                                                               >  1:288323/1‑300 (MQ=255)
                                                                                                                                                                          tGATCCCGGCAGGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTc                                                                                                              <  2:84554/301‑1 (MQ=255)
                                                                                                                                                                                   aaagtaTCGCTGTGGTTGGGGTGGCTTACAAGCTGAA‑GAAGGAAGGCCAGGCCAATGGATGATAGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTGGGCGCCGGAGA‑‑GA‑‑TCGGGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgcc                                                                                                     <  1:100320/296‑1 (MQ=255)
                                                                                                                                                                                    aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgcca                                                                                                    >  1:707186/1‑301 (MQ=255)
                                                                                                                                                                                                        tggCTTTCGAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GAGTTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACATGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGt                                                                                <  1:692891/301‑1 (MQ=255)
                                                                                                                                                                                                            tttgaaccTGAA‑GAAGGCCCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTaaaa                                                                            <  2:89591/294‑1 (MQ=255)
                                                                                                                                                                                                                                                               gTAAGGGGTATGCCGGTAGGTTTGCGGCGCC‑‑TGTGAGA‑‑GCGAGCGCCGCGCGGGCGGAGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGc                           <  2:466771/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gggTTTTCGGCGAC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAACTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATg          <  2:696767/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        ggCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCAAACGGCATCCCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAAcc  >  1:665388/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                                                                    
GAGCATCAGCATTACCCACCTCAAATTCCGCAAGGCCAAGCTGGCCGCCGGGATCACCCCCTTCTACAAAAGCCTGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAA‑GAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTT‑GTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGCCAACGGCATACCCATGGTAAAAGCCGCCACATGAAGGCCGTAGGCCAGGTACAGGCAAAGAAACACCAGGCCTTCTGCGCGAGTGATGCGGTAACC  >  AE015451/1069960‑1070536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑