Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,015,808 T→G intergenic (+70/‑252) mocA → / → ygfK CTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,015,8080TG96.5% 203.3 / ‑0.5 57intergenic (+70/‑252)mocA/ygfKCTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (21/34);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

AATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCG  >  NC_000913/3015664‑3015955
                                                                                                                                                |                                                                                                                                                   
aaTGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGa                                                                                                                                               <  2:706736/151‑1 (MQ=255)
 aTGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGAt                                                                                                                                              <  2:352483/151‑1 (MQ=255)
     aGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTg                                                                                                                                          >  1:838940/1‑151 (MQ=255)
     aGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTg                                                                                                                                          <  1:212205/151‑1 (MQ=255)
         ataatTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgaga                                                                                                                                       >  1:688661/1‑150 (MQ=255)
           aatTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgagaaa                                                                                                                                     <  1:731263/150‑1 (MQ=255)
              tttAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgagaaaaga                                                                                                                                  <  1:888777/150‑1 (MQ=255)
               ttAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgagaaaagag                                                                                                                                 <  1:10381/150‑1 (MQ=255)
               ttAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgagaaaag                                                                                                                                   >  1:733541/1‑148 (MQ=255)
               ttAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgagaaaag                                                                                                                                   >  2:672708/1‑148 (MQ=255)
                        gATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGa                                                                                                                        <  2:317965/150‑1 (MQ=255)
                           aCTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATaa                                                                                                                     >  1:288569/1‑150 (MQ=255)
                                    gacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGc                                                                                                           <  1:53258/151‑1 (MQ=255)
                                      cTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCa                                                                                                          >  2:408993/1‑150 (MQ=255)
                                        ttattaCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaat                                                                                                        <  1:122304/150‑1 (MQ=255)
                                          attaCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaataat                                                                                                     <  2:901340/151‑1 (MQ=255)
                                          attaCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaataa                                                                                                      <  2:154852/150‑1 (MQ=255)
                                            taCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaataata                                                                                                    >  2:315166/1‑150 (MQ=255)
                                              cAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAAtt                                                                                                 <  2:242130/151‑1 (MQ=255)
                                              cAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAAtt                                                                                                 <  1:812373/151‑1 (MQ=255)
                                                       aaaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGt                                                                                        <  1:526711/151‑1 (MQ=255)
                                                        aaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTg                                                                                       <  1:130925/151‑1 (MQ=255)
                                                           gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAtt                                                                                    >  1:534269/1‑151 (MQ=255)
                                                            tataCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAttt                                                                                   >  1:529566/1‑151 (MQ=255)
                                                                cTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAAt                                                                               >  2:813894/1‑151 (MQ=255)
                                                                  gAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAAtt                                                                              <  1:73446/150‑1 (MQ=255)
                                                                  gAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATtc                                                                             <  1:810415/151‑1 (MQ=255)
                                                                     aTTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCa                                                                          <  2:580391/151‑1 (MQ=255)
                                                                         aGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAttt                                                                       <  1:782697/150‑1 (MQ=255)
                                                                            aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGc                                                                    >  2:399176/1‑150 (MQ=255)
                                                                            aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCa                                                                   <  1:210817/151‑1 (MQ=255)
                                                                              aCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCAcc                                                                 >  1:670811/1‑151 (MQ=255)
                                                                                  aGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCtgt                                                              <  2:278106/150‑1 (MQ=255)
                                                                                  aGTTAAACGAATGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCGCCtgt                                                              <  2:398013/150‑1 (MQ=255)
                                                                                       aaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttc                                                        >  1:225900/1‑151 (MQ=255)
                                                                                        aCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttc                                                        >  1:755437/1‑150 (MQ=255)
                                                                                           aTTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcattc                                                    >  1:536503/1‑151 (MQ=255)
                                                                                              gCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcattcata                                                 <  2:769843/151‑3 (MQ=255)
                                                                                                 ttCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaa                                               >  2:334317/1‑150 (MQ=255)
                                                                                                  tCAAATATTTAAACCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaaa                                              <  2:510754/150‑1 (MQ=255)
                                                                                                     aaTATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGcc                                          <  2:786209/151‑1 (MQ=255)
                                                                                                            aaTCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAg                                    <  1:788179/150‑1 (MQ=255)
                                                                                                            aaTCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAg                                    <  2:460934/150‑1 (MQ=255)
                                                                                                            aaTCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAg                                    <  1:226264/150‑1 (MQ=255)
                                                                                                            aaTCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                   >  2:365537/1‑151 (MQ=255)
                                                                                                             aTCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGtt                                  >  1:568019/1‑151 (MQ=255)
                                                                                                               ccTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGAATTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                 <  2:194671/150‑1 (MQ=255)
                                                                                                                   ccGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                            >  2:521390/1‑151 (MQ=255)
                                                                                                                   ccGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                            <  1:178191/151‑1 (MQ=255)
                                                                                                                   ccGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAg                             <  1:111634/150‑1 (MQ=255)
                                                                                                                       cACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgccgc                         <  1:315469/150‑1 (MQ=255)
                                                                                                                          tCATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                      >  1:593632/1‑150 (MQ=255)
                                                                                                                           cATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                      <  2:414485/149‑1 (MQ=255)
                                                                                                                           cATTATCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                      <  1:421850/149‑1 (MQ=255)
                                                                                                                               aTCATAATGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                 <  1:203074/150‑1 (MQ=255)
                                                                                                                                      tGATAAGGAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat          >  2:320400/1‑150 (MQ=255)
                                                                                                                                             gAAGCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCg  >  1:615508/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCG  >  NC_000913/3015664‑3015955

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