Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,169,169 C→A 12.5% Q69K (CAA→AAA)  CJ019_03921 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,169,1690CA12.5% 128.0 / 2.4 57Q69K (CAA→AAA) CJ019_03921hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/17);  new base A (2/6);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 5.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTAAAACTCGCA  >  ACB122/4169040‑4169295
                                                                                                                                 |                                                                                                                              
cGCTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAAtt                                                                                                                           >  8:112782/1‑135 (MQ=255)
cGCTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAAtt                                                                                                                           >  2:218419/1‑135 (MQ=255)
cGCTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGAAGGCCGCCAAAtt                                                                                                                           >  1:25053/1‑135 (MQ=255)
 gCTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAACGCCTAACGCATTATTTTTATGCAGGCCGACAAAtta                                                                                                                          >  5:71459/1‑134 (MQ=255)
  cTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAAtttt                                                                                                                         >  7:289457/1‑135 (MQ=255)
  cTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAAtttt                                                                                                                         >  7:56824/1‑135 (MQ=255)
  cTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCAGTATTTTTATGCAGGCCGCCAAAtttt                                                                                                                         >  2:333487/1‑135 (MQ=255)
  cTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCAATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAAtttt                                                                                                                         >  7:152161/1‑135 (MQ=255)
      cGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTga                                                                                                                      <  8:264338/134‑1 (MQ=255)
          aGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAAt                                                                                                                 <  7:50230/135‑1 (MQ=255)
                 ccTCCATTCTTGTGGCGAAACTTGGTTAAGCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTaaaaa                                                                                                          <  4:26258/135‑1 (MQ=255)
                 ccTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTaaaaa                                                                                                          <  3:4272/135‑1 (MQ=255)
                      aTTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCaaa                                                                                                     >  7:191642/1‑135 (MQ=255)
                        tCTTGTGGCTAAACTTGTATAATCAATTCTATAGCACAAATAATATTATCTATATGCAACCAACTATCCAAAGAAAGCATAAATCATTATTTTTATGCAGTCCTAAAAATTTTTGAAGAATTTAAAAATCAAAg                                                                                                    >  6:339950/1‑134 (MQ=255)
                             tgGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTaaaa                                                                                               <  1:215129/134‑1 (MQ=255)
                              ggCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAt                                                                                             <  1:218419/135‑1 (MQ=255)
                              ggCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAt                                                                                             <  6:83171/135‑1 (MQ=255)
                              ggCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAt                                                                                             <  2:162368/135‑1 (MQ=255)
                                 gtaaCTTGGATTATGAATGCTATAGCAGTAATAATATTATCGATATGCAACCTACTATGCTAAGAAAGCCTTACGCATTATTTTTATGCAGGCAGCCTTTTTTTTGAAGAATTTAAATATCAATTTAAAAAtctc                                                                                          <  2:84881/133‑1 (MQ=255)
                                 gAAACTTGGATAATCAATGCTATATCACAAATAATATTTTCGATATGCAACCAACTATCCAAAGAAAGCCTAAAGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAtctc                                                                                          <  8:172744/135‑1 (MQ=255)
                                 gAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAtctc                                                                                          <  6:328674/135‑1 (MQ=255)
                                  aaaCTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAAACAAATATACAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTTAAGAATTTAAAAATCAAAGTAAAAAtctc                                                                                          <  3:176224/134‑1 (MQ=255)
                                   atcTTGGATAATGAATGCTATAGCACAAATAATATTATCGATATGCCACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATctccc                                                                                        <  6:65618/133‑1 (MQ=255)
                                   aaCTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATctccc                                                                                        <  8:249921/135‑1 (MQ=255)
                                                   gCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAAAgacgac                                                                        >  7:379321/1‑135 (MQ=255)
                                                       taGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTggg                                                                    >  6:116962/1‑135 (MQ=255)
                                                          cacaAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGa                                                                 >  3:49786/1‑135 (MQ=255)
                                                          cacaAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGa                                                                 >  1:281700/1‑135 (MQ=255)
                                                               ataataTTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACttt                                                            >  7:199319/1‑135 (MQ=255)
                                                                         caatatGCAACCAACAATCCAAACAAAGCAAAACGAATTATTTTTATGCAGGCCGCAAAATTTTTTAATAATTAAAAAATAAAACAAAAAATCTCCCTCCAAGTAGACACGACTGGGGGAACTTTATACTTTGcc                                                  <  8:237063/133‑1 (MQ=255)
                                                                          gATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACTACTGGGGGAACTTTATAATTTGCCa                                                 >  6:40975/1‑135 (MQ=255)
                                                                          gATATGCAACCAACTATCAAAAGAAATAATAACGAATTATTATTATGAAGGCAGCAAAATTTTTGAAGAATTTAAAAATCAAATTAAAAATCTCCCTCAAAGTAGAGACGACTGGGGGAACATTAATCTTTGcaa                                                 >  1:342877/1‑133 (MQ=255)
                                                                           atatGCAACCAACTATCCAAAGAAAGCCTAACTCATTATTTTTATGCAGGCCGCCAAATTTTTGATGAATTTAAAAATCAAAGTAAATATCTCACTCCAAGTAGAGACTACTGGGGGAAGATTATACTTTGCCa                                                 >  8:351421/1‑134 (MQ=255)
                                                                           atatGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCa                                                 >  7:123735/1‑134 (MQ=255)
                                                                              tGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTa                                              >  4:228906/1‑134 (MQ=255)
                                                                               ccaaacaacTATCCAAACAAAGCATAACGAATTATTTTTATGCAGGCCGAAAAATTTTTGAATAATTAAAAAATAAAAGAAAAAATCACCCACCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATg                                            <  4:49786/134‑1 (MQ=255)
                                                                                caaccaacTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATg                                            >  5:379226/1‑134 (MQ=255)
                                                                                caaccaacTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATg                                            >  2:248829/1‑134 (MQ=255)
                                                                                      acTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACa                                     >  6:374947/1‑135 (MQ=255)
                                                                                         aTCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAtctc                                                                                          >  1:200458/1‑79 (MQ=255)
                                                                                         aTCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAAtctc                                                                                          <  2:200458/79‑1 (MQ=255)
                                                                                         aTCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACAtc                                   >  6:168120/1‑134 (MQ=255)
                                                                                             aaagaaagCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAAc                              >  8:104863/1‑135 (MQ=255)
                                                                                                       taactaatttttTTTAGGCAGGCCGCAAAATTTTTGAATAATTTAAAAAAAAAAGTAAAAATCTCCATCCCTGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAtt                    <  8:182721/129‑1 (MQ=255)
                                                                                                       tAACTCTTTATTTTTATGCAGGCCGAAAAATTTTTTAATAATTTAAAAATCAAAGTAAAAATCTCCCACCAATTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAtt                    <  2:25654/135‑1 (MQ=255)
                                                                                                       tAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAt                     >  4:34862/1‑134 (MQ=255)
                                                                                                       tAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAt                     >  7:90171/1‑134 (MQ=255)
                                                                                                       tAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAt                     >  2:230997/1‑134 (MQ=255)
                                                                                                       tAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAt                     >  2:229376/1‑134 (MQ=255)
                                                                                                       tAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACAt                     >  6:123024/1‑134 (MQ=255)
                                                                                                                 tttttATGCAGGCCTCAAAATTTTTTAATAATATAAAAATAAAATAAAAAATCACCAACCAAGTCCACACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTa           <  3:228906/134‑1 (MQ=255)
                                                                                                                  ttttATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTaaa         >  8:289511/1‑135 (MQ=255)
                                                                                                                   tttATCCAGCCCTCCAAATTTTTTAATAATTAAAAAAAAAAATAAAAAAGCTCCCACCAAGTAGCCACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTaaa         <  1:333487/134‑1 (MQ=255)
                                                                                                                    ttAGGCAGGCCCCAAAATTTTTGAAGAATTAAAAAATCAAACAAAAAATCACCCTCCAAGTAGCCACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTaaaa        <  2:281700/134‑1 (MQ=255)
                                                                                                                      ctccAGGCCCCCAAATTTTTAAAATATATAAAACTCAAAAAAAAAACCCCCGTCCATTACTAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTAAAACt      <  4:221037/131‑1 (MQ=255)
                                                                                                                          aGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTAAAACTCGCa  >  4:366388/1‑134 (MQ=255)
                                                                                                                          aGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTACCTCCAAGTATAGACGACGGGGGGCACTTTTTACTTTTCCACTTTGAAAAACATAAGCACCGTACACCTTAAACGGATAAAACTCGCa  >  1:174227/1‑134 (MQ=255)
                                                                                                                                 |                                                                                                                              
CGCTACCGAAAGATGAGCCTCCATTCTTGTGGCGAAACTTGGATAATCAATGCTATAGCACAAATAATATTATCGATATGCAACCAACTATCCAAAGAAAGCCTAACGCATTATTTTTATGCAGGCCGCCAAATTTTTGAAGAATTTAAAAATCAAAGTAAAAATCTCCCTCCAAGTAGAGACGACTGGGGGAACTTTATACTTTGCCACTATGAAAAACATCTCAACCTTACACATTCAACGCTTAAAACTCGCA  >  ACB122/4169040‑4169295

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