Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,357,973 A→T 8.2% intergenic (+50/+76) CJ019_04103 → / ← CJ019_04104 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,357,9730AT 8.2% 126.8 / 3.4 49intergenic (+50/+76)CJ019_04103/CJ019_04104hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base T (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

ACAAGCCCTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTCAGCGGCGACACGTCGGG  >  ACB122/4357839‑4358107
                                                                                                                                      |                                                                                                                                      
acaAGCCCTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCACACTCTGTCGGATAGCTGGAAGTGTGGCGa                                                                                                                                        >  3:7370/1‑135 (MQ=255)
    gCCCTGGAAAACACGGTGCGCTATGTCAACGTCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCGGGAAGTGTGGCGACGCa                                                                                                                                    <  4:170986/135‑1 (MQ=255)
    gCCCTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    <  1:200243/135‑1 (MQ=255)
    gCCCTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    <  8:98454/135‑1 (MQ=255)
    gCCCTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    <  5:87672/135‑1 (MQ=255)
    ccccTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    <  1:143094/134‑1 (MQ=255)
          gAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCAt                                                                                                                               >  8:75181/1‑134 (MQ=255)
             aaCACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCa                                                                                                                           >  3:223274/1‑135 (MQ=255)
               cacGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACg                                                                                                                         <  4:97301/135‑1 (MQ=255)
               cacGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACg                                                                                                                         <  4:173914/135‑1 (MQ=255)
                acGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGc                                                                                                                        <  8:47928/135‑1 (MQ=255)
                       gcTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGt                                                                                                                 >  5:109915/1‑135 (MQ=255)
                              cAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCt                                                                                                          <  7:222217/135‑1 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGTCGCAGGCATTTCACGCACTGCGTCTTCGCtgatga                                                                                                     >  7:92858/1‑135 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGTCGCAGGCATTTCACGCACTGCGTCTTCGCtgatg                                                                                                      >  7:92861/1‑134 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    <  3:150229/104‑1 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    >  4:150229/1‑104 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    >  4:145354/1‑104 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCa                                                                                                                                    <  3:145354/104‑1 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCtggtg                                                                                                      >  5:152155/1‑134 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTAACGCACTGCCTCTTCGCggatg                                                                                                      >  3:186443/1‑134 (MQ=255)
                                   gCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCAGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACCCACTGCGTCTTCGCtgatg                                                                                                      >  3:139512/1‑134 (MQ=255)
                                    cTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGTCGCAGGCATTTCACGCACTGCGTCTTCGCTGATGAc                                                                                                    <  8:92858/135‑1 (MQ=255)
                                    cTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGTCGCAGGCATTTCACGCACTGCGTCTTCGCTGATGAc                                                                                                    <  8:92861/135‑1 (MQ=255)
                                          gatgaAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTAGCTGCTgtaagtaa                                                                                               >  5:272583/1‑134 (MQ=255)
                                               aaaaaTCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTcgcc                                                                                         <  3:240151/135‑1 (MQ=255)
                                                 aattCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTcgccgc                                                                                       <  8:59082/132‑1 (MQ=255)
                                                 aaaTCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTAACGCACTGCGTCTTCGCTGATGACAGTAAGTcgccgc                                                                                       >  7:175197/1‑135 (MQ=255)
                                                                 gcgcTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGCAGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCTgcgcgg                                                                       >  5:24018/1‑133 (MQ=255)
                                                                         cAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGGGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTc                                                               <  3:85890/135‑1 (MQ=255)
                                                                                    aCGGGCCTGGGTGAGTTTGCCCACTCTGTCGGATAGGTGGACGTGTGGCTAGGCAGGCTTTTCCCGATCTGCGTGTTCGATGATGAACGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCg                                                     <  4:186443/134‑1 (MQ=255)
                                                                                    aCGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGTATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGa                                                    >  7:18159/1‑135 (MQ=255)
                                                                                          cAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGc                                              <  7:75181/135‑1 (MQ=255)
                                                                                           aGGCTGATGATGCCCACTCTGGCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGcc                                             <  4:202325/135‑1 (MQ=255)
                                                                                              cTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGCCGCAGGCATTTAACGCACTGCGGCTTGGCGGATGACAGTAAGTCGCCGCTCGCCCACCCGCGCGGCGCGCCGCCCCGGTTTAGGGGCCGCCAgg                                          >  5:302906/1‑134 (MQ=255)
                                                                                               tgatgaTGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGaa                                         <  4:223274/135‑1 (MQ=255)
                                                                                                    atgCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTAACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCCGCCAGAAAGCCt                                    >  4:183207/1‑135 (MQ=255)
                                                                                                     tgCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCGATTTCGAGGCAGCCAGCAAGccag                                   >  2:37301/1‑133 (MQ=255)
                                                                                                                   gATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTa                                                                                                >  6:157026/1‑60 (MQ=255)
                                                                                                                   gATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTa                                                                                                <  5:157026/60‑1 (MQ=255)
                                                                                                                           gAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCGCCGCGATTTCGAGGCCGCCAGAAAGCCTGCCGCCGGCGCACCGTCAgcggc             >  3:29642/1‑135 (MQ=255)
                                                                                                                           gAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTAGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGGTATTTCGAGGCAGCCAGAAACCCAGCCGCAGGCCCCCCTTCAccggc             >  7:242357/1‑135 (MQ=255)
                                                                                                                            aaGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTCAgcgg              <  7:87061/133‑1 (MQ=255)
                                                                                                                             aGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGCGGCAGCCAGAAAGCCAGCCGCAGGCGCACCTTCAgcggcgc           >  7:139339/1‑134 (MQ=255)
                                                                                                                             aGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCGCCGCGATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCCCCTTCAgcggcgc           >  3:84174/1‑134 (MQ=255)
                                                                                                                             aGTGTGGCGACGCAGGCATTTAAAGCACTGCGTGTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCACCGCGTATTCGCGGCCGCCAGAAAGCCTTCCGCACGCCCACCTTCAgcggcca           >  6:109194/1‑133 (MQ=255)
                                                                                                                             aGTGTGGCGACGCAGGCAATTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCACCGGCCCTGCGCGCCGCCGCGATTTAGAGGCAGCCACAAAGCCAGCCGCACGCGCCCCGTCAgccgcgc           >  2:172354/1‑134 (MQ=255)
                                                                                                                               tgtgGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTCAGCGGCGaca         <  5:127734/135‑1 (MQ=255)
                                                                                                                                     gACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTCAGCGGCGACACGTCgg   >  7:86482/1‑135 (MQ=255)
                                                                                                                                     gACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCGCCGCGATTTAGAGGCAGCCAGAAAGCCTGCCGCAGGCGGACCTTCAGCGGCGACACGTCgg   >  5:155860/1‑135 (MQ=255)
                                                                                                                                      aCGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCGATTTCGAGGCCGCCAGAAAGCCTGCCGCAGGCGCCCCTTCCGCGGCGCCCCGTCggg  >  3:64355/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                      
ACAAGCCCTGGAAAACACGGTGCGCTATGTCAACGGCTCCTTGATGAAAAAATCCAAAAGTGAGGGCGCTGGCCAGCGAGGATGACGGGCCAGGCTGATGATGCCCACTCTGTCGGATAGCTGGAAGTGTGGCGACGCAGGCATTTCACGCACTGCGTCTTCGCTGATGACAGTAAGTCGCCGCTCGACCTCCGGCGCTGCGCGCCTCCGCTATTTCGAGGCAGCCAGAAAGCCTGCCGCAGGCGCACCTTCAGCGGCGACACGTCGGG  >  ACB122/4357839‑4358107

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: ‑