Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 653,393 T→G 12.9% intergenic (+17/‑5) U712_03195 → / → U712_03200 Uncharacterized protein yeaD/Uncharacterized protein yebA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881653,3930TG12.9% 96.5 / 2.9 50intergenic (+17/‑5)U712_03195/U712_03200Uncharacterized protein yeaD/Uncharacterized protein yebA
Reads supporting (aligned to +/- strand):  ref base T (17/26);  new base G (5/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGGCAAAAGAGAAAGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACAGGCTTTTTTAGTG  >  CP006881/653255‑653530
                                                                                                                                          |                                                                                                                                          
tACTGGCAAAAGAGAAAGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGt                                                                                                                                            >  7:8075/1‑139 (MQ=255)
     gCAAAAGAGAAAGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTcaca                                                                                                                                      >  3:79372/1‑139 (MQ=255)
     gCAAAAGAGAAAGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTcaca                                                                                                                                      >  1:8731/1‑139 (MQ=255)
         aaGAGAAAGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGGGTTGTTTAGAAACGGAGGTGACTGGGTCACATGcc                                                                                                                                  >  2:370762/1‑139 (MQ=255)
              aaaGATGCTGAG‑‑CTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACgat                                                                                                                           <  2:221484/139‑1 (MQ=255)
                aGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACgat                                                                                                                           <  1:207912/139‑1 (MQ=255)
                 gATGGTGATCTCTCTCAAGCTTAAAATGTCATTTTTGATCGCATGTTGTAAGTCCAAATCAGAGATGGGATCATGTGGGGCGGGCGTGTTTCGAGTGTTGTTTTGAAAGTGAGGTGTCTTGTTCACATGCCGCACgat                                                                                                                           <  2:160098/138‑1 (MQ=255)
                             tctcAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAATGAGAGGTGACTGGTTCACATGCCGCACGAGGATCACAAggg                                                                                                                >  4:300753/1‑137 (MQ=255)
                             tctcAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGaa                                                                                                              <  2:8731/139‑1 (MQ=255)
                             tctcAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGGGTGGTTTTGAAAGAGAGGTGCCTGGTTCCCATGCCGCACGATGATCACAAggg                                                                                                                >  6:265166/1‑137 (MQ=255)
                              ctcAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTATCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAg                                                                                                             <  4:175913/139‑1 (MQ=255)
                                       aaaaTGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAg                                                                                                    >  2:26/1‑139 (MQ=255)
                                       aaaaTGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGGTGTTTAGAAAGAGAGGTGACTGGTTGACATGCCGCCCGCTGCTCACAAGGGAAGCCGTTTGAg                                                                                                    >  4:355203/1‑139 (MQ=255)
                                           tGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTg                                                                                                <  1:370762/139‑1 (MQ=255)
                                           tGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTg                                                                                                <  4:348284/139‑1 (MQ=255)
                                               aTGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTg                                                                                             <  1:26/138‑1 (MQ=255)
                                                 gATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGtta                                                                                          >  3:389673/1‑139 (MQ=255)
                                                           atgatgAAACGCCAAATCCGAGTTCGGATCCTGAGGGGCGGGCGAGTTTCGCGTGGTGTTTAGAAAGAGCGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTc                                                                                <  4:870/139‑1 (MQ=255)
                                                             gatgaAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGGTTTTAAAGAGAGGGGACTGGGTGACCAGGCGCCACGAGATGACAATGGAACCCGTTTGAGGTTGGGGTTTTTTTATTTTTTg                                                                              >  4:221419/1‑139 (MQ=255)
                                                               tgaAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGc                                                                            <  4:79372/139‑1 (MQ=255)
                                                                 aaaCGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCtt                                                                          <  4:105297/139‑1 (MQ=255)
                                                                     gCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGGGGTGGTTAGAAATCCAGGTGGATGGGTCACAGGCCGCACGCTGCTAACAACGGGAGCCGGTTTAGGATGGTATGATTTTTTTTTGTGGGgttttt                                                                      >  6:52771/1‑139 (MQ=255)
                                                                      ccAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTa                                                                     <  5:265166/139‑1 (MQ=255)
                                                                       cAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTAt                                                                    <  7:226728/139‑1 (MQ=255)
                                                                       cAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTAt                                                                    <  4:389673/139‑1 (MQ=255)
                                                                               agTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTgg                                                            >  5:172675/1‑139 (MQ=255)
                                                                               agTTCGGATCATGAGGGGCGGGCGAGTTTCGCGGGTTGTTTTGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGACGCCGTTTGAGGTTGCGGTTTTTTTATTTTCTGGCTTTTTTATTGCTGTgg                                                            >  5:35058/1‑139 (MQ=255)
                                                                                 ttCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGa                                                          <  1:67086/139‑1 (MQ=255)
                                                                                 ttCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGa                                                          <  6:172675/139‑1 (MQ=255)
                                                                                 ttCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGa                                                          <  5:403555/139‑1 (MQ=255)
                                                                                 ttCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCAGGCTTTTTTATTGCTGTGGGa                                                          <  5:319928/139‑1 (MQ=255)
                                                                                 gtCGGAGCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGa                                                          <  5:52771/138‑1 (MQ=255)
                                                                                         aTGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGa                                                          >  5:161893/1‑131 (MQ=255)
                                                                                         aTGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGa                                                          <  6:161893/131‑1 (MQ=255)
                                                                                         aTGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCt                                                    >  6:160969/1‑137 (MQ=255)
                                                                                         aTGAGGGGCGGGCGAGTTTCGCGGGTTGTTTAGAAAGAGAGGGGGCTGGGTCACCTGCCGCACGCTGATCACAAGGGAAGCCGTTTGAGCTTGCTGTGATTTTATTTTCTGGCTGTTTTATTGCTGTGGGGGTGGCt                                                    >  4:169004/1‑137 (MQ=255)
                                                                                         aTGAGAGGAGGGCTATTTTCTCGTGTTGTTTATAAACAGAGGTGACTGGTTCACATCCAGCACGATCATCACAAGGGAAGCCGTTATAGCATGCTGTTATATTATTTTCTGTGTTTTTTAATGCTTTGGGAGTGTCt                                                    >  6:175223/1‑137 (MQ=255)
                                                                                                gcgggcgAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTTTTGCGGTGGGAGTGGCTTTGGgcgcg                                           >  5:10176/1‑138 (MQ=255)
                                                                                                gcgggcgAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGGTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCt                                           >  7:321759/1‑139 (MQ=255)
                                                                                                  gggcgAGTTTCGTGTGTTGTTTAGTTCGAGATGTGACTGGTTCACATGCCGCACGATGAGCTAAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGTGATTGGCTTAGGCCGCTTg                                         <  8:329054/139‑1 (MQ=255)
                                                                                                   ggcgAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTa                                                                     >  8:318435/1‑110 (MQ=255)
                                                                                                   ggcgAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTa                                                                     <  7:318435/110‑1 (MQ=255)
                                                                                                                     tttAGAAAGAGAGGTGACTGGGTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTgaaacgaa                      <  8:321759/139‑1 (MQ=255)
                                                                                                                       tAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTgaaacgaaa                     >  7:57199/1‑138 (MQ=255)
                                                                                                                          aaaGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATa                 <  8:8075/139‑1 (MQ=255)
                                                                                                                            agagagGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATa                 >  2:113430/1‑137 (MQ=255)
                                                                                                                            agagagGTGACTGGGTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACa               <  3:300753/139‑1 (MQ=255)
                                                                                                                                     aCTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACAGGCttttt       >  1:18861/1‑138 (MQ=255)
                                                                                                                                         gTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACAGGCTTTTTTAgtg  <  3:355203/139‑1 (MQ=255)
                                                                                                                                         gTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACAGGCTTTTTTAgt   <  3:169004/138‑1 (MQ=255)
                                                                                                                                         gTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACAGGCTTTTTTAgt   <  5:175223/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TACTGGCAAAAGAGAAAGATGCTGAGCTCTCTCAAGCTGAAAATGTCATGATTGAACGCATGATGAAACGCCAAATCAGAGTTCGGATCATGAGGGGCGGGCGAGTTTCGCGTGTTGTTTAGAAAGAGAGGTGACTGGTTCACATGCCGCACGATGATCACAAGGGAAGCCGTTTGAGCTTGCTGTTATTTTATTTTCTGGCTTTTTTATTGCTGTGGGAGTGGCTTAGGCCGCTTGACAGTTTTACTGAAACGAAACATACAGGCTTTTTTAGTG  >  CP006881/653255‑653530

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: