Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,107,469 +C intergenic (‑296/‑101) b2965 ← / → yqgA ornithine decarboxylase isozyme/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,107,4691.C97.3% 386.3 / 5.0 147intergenic (‑296/‑101)b2965/yqgAornithine decarboxylase isozyme/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (58/85);  total (61/86)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 155 total aligned reads displayed.

GAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGA  >  NC_000913/3107394‑3107538
                                                                             |                                                                     
  gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                      <  1:4743166/76‑3 (MQ=255)
  gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                      >  2:5247773/1‑74 (MQ=255)
  gaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                      >  1:3378694/1‑74 (MQ=255)
  gaATGTAAGGAGAGGGCGTGCCTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcg                                                                      >  1:2979834/1‑74 (MQ=255)
   aaTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGCGGAGAAAATTCGCGCAATCATGGCACCcttcgc                                                                     >  2:890803/1‑73 (MQ=255)
   aaTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                     >  2:1621133/1‑73 (MQ=255)
       tAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTaaa                                                                 <  1:3873135/76‑1 (MQ=255)
        aaGGAGATGGCGTGCATCGGGCAAACTCCGTACATAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                <  2:112391/76‑1 (MQ=255)
        aaGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAg                                                                >  1:4493589/1‑76 (MQ=255)
         aGGAGATGGCGTGCATCGGGCGAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                               <  2:4857552/76‑1 (MQ=255)
           gagaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGt                                                             <  2:5131462/76‑1 (MQ=255)
           gagaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGt                                                             <  1:5754290/76‑1 (MQ=255)
            agaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTg                                                            >  2:3240784/1‑76 (MQ=255)
             gaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGt                                                             <  2:4063102/74‑1 (MQ=255)
             gaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGt                                                             >  1:4063103/1‑74 (MQ=255)
             gaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGc                                                           >  2:635450/1‑76 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         <  1:3346547/76‑1 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         <  2:1678277/76‑1 (MQ=255)
                ggCGTGCATCGGGCAACCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAg                                                        <  1:4630343/76‑1 (MQ=255)
                ggCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAg                                                        <  1:825387/76‑1 (MQ=255)
                ggCGGGCATAGGACAAATTCCCGACGGGGAGGAGAAAGTGAGCGCTGTCATGGCACCTTTCGCTAAAGCGTGcgag                                                        >  1:770423/1‑73 (MQ=255)
                 gCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGc                                                       <  1:77019/76‑1 (MQ=255)
                  cGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGcc                                                      >  1:227980/1‑76 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                     >  1:1075402/1‑76 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                     >  1:4128036/1‑76 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTGTCGCTAAAGCGTGCAAGCCa                                                     >  2:3805963/1‑76 (MQ=255)
                     gaaTCGGGCAAACTCCGAACAGAGAGGAGAAAATTAGAGCAATCATGGAACCATTCGCTAAAGCGTGCAAGCCAgg                                                   >  1:5235057/3‑76 (MQ=255)
                     gCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                   >  1:2558897/1‑76 (MQ=255)
                     gCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                   <  1:3524937/76‑1 (MQ=255)
                     gCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                   >  1:1865460/1‑76 (MQ=255)
                       aTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                 >  1:5011154/1‑76 (MQ=255)
                         cGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAGGCGTGCAAGCCAGGACtt                                               <  2:518666/76‑1 (MQ=255)
                         cGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAAtt                                               >  1:5620211/1‑76 (MQ=255)
                          gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                              <  1:3618332/76‑1 (MQ=255)
                          gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                              <  2:1222898/76‑1 (MQ=255)
                          gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                              >  2:3878164/1‑76 (MQ=255)
                           ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCACGCCAGGACTTTg                                             >  2:1752074/1‑76 (MQ=255)
                           ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTg                                             <  2:1646012/76‑1 (MQ=255)
                            gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGc                                            >  1:5346383/1‑76 (MQ=255)
                            gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGc                                            >  1:3271586/1‑76 (MQ=255)
                             cAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCg                                           >  2:1891781/1‑76 (MQ=255)
                              aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGa                                          <  2:358138/76‑1 (MQ=255)
                              aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGa                                          <  2:3191627/76‑1 (MQ=255)
                              aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGa                                          >  1:3436186/1‑76 (MQ=255)
                               aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaa                                         >  1:21441/1‑76 (MQ=255)
                               aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaa                                         <  1:271958/76‑1 (MQ=255)
                                 cTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaaca                                       <  1:2998149/76‑1 (MQ=255)
                                 cTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaaca                                       >  1:3180115/1‑76 (MQ=255)
                                  tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaa                                      <  1:1628431/76‑1 (MQ=255)
                                  tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaa                                      <  2:4128035/76‑1 (MQ=255)
                                  tCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaa                                      >  2:2100928/1‑76 (MQ=255)
                                  tCCGTACAGAAAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaa                                      <  1:573661/76‑1 (MQ=255)
                                   ccGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaa                                     <  2:4194814/76‑1 (MQ=255)
                                   ccGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaa                                     <  2:2299635/76‑1 (MQ=255)
                                    cGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacacac                                    >  2:1034052/1‑76 (MQ=255)
                                    cGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaac                                    >  1:1483481/1‑76 (MQ=255)
                                    cGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaac                                    <  1:5401460/76‑1 (MQ=255)
                                     gTACAGAGAGGAGAAAATTCGCGCAATCGTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                   <  1:3772332/76‑1 (MQ=255)
                                     gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                   <  2:3875942/76‑1 (MQ=255)
                                     gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                   <  2:2541135/76‑1 (MQ=255)
                                      tACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAacaaacaa                                  <  2:4778852/76‑1 (MQ=255)
                                      tACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAacaaacaa                                  <  1:4035892/76‑1 (MQ=255)
                                       aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGCACAAACAAg                                 >  2:1896636/1‑76 (MQ=255)
                                       aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                 >  1:2515455/1‑76 (MQ=255)
                                        cAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGc                                >  1:2308904/1‑76 (MQ=255)
                                         agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                               <  2:940142/76‑1 (MQ=255)
                                         agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                               <  2:696327/76‑1 (MQ=255)
                                         agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                               >  1:2466121/1‑76 (MQ=255)
                                           agagGAGAAAATTCGCGTAATCATGGCCCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                             <  2:5671227/76‑1 (MQ=255)
                                              ggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGa                          <  2:1034257/76‑1 (MQ=255)
                                               gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                         <  1:5566194/76‑1 (MQ=255)
                                                agaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                        <  2:2935526/76‑1 (MQ=255)
                                                agaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                        >  1:442066/1‑76 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      <  2:3217119/76‑1 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      <  1:4122168/76‑1 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      <  2:2064026/76‑1 (MQ=255)
                                                     aTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAcc                   >  1:4552545/1‑76 (MQ=255)
                                                     aTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAcc                   <  2:4478204/76‑1 (MQ=255)
                                                     aTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAcc                   <  1:788593/76‑1 (MQ=255)
                                                      ttCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCa                  <  1:3228010/76‑1 (MQ=255)
                                                          cgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACg              <  1:170687/76‑1 (MQ=255)
                                                           gcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             >  2:4490983/1‑76 (MQ=255)
                                                           gcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             <  1:5095094/76‑1 (MQ=255)
                                                            cAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGt            <  1:723850/76‑1 (MQ=255)
                                                             aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTc           <  1:2139847/76‑1 (MQ=255)
                                                             aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTc           >  2:670790/1‑76 (MQ=255)
                                                             aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTc           <  1:1785901/76‑1 (MQ=255)
                                                              aTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCa          >  1:3619277/1‑76 (MQ=255)
                                                               tgaTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat         <  2:1804686/74‑1 (MQ=255)
                                                               tCATGGCGCCTTTCGCTAAAGCGTGCACGCCAGGCCTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat         <  2:81423/76‑1 (MQ=255)
                                                                cATGGCTCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata        <  1:3691310/76‑1 (MQ=255)
                                                                cATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCata        >  2:3589741/1‑76 (MQ=255)
                                                                 aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataa       >  1:1443629/1‑76 (MQ=255)
                                                                  tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGGCTTTGCGAACAAACAGGCCATCGACGGCACCAAACGTTCgtaat      <  2:395494/76‑1 (MQ=255)
                                                                  gggCGCCTTTCGCTAAAGCGTGCAAGCCGGGACTTTGCGAACAACCAAGCCAGCGACGGCACCAAACGGTCataat      <  2:2979834/75‑1 (MQ=255)
                                                                   ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata     <  1:2703466/76‑1 (MQ=255)
                                                                    gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa    <  2:3116888/76‑1 (MQ=255)
                                                                    gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa    >  2:2371219/1‑76 (MQ=255)
                                                                     cACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAg   >  1:117734/1‑76 (MQ=255)
                                                                      aCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGa  >  1:4662876/1‑76 (MQ=255)
                                                                             |                                                                     
GAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGA  >  NC_000913/3107394‑3107538

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