Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,669,340 A→G intergenic (‑152/‑252) yhjA ← / → treF putative cytochrome C peroxidase/cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,669,3400AG84.8% 194.9 / 24.9 79intergenic (‑152/‑252)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (8/4);  new base G (34/33);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCG  >  NC_000913/3669244‑3669435
                                                                                                |                                                                                               
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                              >  2:1779481/1‑100 (MQ=255)
gcAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAg                                                                                              >  1:1383131/1‑100 (MQ=255)
     cGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgt                                                                                         <  2:1680836/100‑1 (MQ=255)
      gTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCtgtg                                                                                        <  2:1813765/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                      <  2:1570470/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                      <  1:1211069/100‑1 (MQ=255)
        ctctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGcc                                                                                      >  2:85635/1‑100 (MQ=255)
          ctGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTa                                                                                    <  2:21363/100‑1 (MQ=255)
          ctGCGAGACAAAGACATCATAAAAAACAGGGAGAAACGATTTTTGATAGGAATCAATAAATTACATGACTATAAACTATCACTAAGGTAGCCTGTGCCt                                                                                     >  1:1203117/1‑99 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                  >  2:941120/1‑100 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTAtt                                                                                  >  2:1873685/1‑100 (MQ=255)
            gCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTAtt                                                                                  >  1:1829138/1‑100 (MQ=255)
                gaCAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                              <  2:900723/100‑1 (MQ=255)
                 aCAAAGACAACATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTAccc                                                                              <  1:1274424/99‑1 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTTTTTTTTTTTAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                            >  1:1878372/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                            >  2:754497/1‑100 (MQ=255)
                  cAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAg                                                                            >  2:1003903/1‑100 (MQ=255)
                   aaaGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGa                                                                           >  1:566505/1‑100 (MQ=255)
                      gCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCa                                                                        >  1:478514/1‑100 (MQ=255)
                       ccATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAg                                                                       >  2:1825332/1‑100 (MQ=255)
                         atcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATGAAATATCAAGAAGGTAGCCTGTGCCTATTACCCAGATCAGGa                                                                     >  1:1364084/1‑100 (MQ=255)
                         atcatAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGa                                                                     >  2:1275472/1‑100 (MQ=255)
                          tcatAAAAAACAGGGAGATACGCTTTTTGATAGGATTCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGaa                                                                    >  1:1104631/1‑100 (MQ=255)
                           catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                   >  2:215804/1‑100 (MQ=255)
                           catAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAAt                                                                   >  2:33489/1‑100 (MQ=255)
                            atgaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                  <  2:1450039/97‑1 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                  <  2:1470459/100‑1 (MQ=255)
                            atAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                  <  1:577548/100‑1 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTTTTTTTTTTTAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATg                                                                  >  1:202718/1‑99 (MQ=255)
                             tAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATgg                                                                 >  1:1863775/1‑100 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                <  1:366881/100‑1 (MQ=255)
                              aaaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGc                                                                >  1:1716419/1‑100 (MQ=255)
                               aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                               <  1:538740/100‑1 (MQ=255)
                               aaaaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCa                                                               <  1:777148/100‑1 (MQ=255)
                                 aaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCata                                                             <  2:25101/100‑1 (MQ=255)
                                 aaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCata                                                             <  2:25046/100‑1 (MQ=255)
                                 aaaCAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCata                                                             <  1:1688743/100‑1 (MQ=255)
                                     aGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCAt                                                         <  2:672096/100‑1 (MQ=255)
                                      gggAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATg                                                        >  1:1126188/1‑100 (MQ=255)
                                        gagaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGAt                                                      >  1:1830958/1‑100 (MQ=255)
                                          gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTg                                                    >  2:673408/1‑100 (MQ=255)
                                          gaTACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTg                                                    >  2:1890825/1‑100 (MQ=255)
                                             aCGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTgcgg                                                 <  2:741301/100‑1 (MQ=255)
                                                cTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCga                                              >  2:453646/1‑100 (MQ=255)
                                                 tttttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCgag                                             <  1:246302/100‑1 (MQ=255)
                                                   tttGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGagag                                           <  2:521964/100‑1 (MQ=255)
                                                      gATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAg                                        >  1:834841/1‑100 (MQ=255)
                                                       aaagcaaccaaacaacTACAGTAACAGAAACTATCAAGAAAACATCCTTTGCCTACTACCCAGAGCAGTAATATCATATCATGATTGCGGCgaga                                            <  2:848048/87‑1 (MQ=255)
                                                       aTAGGAATCAATCAATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAgg                                       <  2:1814956/100‑1 (MQ=255)
                                                             atcaatcaATTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGaa                                 >  1:1417205/1‑100 (MQ=255)
                                                              tcaatcaaTTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAAt                                <  2:1238143/100‑1 (MQ=255)
                                                               caatcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATg                               <  2:850129/100‑1 (MQ=255)
                                                                aatcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGt                              >  1:1869077/1‑100 (MQ=255)
                                                                  tcaatTACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTc                            <  1:1890812/100‑1 (MQ=255)
                                                                      ttACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaa                        <  2:1830970/100‑1 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                       <  2:1110509/100‑1 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATGACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                       <  1:193790/100‑1 (MQ=255)
                                                                       tACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaa                       >  1:465898/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                     >  1:1024411/1‑100 (MQ=255)
                                                                         cATGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaat                     <  2:1500212/100‑1 (MQ=255)
                                                                          aTGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatg                    <  2:1383139/100‑1 (MQ=255)
                                                                           tGAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                   >  1:1711182/1‑100 (MQ=255)
                                                                           tGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatga                   <  2:624184/100‑1 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaa                  <  1:146077/100‑1 (MQ=255)
                                                                            gAATATAAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaa                  <  2:878002/100‑1 (MQ=255)
                                                                             aaTATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCtgaaaatgaaa                 >  2:980650/1‑100 (MQ=255)
                                                                               tataAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGaaaaa               >  2:1880091/1‑100 (MQ=255)
                                                                               tataAACTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGaaaaa               >  2:687995/1‑100 (MQ=255)
                                                                                 taAACTATCAATAAGGTAGCCTGTGCCTATTACCAAGATAAGGAATGGCATATCATGCTTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAaaca             >  1:1249791/1‑100 (MQ=255)
                                                                                    aCTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAacaaca          >  2:76623/1‑100 (MQ=255)
                                                                                     cTATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAacac         <  2:443252/100‑1 (MQ=255)
                                                                                      tATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCTGCGAGAGCAGAATGGAATGTTCTGAAAATGAAAAACAACAcc        >  2:1801591/1‑100 (MQ=255)
                                                                                      tATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACAcc        <  2:1878384/100‑1 (MQ=255)
                                                                                      tATCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACAcc        >  1:33611/1‑100 (MQ=255)
                                                                                       aTCAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCa       <  1:1483146/100‑1 (MQ=255)
                                                                                         cAATAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCaaa     >  1:515571/1‑100 (MQ=255)
                                                                                          aaTAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAAt    >  2:15668/1‑100 (MQ=255)
                                                                                            tAAGGTAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCg  <  1:85635/100‑1 (MQ=255)
                                                                                            tAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGTGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCg  >  1:1816369/1‑100 (MQ=255)
                                                                                                |                                                                                               
GCAAACGTCTCTGCGAGACAAAGCCATCATAAAAAACAGGGAGATACGCTTTTTGATAGGAATCAATCAATTACATGAATATAAACTATCAATAAGATAGCCTGTGCCTATTACCCAGATCAGGAATGGCATATCATGATTGCGGCGAGAGCAGGATTGAATGTTCTGAAAATGAAAAACAACACCAAATCG  >  NC_000913/3669244‑3669435

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