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.C98.8% 452.6 / ‑0.8 169intergenic (‑296/‑101)b2965/yqgAornithine decarboxylase isozyme/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (81/86);  total (83/86)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 177 total aligned reads displayed.

GAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGA  >  NC_000913/3107394‑3107538
                                                                             |                                                                     
gagaAAGTAAGGAGAGGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCttt                                                                        <  2:380969/76‑1 (MQ=255)
 agaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                       >  2:401770/1‑75 (MQ=255)
 agaATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttc                                                                       >  2:3189758/1‑75 (MQ=255)
   aaTGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCtttcgc                                                                     >  2:2720495/1‑73 (MQ=255)
         aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                               <  1:1422885/76‑1 (MQ=255)
         aGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGc                                                               <  2:2764215/76‑1 (MQ=255)
          ggAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCg                                                              >  1:1999752/1‑76 (MQ=255)
             gaTGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCCTGc                                                           >  2:4422811/1‑76 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         <  2:770119/76‑1 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         <  1:2223835/76‑1 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         <  1:35530/76‑1 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         >  1:3328715/1‑76 (MQ=255)
               tGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCaa                                                         <  1:2500196/76‑1 (MQ=255)
                ggCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAg                                                        >  1:3303219/1‑76 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                     >  2:1883079/1‑76 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                     <  2:1845047/76‑1 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCa                                                     >  1:2729445/1‑76 (MQ=255)
                   gTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGCGCAAGCCa                                                     >  2:4713163/1‑76 (MQ=255)
                     gCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                   <  1:968887/76‑1 (MQ=255)
                     gCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAgg                                                   >  1:3252808/1‑76 (MQ=255)
                       aTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                 >  1:1152568/1‑76 (MQ=255)
                       aTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                 >  1:3903650/1‑76 (MQ=255)
                       aTCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                 >  1:2768044/1‑76 (MQ=255)
                       aTCGGACAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGAc                                                 >  1:2965161/1‑76 (MQ=255)
                          gggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACttt                                              >  1:3159706/1‑76 (MQ=255)
                           ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCGCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTg                                             <  2:520307/76‑1 (MQ=255)
                           ggCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTg                                             >  2:4773468/1‑76 (MQ=255)
                           gccAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTg                                             >  1:3946911/3‑76 (MQ=255)
                            gCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGc                                            <  1:2518658/76‑1 (MQ=255)
                             cAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTCAAGCGTGCAAGCCAGGACTTTGct                                           >  1:3983711/1‑75 (MQ=255)
                             cAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCg                                           >  1:3972577/1‑76 (MQ=255)
                              aaaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGa                                          >  2:3500835/1‑76 (MQ=255)
                              aaaCTCAGTACAGAGAGCAGAAAATGCGCGAAATCACGGCACCGGTAGATAAAGCACGCAAGCCAGGACGTTGcta                                          >  1:1667708/1‑74 (MQ=255)
                               aaCTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaa                                         >  1:550325/1‑76 (MQ=255)
                                     gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                   <  2:487707/76‑1 (MQ=255)
                                     gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                   >  2:4794411/1‑76 (MQ=255)
                                     gTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGaacaaaca                                   >  2:3524034/1‑76 (MQ=255)
                                      tACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAacaaacaa                                  >  2:649763/1‑76 (MQ=255)
                                       aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAacaaacaac                                 >  1:365994/1‑75 (MQ=255)
                                       aCAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAg                                 <  1:4123378/76‑1 (MQ=255)
                                        cAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGc                                >  1:3530256/1‑76 (MQ=255)
                                         agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                               >  1:3369980/1‑76 (MQ=255)
                                         agagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGcc                               <  2:2124905/76‑1 (MQ=255)
                                          gagagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCa                              >  2:1488437/1‑76 (MQ=255)
                                           agagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                             >  1:3418287/1‑76 (MQ=255)
                                           agagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                             <  2:4042283/76‑1 (MQ=255)
                                           agagGAGAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                             <  2:1088939/76‑1 (MQ=255)
                                           agagGAAAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCAt                             <  1:786789/76‑1 (MQ=255)
                                            gaggagAAAATTCGCGCAATCATGGCACCGTTCGCTAAAGCTTGCAAGCTAGGACATTGCGAACAAACAAGCCATc                            >  1:1370602/1‑76 (MQ=255)
                                            gaggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                            <  2:2768044/76‑1 (MQ=255)
                                            gaggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATc                            <  2:922020/76‑1 (MQ=255)
                                             aggagAAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATca                           >  1:2582258/1‑75 (MQ=255)
                                               gagaAAATTCTCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                         >  2:1932067/1‑76 (MQ=255)
                                               gagaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGAc                         <  1:3061841/76‑1 (MQ=255)
                                               aaaaaaaTCCGCGCAATCATGGCGCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACCATCAAGCCATCGAc                         <  1:60853/73‑1 (MQ=255)
                                                agaAAATTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACg                        <  1:1501563/76‑1 (MQ=255)
                                                 gaAAATTCGCGCAATCATGGCACCTTGCGCTAAAGTGGGCAAGCCAGGACTGGGCGAACAAACAAGCCATCGACgg                       <  1:4183345/76‑1 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      <  2:1786309/76‑1 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      >  2:4158891/1‑76 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      >  2:694930/1‑76 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      <  1:2423442/76‑1 (MQ=255)
                                                  aaaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGc                      >  1:3600856/1‑76 (MQ=255)
                                                   aaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCa                     <  1:4160670/76‑1 (MQ=255)
                                                   aaaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCa                     >  1:2374405/1‑76 (MQ=255)
                                                    aaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAc                    <  1:1656799/76‑1 (MQ=255)
                                                    aaTTCGCGCAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCAc                    <  1:3696897/76‑1 (MQ=255)
                                                         gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc               <  2:4558353/76‑1 (MQ=255)
                                                         gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc               <  1:4832353/76‑1 (MQ=255)
                                                         gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc               <  2:3949284/76‑1 (MQ=255)
                                                         gcgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAAc               >  1:2475597/1‑76 (MQ=255)
                                                          cgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACg              <  1:287015/76‑1 (MQ=255)
                                                          cgcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACg              <  2:3337766/76‑1 (MQ=255)
                                                           gcAATCATGTCACCTTGCGCTAAAGCGTGCAAGCCAGGACGTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             <  1:2270814/76‑1 (MQ=255)
                                                           gcAATCATGGCGCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             <  2:2292533/76‑1 (MQ=255)
                                                           gcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             <  1:2266557/76‑1 (MQ=255)
                                                           gcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             >  2:2599880/1‑76 (MQ=255)
                                                           gcAATCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACgg             <  1:3980720/76‑1 (MQ=255)
                                                             aaTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTc           >  2:4409586/1‑76 (MQ=255)
                                                              aTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAGCAAGCCATCGACGGCACCAAACGGTCa          <  2:4490262/76‑1 (MQ=255)
                                                              aTCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCa          <  2:3252808/76‑1 (MQ=255)
                                                               tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat         <  1:1529081/76‑1 (MQ=255)
                                                               tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat         <  1:3126954/76‑1 (MQ=255)
                                                               tCATGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCat         <  1:2187494/76‑1 (MQ=255)
                                                                 aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataa       >  2:4372417/1‑76 (MQ=255)
                                                                 aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataa       <  1:1193308/76‑1 (MQ=255)
                                                                 aTGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataa       <  2:1086949/76‑1 (MQ=255)
                                                                  tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGGCGGCACCAAACGGTCataat      >  1:2854070/1‑76 (MQ=255)
                                                                  tGGCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataat      >  1:4345817/1‑76 (MQ=255)
                                                                  tGGCACCTTTCGCTAAAGCGTGCAACCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataat      >  2:2410540/1‑76 (MQ=255)
                                                                   ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata     <  2:1994597/76‑1 (MQ=255)
                                                                   ggCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata     >  2:1689574/1‑76 (MQ=255)
                                                                   ggCACCTTTCGCGAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCataata     <  1:1308442/76‑1 (MQ=255)
                                                                    tcACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGAACCAAACGGTCataatca    >  1:2656094/2‑74 (MQ=255)
                                                                    gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACCACCCATCGACGGCACCAAACGGTCAaaaaaa    >  2:143047/1‑76 (MQ=255)
                                                                    gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa    >  2:1532371/1‑76 (MQ=255)
                                                                    gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa    >  1:2625563/1‑76 (MQ=255)
                                                                    gCACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCAtaataa    >  1:1422048/1‑76 (MQ=255)
                                                                     cACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAg   >  2:1244567/1‑76 (MQ=255)
                                                                     cACCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAg   >  1:3193723/1‑76 (MQ=255)
                                                                      aCCTTTCGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCCCCAAACGGTCATAATAAGa  >  2:851147/1‑76 (MQ=255)
                                                                             |                                                                     
GAGAATGTAAGGAGATGGCGTGCATCGGGCAAACTCCGTACAGAGAGGAGAAAATTCGCGCAATCATGGCACCTTTGCTAAAGCGTGCAAGCCAGGACTTTGCGAACAAACAAGCCATCGACGGCACCAAACGGTCATAATAAGA  >  NC_000913/3107394‑3107538

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