Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,986,622 T→G 38.9% D31A (GAT→GCT)  CJ019_03763 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,986,6220TG38.9% 11.9 / 6.2 37D31A (GAT→GCT) CJ019_03763hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/12);  new base G (9/6);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAG  >  ACB122/3986496‑3986749
                                                                                                                              |                                                                                                                               
gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGAGCCGGcgcg                                                                                                                         >  7:52244/1‑135 (MQ=255)
     cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGGGCGGGGgccggcgccggccc                                                                                                                    >  7:269135/1‑133 (MQ=255)
      aTCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCggcgcgggcgcg                                                                                                                   >  5:44807/1‑135 (MQ=255)
      aTCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCggcgcgggcgcg                                                                                                                   >  3:65773/1‑135 (MQ=255)
         cgCACGTCAGCCAGCAAGTGTTCCGGGTACACGGAATTGGGGGCCTGCATCTGGGCCCGTTGGTCCGGGGAGAACCCGGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGt                                                                                                                  <  4:236227/132‑1 (MQ=255)
            cccTCAGCAAGCAACTGATCCGGGTACATGGCAAGGGGGGCCTCCATCTGGTCCAGTTGCGCCTGGGAGAACCCGGGGTCCTTGGCTACGAGCGGGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGTCTGCg                                                                                                             <  6:90973/134‑1 (MQ=255)
                 aGCCAGCAACTGAGCAGGGTACAAGGCAATGGGGGCCAGCATCTGGGCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                        <  8:93094/135‑1 (MQ=255)
                    cAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGGGGCGCCGGCGCGGGCCCGGCGGcggccggccc                                                                                                     >  5:209087/1‑133 (MQ=255)
                       caccTGATCCGGTTACACGCAAAGGGGGGCCAGCATCGGGTCCAGTGGCACCTGGGACACCACCGGGTCCTTGGCTACGAGCGGGGCCGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCAccag                                                                                                  <  4:245715/132‑1 (MQ=255)
                        aaCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCCGAGCCGGCGCGGGGCCGTCGGCGGCCCGccccccc                                                                                                 >  6:123372/1‑132 (MQ=255)
                          cTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCGGCCCGGGCCCGGCGGCGGCCGGGAccacg                                                                                                 >  3:56803/1‑131 (MQ=255)
                                          gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCAGGCCCGGGCCCGTCTGCCGCCGGCCCCACACCCCCGGCAgcgcgcccg                                                                               >  7:63490/1‑130 (MQ=255)
                                          gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCCGCGCCGGCCCCGCAGCGGCCGGCACCACGCCCCGGGCGGCGCGCAgg                                                                               >  7:22348/1‑135 (MQ=255)
                                               gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGGGGGTGCGGCGCGGGCGGGGCGGCGGGCGGCCCCAGCCCCCGGGCCGCGCGCCGGcgccg                                                                          >  5:298840/1‑131 (MQ=255)
                                               gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGCGGCGGCCCGGGGCCGGCTGCGGCCGGCCCCCCGCCCCGGGGCGCGCGCCGGCaaac                                                                          >  1:167676/1‑134 (MQ=255)
                                               gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTTCGGCCTGCCCCAGGCCAGGGGCAGCGCTCAGGCaaaa                                                                          >  4:28074/1‑134 (MQ=255)
                                               gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCCGCGCGGGGACCTCTGCGGCCGGCACCACGCCAGGGGCAGCGCGCAGGCaccg                                                                          >  1:164841/1‑132 (MQ=255)
                                                 gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCCGCGCCGGGCCGGCGGCGGGCGGCACCACGCCCCGGGCGGCGCGCCGGCCAAGCa                                                                        >  3:275476/1‑135 (MQ=255)
                                                      aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGGGCCGGCGCGGGGGCGGCGGCGGGCGGCGCCCGGCCCCGGGCCGCGCGCCGGCCCAGCacgggc                                                                   >  1:13850/1‑130 (MQ=255)
                                                           cTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGCTCCCGCGCGGGCCCGTCGGCGGCCGGGCGCAGGCCCGGGGCCCCGCTCCGGCCAACCAGTGCCggccc                                                              >  2:291867/1‑131 (MQ=255)
                                                                ccAGTTGCGCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAg                                                         <  8:120313/135‑1 (MQ=255)
                                                                     tGCTCCTGTCACAACACCGGGTCGTTGGCTACGCGGGGGGCAGGGCAGGGGGCGGGAGCGGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCaa                                                      <  1:51440/133‑1 (MQ=255)
                                                                       cTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACaa                                                  >  2:269888/1‑135 (MQ=255)
                                                                        tCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAAt                                                 <  6:229373/135‑1 (MQ=255)
                                                                        tCCTGGGAGAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAAt                                                 <  8:74884/135‑1 (MQ=255)
                                                                             ggccaaCACCGGGTCCTTGGCTACGAGGGGGGCAGCGGACGGCGCCGTATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                            <  7:147115/131‑1 (MQ=255)
                                                                              gagaACACCGGGTCCTGGGATACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                            <  6:209087/134‑1 (MQ=255)
                                                                                           ccTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCAGCACCCCGCCAGGGGCAGCGCTCAGGCAAACCAGTGCCGCCAGGAAAGCAAACAATGGCAAGCGCCTCGCTTAta                              >  7:131417/1‑134 (MQ=255)
                                                                                            cTTGGCTACGCGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATcc                             <  1:269888/135‑1 (MQ=255)
                                                                                                           ggCCGCGGCCGGCGCCGCAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc              <  3:159909/135‑1 (MQ=255)
                                                                                                           ggCAGCGGCCGGCGCGGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc              <  2:13850/135‑1 (MQ=255)
                                                                                                           ggCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc              <  4:56803/135‑1 (MQ=255)
                                                                                                              aGCGGAAGGCGCCGGATACCGCGCGGGCACGTCTGCGGCCTGCACCCGGCCAGGGGCCGCCCGCAGGCAAAGCAGGGCCGCCAGGAAACCAAACAATGGCATGCGCATCGCTTATCCCCAAGAGGTCACACcgcg           >  6:44228/1‑135 (MQ=255)
                                                                                                                 ggACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg        <  8:22348/135‑1 (MQ=255)
                                                                                                                 gctcggcggcggcGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg        <  5:274657/131‑1 (MQ=255)
                                                                                                                     ggCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACa                                                               <  7:272366/76‑1 (MQ=255)
                                                                                                                     ggCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACa                                                               >  8:272366/1‑76 (MQ=255)
                                                                                                                       cgcCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAg  <  8:20096/135‑1 (MQ=255)
                                                                                                                              |                                                                                                                               
GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAG  >  ACB122/3986496‑3986749

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