Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,183,426 G→A Q76* (CAG→TAG)  yegX ← putative family 25 glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,183,4260GA96.9% 117.8 / ‑5.9 32Q76* (CAG→TAG) yegXputative family 25 glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/18);  minor base C (0/1);  total (13/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 = 1.00e+00

AGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTGCCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGACGCGG  >  NC_000913/2183281‑2183570
                                                                                                                                                 |                                                                                                                                                
agGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAgc                                                                                                                                             <  1:185278/151‑1 (MQ=255)
 ggCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCTCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTCCCAgc                                                                                                                                             <  2:177688/150‑1 (MQ=255)
    cATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGa                                                                                                                                         <  2:30904/151‑1 (MQ=255)
                               gcgAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACt                                                                                                              <  2:234551/151‑1 (MQ=255)
                                   aaaaTAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCa                                                                                                          <  2:259184/151‑1 (MQ=255)
                                             tCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGa                                                                                                 >  1:653783/1‑150 (MQ=255)
                                              ccaccaGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGaa                                                                                                <  2:1902/150‑1 (MQ=255)
                                                accaGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAg                                                                                             <  2:332545/151‑1 (MQ=255)
                                                  caGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCg                                                                                           >  1:238896/1‑151 (MQ=255)
                                                         tCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGAt                                                                                     <  1:453504/150‑1 (MQ=255)
                                                              ttCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAAcc                                                                               >  2:156876/1‑151 (MQ=255)
                                                                  gTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGc                                                                           >  1:667097/1‑151 (MQ=255)
                                                                       cTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATa                                                                       >  2:244909/1‑150 (MQ=255)
                                                                        ttAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATaa                                                                      >  1:118441/1‑150 (MQ=255)
                                                                                 gCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGt                                                            <  1:567852/151‑1 (MQ=255)
                                                                                                    cATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTg                                         >  2:562361/1‑151 (MQ=255)
                                                                                                    cATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTg                                         >  2:538715/1‑151 (MQ=255)
                                                                                                     aTTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTg                                         >  2:174102/1‑150 (MQ=255)
                                                                                                       tGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGaa                                       <  1:428653/150‑1 (MQ=255)
                                                                                                       tGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAAc                                      >  1:658815/1‑151 (MQ=255)
                                                                                                            cgcATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAActgctg                                 <  1:444754/151‑1 (MQ=255)
                                                                                                              cATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTc                               <  2:326776/151‑1 (MQ=255)
                                                                                                               aTTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCg                              <  2:353729/151‑1 (MQ=255)
                                                                                                                 ttttGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGa                             <  2:55756/150‑1 (MQ=255)
                                                                                                                    tGCCACCCGCTGCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAAt                           >  2:597037/1‑149 (MQ=255)
                                                                                                                               gCCAGTCGATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCAc              <  2:10393/151‑1 (MQ=255)
                                                                                                                                      gATCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTg        >  2:172002/1‑150 (MQ=255)
                                                                                                                                       aTCCGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGAc      <  2:365979/151‑1 (MQ=255)
                                                                                                                                         ccGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGAct     <  1:224334/150‑2 (MQ=255)
                                                                                                                                         ccGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGAcg     <  1:631003/150‑1 (MQ=255)
                                                                                                                                          cGCTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGAcgcg   <  2:31043/151‑1 (MQ=255)
                                                                                                                                            cTCCTACCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGACGCgg  >  2:574604/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATCCGCTCCTGCCAGCGCGAAACATCTATTCCGTGAATGGTGTAACTGGCAGGAATGCGAATAGCGAAAGATTTAACCGGGCGATAACCGTAGAAGTGGATGTAGTCTGTAATCTGAACTGCTGTCGAATAGAAAAAATTCACCGTCTGACGCGG  >  NC_000913/2183281‑2183570

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