Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 1,168,401 C→A S606* (TCG→TAG)  LY180_05445 → glucosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,168,4010CA95.7% 65.4 / NA 24S606* (TCG→TAG) LY180_05445glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (14/8);  minor base G (0/1);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCGGTGATTTCCAGCCG  >  CP006584/1168112‑1168692
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                   
ctgggttcGGTGTGCTTGGCGGGGGTGTTGTGTTATCTTGCCGCCCCGCTGTGGTTTATGTTCCTAGCGCTCTCTAATGAATTGGCGTGTGTACATTGGGTTACCTAACCCGAATACTTCCTGCAACCACGGCAGTTGTTCCCGGGTGGGCCGCTGTGGAGTCCTGAGCTGGCGATTGCACTTTTTGCTTCTACTATGGTGCTGTTGTTCCTGCCGAAGTGATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCGTCTGGCGCGTTACATTATAGTtgctgc                                                                                                                                                                                                                                                                                             <  2:289547/291‑1 (MQ=255)
   cGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTgg                                                                                                                                                                                                                                                                                          >  1:49027/1‑298 (MQ=255)
           gtTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCtt                                                                                                                                                                                                                                                                                  <  1:159467/298‑1 (MQ=255)
             tCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCtttt                                                                                                                                                                                                                                                                                >  1:46880/1‑298 (MQ=255)
                    gggCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGtg                                                                                                                                                                                                                                                                         >  1:301587/1‑298 (MQ=255)
                            tGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGc                                                                                                                                                                                                                                                                 >  1:259535/1‑298 (MQ=255)
                                      ctccgctccgCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGc                                                                                                                                                                                                                                                       >  1:293530/1‑298 (MQ=255)
                                                         aTGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCt                                                                                                                                                                                                                                    <  1:57229/298‑1 (MQ=255)
                                                              ccTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCCGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGtt                                                                                                                                                                                                                               <  1:292452/298‑1 (MQ=255)
                                                                  gcgcTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTGCCATACGGTCTTCGTTGTCa                                                                                                                                                                                                                           >  1:54721/1‑298 (MQ=255)
                                                                                                                                                            gtgcggctggAGCGGGCGATTCCACTTTTTGCTTCGACTATGTGCGTTTTGTTCCTGCCGCAGTTATTGAGCATTTTGCTTAGTTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCt                                                                                                                                 <  2:293530/289‑1 (MQ=255)
                                                                                                                                                                        cTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGGGGTGTGGAATTAACAGCAGCGTGATGATGAGTCCCCTTCCTGGGGGGAAGCGTTCAAACGCCCCGGcgccaaggtgttgc                                                                                                                     >  2:242869/1‑285 (MQ=255)
                                                                                                                                                                         tGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGtt                                                                                                                    >  1:299036/1‑298 (MQ=255)
                                                                                                                                                                             gATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCGGTTAggg                                                                                                                >  1:292704/1‑298 (MQ=255)
                                                                                                                                                                             gATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAAGTCACCGCAGCGTGATGATGACTCCACTTCCTGGTGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAggg                                                                                                                >  1:299709/1‑298 (MQ=255)
                                                                                                                                                                                     tttttGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGtgtg                                                                                                        <  1:242869/298‑1 (MQ=255)
                                                                                                                                                                                                  tATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATg                                                                                           >  1:315814/1‑298 (MQ=255)
                                                                                                                                                                                                                   cTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGTGTGGCTGGATCTGc                                                                           >  1:256740/1‑298 (MQ=255)
                                                                                                                                                                                                                              ttggggattttttgttATCTGGGGAAAAGGAACGAAAGAATTGGGGGTGTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTGCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTc                                                                <  2:299036/284‑1 (MQ=255)
                                                                                                                                                                                                                                           ttttgtgggggCAAAGGAACGGAGGTATAGGGCGGCTTCTGGCTCGTGTCGTTATGGGTGGTGCTGGAAGTGTTTTGTTCCGTGCTGGTGGCTCCGGTACGCATGCTGTTGCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGa                                                   <  2:181569/288‑1 (MQ=255)
                                                                                                                                                                                                                                           cTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGa                                                   >  1:181569/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                        cggcggCTTCTGGCGCGTAACATTATATTTGCTGCTGGAAGTGCTTTTTTCCTTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCgttt                      <  1:126392/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gcTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCCTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGCGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTCCTGATAGGGTTAGTGGGGGCTGTTTGGCTGGCGTGGCCGCATCTGCGTCTCCTGTTCTGGCTGGCCCCGGTTTTCTTCTCGTTGTCCCTTCCACCgtttgtttt                 >  2:129470/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                                            cATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGCAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCTGGGCTGGCTCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCCGGTTGATCCTGTCACCTTTTGTTTCGGTTGTTTCCAGCCg  >  2:286796/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                   
CACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGGTGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCGGTGATTTCCAGCCG  >  CP006584/1168112‑1168692

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: