Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,371,271 Δ2 bp 100% intergenic (‑6/+299) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,371,2710A.98.2% 224.4 / ‑1.2 56intergenic (‑6/+300)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (25/30);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*AP0090484,371,2720A.98.2% 222.4 / ‑1.2 56intergenic (‑7/+299)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (25/30);  minor base G (0/1);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTA  >  AP009048/4371203‑4371334
                                                                    ||                                                              
gTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                               <  1:5010206/71‑1 (MQ=255)
 ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                              <  1:2965827/71‑2 (MQ=255)
 ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                              <  1:2430263/71‑2 (MQ=255)
 ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                               >  1:4218641/1‑70 (MQ=255)
    cGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                               <  1:2425780/67‑1 (MQ=255)
         cagcagCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa                                                                <  1:485685/61‑1 (MQ=255)
         cagcagCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                     <  1:3162383/70‑1 (MQ=255)
          agcagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGa                                                    <  1:1500726/70‑1 (MQ=255)
          agcagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                   <  1:222372/71‑1 (MQ=255)
             agcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc                                                <  1:703562/71‑1 (MQ=255)
              gcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCCACCTGCATGGCAG‑‑AAATAGCGGCGATAAcc                                               <  1:2830235/71‑1 (MQ=255)
               cTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCg                                              >  1:4917176/1‑71 (MQ=255)
                ttttCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGAGGCGATAACCGc                                             >  1:2510783/1‑71 (MQ=255)
                 tttCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGc                                             >  1:2001963/1‑70 (MQ=255)
                       tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGAt                                       >  1:4896991/1‑70 (MQ=255)
                       tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATa                                      <  1:1752987/71‑1 (MQ=255)
                        ctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                     >  1:2076693/1‑71 (MQ=255)
                         tgatgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                     >  1:83617/1‑70 (MQ=255)
                           atgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg                                   >  1:3757946/1‑70 (MQ=255)
                           atgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCAATGATAATg                                   >  1:3280303/1‑70 (MQ=255)
                            tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga                                 <  1:301849/71‑1 (MQ=255)
                            tgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCCATAACCGCCATGATAATGga                                 <  1:3694452/71‑1 (MQ=255)
                             gaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga                                 >  1:237792/1‑70 (MQ=255)
                              aaCCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                     <  1:5359522/49‑1 (MQ=255)
                              aaCCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag                                <  1:930103/70‑1 (MQ=255)
                               aCCAGATAGTCCAGACCGCCAGCAACCTGCATGG‑‑GAGAAATAGCGGCGATAACCGCCATGATAATGGAGAt                              <  1:1596109/71‑1 (MQ=255)
                                 cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGa                                                                 <  1:569851/36‑1 (MQ=255)
                                 cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa                            <  1:2608270/71‑1 (MQ=255)
                                 cAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa                            <  1:2224253/71‑1 (MQ=255)
                                  aGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa                                                                >  1:3033125/1‑36 (MQ=255)
                                  aGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGa                            >  1:1445194/1‑70 (MQ=255)
                                   gATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa                                        >  1:478540/1‑57 (MQ=255)
                                    aTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                     <  1:4957307/59‑1 (MQ=255)
                                    aTAGTCCAGACAGCCAGCAACCTGCATGGCAGaaa                                                               >  1:4436122/1‑35 (MQ=255)
                                      aGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa                                        >  1:4157694/1‑54 (MQ=255)
                                      aGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                     >  1:4434280/1‑57 (MQ=255)
                                        tCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCg                       <  1:2624536/69‑1 (MQ=255)
                                          cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAAc                    <  1:767584/70‑1 (MQ=255)
                                           aGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGa                      <  1:4028718/67‑1 (MQ=255)
                                            gACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                  >  1:4161092/1‑70 (MQ=255)
                                            gACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                  >  1:1030668/1‑70 (MQ=255)
                                             aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa                                        >  1:199372/1‑47 (MQ=38)
                                              ccGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                              >  1:64943/1‑56 (MQ=255)
                                               cGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                              <  1:2575728/55‑1 (MQ=38)
                                                gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt             <  1:4641476/71‑1 (MQ=255)
                                                 ccAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt             >  1:2604381/1‑70 (MQ=255)
                                                  cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                              <  1:5344189/52‑1 (MQ=38)
                                                  cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                              <  1:3090219/52‑1 (MQ=38)
                                                  cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                              <  1:5078904/52‑1 (MQ=38)
                                                  cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGtt            <  1:173872/70‑1 (MQ=255)
                                                   aGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACggg                 >  1:4197969/1‑64 (MQ=255)
                                                   aGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                  >  1:3740688/1‑63 (MQ=255)
                                                    gCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTAcc         >  1:4752602/1‑71 (MQ=255)
                                                    gCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTAcc         >  1:4912293/1‑71 (MQ=255)
                                                      aaCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag                                <  1:5325712/46‑1 (MQ=37)
                                                        ccTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt     >  1:2496330/1‑71 (MQ=255)
                                                        ccTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt     <  1:4707231/71‑1 (MQ=255)
                                                         cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg                                   <  1:1859394/40‑1 (MQ=37)
                                                          tGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGa                >  1:2746987/1‑58 (MQ=255)
                                                          tGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGtt    >  1:3566527/1‑70 (MQ=255)
                                                            cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt   <  1:4355146/69‑1 (MQ=255)
                                                            cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt   <  1:2530449/69‑1 (MQ=255)
                                                            cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa  <  1:522370/70‑1 (MQ=255)
                                                            cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa  <  1:206374/70‑1 (MQ=255)
                                                            cATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa  <  1:1184676/70‑1 (MQ=255)
                                                                    ||                                                              
GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTA  >  AP009048/4371203‑4371334

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: