Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,862,504 Δ1 bp intergenic (‑75/‑267) msrB ← / → gapA methionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,862,5010A.98.2% 248.4 / ‑1.4 54intergenic (‑72/‑270)msrB/gapAmethionine sulfoxide reductase B/glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (25/28);  total (25/29)
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.89e-01

GATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGC  >  NC_000913/1862357‑1862637
                                                                                                                                                |                                                                                                                                        
gATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAAtttt                                                                                                                                    <  2:353342/151‑1 (MQ=255)
      gCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCt                                                                                                                              <  2:540835/150‑1 (MQ=255)
       cGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTa                                                                                                                             >  1:295736/1‑150 (MQ=255)
          cACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAAc                                                                                                                          <  1:923036/150‑1 (MQ=255)
          cACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAAc                                                                                                                          <  1:303745/150‑1 (MQ=255)
            cGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTg                                                                                                                       <  1:147366/151‑1 (MQ=255)
               aaaaCTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCTCTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACCCCAACTGGCAAATTTTGTCCTAAACTTGAtc                                                                                                                    <  1:515557/151‑1 (MQ=255)
                 aaCTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGAtctc                                                                                                                  <  2:657538/151‑1 (MQ=255)
                      cATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTcgacg                                                                                                              >  1:673224/1‑150 (MQ=255)
                              aCAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTg                                                                                                     >  1:228384/1‑151 (MQ=255)
                                 aaTTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCAc                                                                                                  >  1:103433/1‑151 (MQ=255)
                                  aTTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCAcc                                                                                                 <  1:971352/151‑1 (MQ=255)
                                        tttCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAAt                                                                                            <  1:677176/150‑1 (MQ=255)
                                         ttCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATc                                                                                           <  1:526510/150‑1 (MQ=255)
                                           cAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCg                                                                                          <  1:903258/149‑1 (MQ=255)
                                              ttcttcTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGtgatg                                                                                     >  2:735043/1‑151 (MQ=255)
                                               tcttctGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGtgatg                                                                                     <  1:314428/150‑1 (MQ=255)
                                                  tctGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGaaa                                                                                  >  1:573876/1‑150 (MQ=255)
                                                    tGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATc                                                                               >  2:279185/1‑151 (MQ=255)
                                                     gCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATca                                                                              <  1:394849/151‑1 (MQ=255)
                                                      ccGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATcac                                                                             >  1:570819/1‑151 (MQ=255)
                                                            gTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACAttttt                                                                       <  2:695963/151‑1 (MQ=255)
                                                               tattaGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTAt                                                                     <  2:119710/150‑1 (MQ=255)
                                                               tattaGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTAt                                                                     <  1:122528/150‑1 (MQ=255)
                                                                 ttaGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGt                                                                  >  2:294675/1‑151 (MQ=255)
                                                                    gCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAAt                                                               >  2:949296/1‑151 (MQ=255)
                                                                    gCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAAt                                                               <  1:950389/151‑1 (MQ=255)
                                                                       aTTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGc                                                            >  1:804751/1‑151 (MQ=255)
                                                                         ttGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGccc                                                          >  2:789412/1‑151 (MQ=255)
                                                                             tCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCttt                                                       >  1:188010/1‑150 (MQ=255)
                                                                              cacaTCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCttt                                                       <  1:337977/149‑1 (MQ=255)
                                                                              cacaTCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTaa                                                     <  2:155922/151‑1 (MQ=255)
                                                                               acaTCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTaaa                                                    <  2:6151/151‑1 (MQ=255)
                                                                                caTCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTaaa                                                    >  1:906055/1‑150 (MQ=255)
                                                                                caTCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTaaa                                                    >  2:792396/1‑150 (MQ=255)
                                                                                 aTCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAAt                                                  >  2:467330/1‑151 (MQ=255)
                                                                                   ctcACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTc                                                >  2:360860/1‑151 (MQ=255)
                                                                                     cACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGTAAATCACATTTTTATCGTAATTGCCCTTTAAAATTc                                                <  1:444094/149‑1 (MQ=255)
                                                                                     cACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAAt                                                  <  2:295702/147‑1 (MQ=255)
                                                                                        tttAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGgc                                           >  1:577526/1‑151 (MQ=255)
                                                                                         ttAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCgggg                                            <  2:270747/149‑1 (MQ=255)
                                                                                          tAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAACCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGgcg                                          <  1:68338/150‑1 (MQ=255)
                                                                                                tGCTCACATTACGGGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTAGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATCGGGGGCCCGGcccc                                  >  2:683904/1‑151 (MQ=255)
                                                                                                   tCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCAt                                >  2:979016/1‑151 (MQ=255)
                                                                                                    cacaTTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCAtg                               >  2:832225/1‑151 (MQ=255)
                                                                                                     acaTTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCAtgt                              >  2:881341/1‑151 (MQ=255)
                                                                                                         tACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGt                           <  2:224073/150‑1 (MQ=255)
                                                                                                            gTGACTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCa                       <  1:467378/151‑1 (MQ=255)
                                                                                                                cTGATTCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGcc                   >  2:935153/1‑151 (MQ=255)
                                                                                                                    ttCTAACAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCaa                >  1:767830/1‑150 (MQ=255)
                                                                                                                         aCAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAg              <  1:310154/147‑1 (MQ=255)
                                                                                                                         aCAAAACATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAag          <  1:913697/151‑1 (MQ=255)
                                                                                                                               cATTAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTg      >  1:440032/1‑149 (MQ=255)
                                                                                                                                  tAACACCAACTGGCAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGc  <  2:577462/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
GATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAGGTTTATTAGCCATTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACTGATTCTAACAAAACATTAACACCAACTGGCAAAATTTTGTCCTAAACTTGATCTCGACGAAATGGCTGCACCTAAATCGTGATGAAAATCACATTTTTATCGTAATTGCCCTTTAAAATTCGGGGCGCCGACCCCATGTGGTCTCAAGCCCAAAGGAAGAGTGAGGC  >  NC_000913/1862357‑1862637

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