Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,994,820 G→T 8.2% intergenic (+322/‑62) gabP_5 → / → CJ019_04690 GABA permease/IS5 family transposase ISPsp10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,994,8200GT 8.2% 131.9 / 4.8 49intergenic (+322/‑62)gabP_5/CJ019_04690GABA permease/IS5 family transposase ISPsp10
Reads supporting (aligned to +/- strand):  ref base G (34/11);  new base T (2/2);  total (36/13)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-02

AGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGATCTCTTCATCGAAA  >  ACB122/4994688‑4994945
                                                                                                                                    |                                                                                                                             
agTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGAc                                                                                                                             >  1:3636/1‑135 (MQ=255)
agTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACg                                                                                                                            >  7:340483/1‑136 (MQ=255)
     ttttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGcccc                                                                                                                        >  3:61841/1‑135 (MQ=255)
     ttttCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGcccc                                                                                                                        >  3:145504/1‑135 (MQ=255)
               gAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCt                                                                                                               >  6:49626/1‑134 (MQ=255)
               gAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCt                                                                                                               >  6:170351/1‑134 (MQ=255)
                  tCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGaaa                                                                                                           >  3:355087/1‑135 (MQ=255)
                       ccTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGCGCTGAAATTCTg                                                                                                      >  8:390841/2‑135 (MQ=255)
                           cTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAgg                                                                                                  <  3:182793/135‑1 (MQ=255)
                              cAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCTACGCCCCTCGGGAGCTg                                                                                                              >  1:374446/1‑120 (MQ=255)
                              cAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCTACGCCCCTCGGGAGCTg                                                                                                              >  1:379208/1‑120 (MQ=255)
                              cAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCTACGCCCCTCGGGAGCTg                                                                                                              <  2:374446/120‑1 (MQ=255)
                              cAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCTACGCCCCTCGGGAGCTg                                                                                                              <  2:379208/120‑1 (MQ=255)
                              cAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTtc                                                                                              >  8:257138/1‑136 (MQ=255)
                                        cATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCAt                                                                                      >  3:266763/1‑134 (MQ=255)
                                        cATCTCAAGCGACTGCGTGCTCTGGCTGTCTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCAt                                                                                      >  3:323949/1‑134 (MQ=255)
                                         aTCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCt                                                                                    >  2:53889/1‑135 (MQ=255)
                                                       cGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCGACGGAAAGGAATTc                                                                      >  1:88418/1‑135 (MQ=255)
                                                             ctGGCTGTGTTAACAGCTCCCATCCTTGCGCCACTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGAt                                                                <  2:38143/135‑1 (MQ=255)
                                                                gCTGTGTTAACAGCTCCCATCCTTGCGCCGCGCCCTAGGTGCTGTACGAAATCCCGATAAAGTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATg                                                               <  7:248881/133‑1 (MQ=255)
                                                                    tgtTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAg                                                         >  1:241086/1‑135 (MQ=255)
                                                                       tAACAGCTCCCATCCTTGCTCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCg                                                       >  3:243320/1‑134 (MQ=255)
                                                                       tAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCgg                                                      >  7:67919/1‑135 (MQ=255)
                                                                       tAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCg                                                       >  5:119361/1‑134 (MQ=255)
                                                                         aCAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTGCGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTAc                                                   >  7:348415/1‑136 (MQ=255)
                                                                         aCAGCTCCCATCCTTGCGCCGCGCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGt                                                     <  5:170351/134‑1 (MQ=255)
                                                                                 cATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGATAGGAATTCTGTGATGGCAAAGCGGTACGAActct                                            >  5:199759/1‑135 (MQ=255)
                                                                                 cATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAActct                                            >  1:232091/1‑135 (MQ=255)
                                                                                   tCCTTGCGCCGCTCCCTAGGTTCTGTACTAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCgg                                         <  8:264393/136‑1 (MQ=255)
                                                                                   tCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCg                                          >  2:219087/1‑135 (MQ=255)
                                                                                   tCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGAAGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCg                                          >  3:303148/1‑135 (MQ=255)
                                                                                     cTTGCGCCGATCCCTAGGTTCTGTACGAAATCACGAGAAACTGAATCGACGCCCCTCGCGAGCTGAAATTCTTGAGGGTTCTCCGCCATCTACGGAAACGAATTATGTTATGGCAAAGCGGTACGAACTCTCGGa                                        >  2:212808/1‑135 (MQ=17)
                                                                                         cgccgcTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGt                                     <  2:232091/134‑1 (MQ=255)
                                                                                           ccgcTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAg                                                                                                                 <  8:211321/56‑1 (MQ=255)
                                                                                           ccgcTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAg                                                                                                                 >  7:211321/1‑56 (MQ=255)
                                                                                            cgcTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGcc                                 <  4:61841/135‑1 (MQ=255)
                                                                                            cgcTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCt                                >  8:284363/1‑136 (MQ=255)
                                                                                                 ccTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGa                            >  3:331496/1‑135 (MQ=255)
                                                                                                 ccTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGa                            >  2:165593/1‑135 (MQ=255)
                                                                                                            tACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCg                 >  5:5810/1‑135 (MQ=40)
                                                                                                            tACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCg                 >  3:79782/1‑135 (MQ=40)
                                                                                                              cGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGACCTCTCGGACGTGGCCTGGGATGTGGTCTCCGAt               >  6:321018/1‑135 (MQ=255)
                                                                                                              cGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGAt               >  5:92034/1‑135 (MQ=40)
                                                                                                              cGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGAt               >  3:185742/1‑135 (MQ=40)
                                                                                                              cGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGAt               >  2:311026/1‑135 (MQ=40)
                                                                                                              cGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTCGCCTGGGATGTGGTCTCCGAt               >  3:372345/1‑135 (MQ=255)
                                                                                                                      cGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGATCTCTTCAt       <  4:266763/135‑1 (MQ=40)
                                                                                                                           aagtcaaTCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGATCTCTTCATCGaaa  <  1:24048/130‑1 (MQ=255)
                                                                                                                           aagtGAAGCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGATCTCTTCATCGaaa  <  3:304510/132‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
AGTTGTTTTCCGCAGGAGTCTGGACTCCTGCAGAAGCTGCCATCTCAAGCGACTGCGTGCTCTGGCTGTGTTAACAGCTCCCATCCTTGCGCCGCTCCCTAGGTTCTGTACGAAATCCCGAGAAACTGAATCGACGCCCCTCGGGAGCTGAAATTCTGTAGGGTTCTCCGCCATCTACGGAAAGGAATTCTGTGATGGCAAAGCGGTACGAACTCTCGGACGTGGCCTGGGATGTGGTCTCCGATCTCTTCATCGAAA  >  ACB122/4994688‑4994945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑