Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,002,486 T→C 8.1% M1T (ATG→ACG) † CJ019_04695 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,002,4860TC 8.1% 140.3 / 2.9 49M1T (ATG→ACG) †CJ019_04695hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/17);  new base C (2/2);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

TCTGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGGTGATGACCCAAGCTGTCAAACAGACC  >  ACB122/5002352‑5002619
                                                                                                                                      |                                                                                                                                     
tCTGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCAt                                                                                                                                       <  8:161260/135‑1 (MQ=255)
tCTGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCAt                                                                                                                                       <  4:252442/135‑1 (MQ=255)
  tGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGCGGAGAAATACGCATGt                                                                                                                                     >  4:31003/1‑135 (MQ=255)
  tGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGt                                                                                                                                     >  6:98187/1‑135 (MQ=255)
  tGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGt                                                                                                                                     >  8:64600/1‑135 (MQ=255)
     ttCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGa                                                                                                                                  <  7:287556/135‑1 (MQ=255)
     ttCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGa                                                                                                                                  <  1:266567/135‑1 (MQ=255)
     ttCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGa                                                                                                                                  <  4:163904/135‑1 (MQ=255)
      tCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGa                                                                                                                                  <  5:98187/134‑1 (MQ=255)
          ccAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCa                                                                                                                             <  8:209436/135‑1 (MQ=255)
            agcagAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGc                                                                                                                           >  2:113719/1‑135 (MQ=255)
            agcagAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGAACCCCAGAGGAGAAATACGCCTCTAGACTCCAGc                                                                                                                           >  6:195392/1‑135 (MQ=255)
                gagaTTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTg                                                                                                                       <  6:131702/135‑1 (MQ=255)
                     gacaggggTATGGTCGTAATATACTTTTTTTGCTAGATTATATTTTCTTGCCGCCCTTGACGTGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTcc                                                                                                                   <  6:164384/129‑1 (MQ=255)
                      cAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGc                                                                                                                 >  3:146021/1‑135 (MQ=255)
                           ggTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGa                                                                                                            >  6:292450/1‑135 (MQ=255)
                           ggTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCGGa                                                                                                            >  6:51625/1‑135 (MQ=255)
                                  tCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTg                                                                                                         >  3:145361/1‑131 (MQ=255)
                                  tCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTg                                                                                                         <  4:145361/131‑1 (MQ=255)
                                  tCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGAt                                                                                                     >  6:280114/1‑135 (MQ=255)
                                        tataCTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATctgtct                                                                                               >  5:164384/1‑130 (MQ=255)
                                                 tttGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTACTCTTACCACAGAGCAGAAATACGCAGCACGACTCAAGCTCTGTTCCGCACTGCGTGCGAACCTGTGAAggtggac                                                                                      >  4:236239/1‑133 (MQ=255)
                                                  ttGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCg                                                                                     >  8:240076/1‑135 (MQ=255)
                                                  ttGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCg                                                                                     >  2:173328/1‑135 (MQ=255)
                                                  ttGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCACGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCg                                                                                     >  1:239775/1‑135 (MQ=255)
                                                  ttGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCACGTCGACTCAAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCg                                                                                     >  1:235873/1‑135 (MQ=255)
                                                  ttGCTAGATAATATTTTCTTGCCGCCCTTGACGAGACCCATTTTTGCGTCTTCATTTACTCGTTCCCCAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCg                                                                                     >  8:359/1‑135 (MQ=255)
                                                         ataataTGTTCTTGCCGCCCGTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGATAAAAACGCATGTCGACGCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGcc                                                                              <  2:75377/135‑1 (MQ=255)
                                                           aataTTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCTTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCACGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCg                                                                            <  2:235873/135‑1 (MQ=255)
                                                           aataTTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCACGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCg                                                                            <  2:239775/135‑1 (MQ=255)
                                                              atTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGAtt                                                                         <  5:280114/135‑1 (MQ=255)
                                                              atTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGAtt                                                                         <  1:113719/135‑1 (MQ=255)
                                                              atTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGAtt                                                                         <  3:188965/135‑1 (MQ=255)
                                                              atTTTCTTGCCGCCCTGGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACCGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGAtt                                                                         <  5:5504/135‑1 (MQ=255)
                                                                     tGCCGCCCTTGACGAGCCGCATTTTTGCGTCTTCATTTACTCGTAACACAGAGGAGAAATAAGCATGTCGACTACAGCTCTGTTCCGCCCTGAGTGAGATCCTGTGACGGTGGTCGGTAGGACCGATTACGTCag                                                                  >  3:198734/1‑134 (MQ=255)
                                                                     tGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCAt                                                                  >  3:111814/1‑135 (MQ=255)
                                                                          cccTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACa                                                             >  5:34173/1‑135 (MQ=255)
                                                                            cTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGAtt                                                                         <  6:109715/121‑1 (MQ=255)
                                                                            cTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGAtt                                                                         >  5:109715/1‑121 (MQ=255)
                                                                                       ccATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCAAGCTCTGTTCCGACCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACATCTGCCTGCAGtt                                                >  6:161394/1‑134 (MQ=255)
                                                                                        cATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTg                                                >  4:265169/1‑134 (MQ=255)
                                                                                             ttGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGa                                          >  1:200164/1‑135 (MQ=255)
                                                                                                       aTTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGa                                                                                                            >  7:188154/1‑59 (MQ=255)
                                                                                                       aTTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGa                                                                                                            <  8:188154/59‑1 (MQ=255)
                                                                                                                 aCCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGGtgatga                      >  1:81754/1‑135 (MQ=255)
                                                                                                                   cacaGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGt                                                                                   >  5:128830/1‑72 (MQ=255)
                                                                                                                   cacaGAGGAGAAATACGCATGTCGACTCCAGCTCTCTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGt                                                                                   <  6:128830/72‑1 (MQ=255)
                                                                                                                         ggagAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGGTGATGACCCAAGCt              >  3:196277/1‑135 (MQ=255)
                                                                                                                           agaAATAAGCATGTCGACTCCAGCTCTGTTCAGACCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATATACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGCTGATGACCCAAGCtc             >  8:65039/1‑133 (MQ=255)
                                                                                                                                    cATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGGTGATGACCCAAGCTGTCAAACAGAc   <  8:147665/135‑1 (MQ=255)
                                                                                                                                     aTGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGGTGATGACCCAAGCTGTCAAACAGAcc  >  4:17494/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TCTGTTTCGACCAGCAGAGATTCAAAGGGTATGCTCGTAATATACTGATTTTGCTAGATAATATTTTCTTGCCGCCCTTGACGAGCCCCATTTTTGCGTCTTCATTTCCTCGTACCACAGAGGAGAAATACGCATGTCGACTCCAGCTCTGTTCCGCCCTGAGTGCGATCCTGTGACGGTGGTCGGTAGGCCCGATTACGTCATCTACAGCTGCCTGCAGTGGCAAGACCAAAATCTCCGGGTGATGACCCAAGCTGTCAAACAGACC  >  ACB122/5002352‑5002619

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: