Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,856,453 C→T *118* (TAG→TAA ygbA ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,4530CT97.0% 94.0 / ‑4.4 33*118* (TAG→TAAygbAuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (15/17);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856336‑2856583
                                                                                                                     |                                                                                                                                  
gcgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                               <  1:402954/139‑1 (MQ=255)
gcgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                               <  1:906814/139‑1 (MQ=255)
               cTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATCTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                                                                <  1:358874/139‑1 (MQ=255)
                       cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTc                                                                                           >  1:348179/1‑136 (MQ=255)
                       cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTc                                                                                           <  2:348179/136‑1 (MQ=255)
                         cGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCa                                                                                      <  2:428569/139‑1 (MQ=255)
                                      aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  2:300468/1‑139 (MQ=255)
                                           tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 <  2:1132604/126‑1 (MQ=255)
                                           tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 >  1:1132604/1‑126 (MQ=255)
                                           tGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCga                                                                    <  1:300468/139‑1 (MQ=255)
                                                cGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgaga                                                               >  2:429943/1‑139 (MQ=255)
                                                  gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagatg                                                             >  1:169653/1‑139 (MQ=255)
                                                       gAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                        <  2:416328/139‑1 (MQ=255)
                                                       gAATGGATACATCGGCATTGATTTAACTGCAATTTGCCGGACAGATCTGCCTGTCCGGCATACTATTAATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAAGCTGGTCATCGATGAGATGACgaa                                                        <  1:303800/139‑3 (MQ=255)
                                                             aTAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa                                                  >  1:1116367/1‑139 (MQ=255)
                                                                                aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGTATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATTAGATGACGCACGGTTAAGATCGGATGACGCCACAg                               <  1:429943/139‑1 (MQ=255)
                                                                                aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAg                               <  2:1116367/139‑1 (MQ=255)
                                                                                                   aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGat                                                                                                           <  2:956313/44‑1 (MQ=255)
                                                                                                   aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGat                                                                                                           >  1:956313/1‑44 (MQ=255)
                                                                                                     cTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTACGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACg          >  1:217887/1‑139 (MQ=255)
                                                                                                        ccTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         >  2:376447/1‑89 (MQ=255)
                                                                                                        ccTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  1:376447/89‑1 (MQ=255)
                                                                                                        ccTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAg                               >  2:1038632/1‑115 (MQ=255)
                                                                                                        ccTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAg                               <  1:1038632/115‑1 (MQ=255)
                                                                                                           gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  1:728856/1‑70 (MQ=255)
                                                                                                           gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         <  2:728856/70‑1 (MQ=255)
                                                                                                           gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt    >  1:90839/1‑139 (MQ=255)
                                                                                                           gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt    <  2:169653/139‑1 (MQ=255)
                                                                                                             ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         >  2:1145978/1‑84 (MQ=255)
                                                                                                             ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  1:1145978/84‑1 (MQ=255)
                                                                                                             ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCt  >  2:400956/1‑139 (MQ=255)
                                                                                                              cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTATCATCgatg                                                                  <  1:625968/74‑1 (MQ=255)
                                                                                                              cGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTATCATCgatg                                                                  >  2:625968/1‑74 (MQ=255)
                                                                                                                     |                                                                                                                                  
GCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856336‑2856583

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