Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,015,583 A→C 8.3% M142L (ATG→CTG) ‡ mocA → CTP:molybdopterin cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,015,5830AC 8.3% 159.1 / 5.2 48M142L (ATG→CTG) ‡mocACTP:molybdopterin cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base C (2/2);  total (22/26)
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.96e-01

AACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATA  >  NC_000913/3015437‑3015727
                                                                                                                                                  |                                                                                                                                                
aaCATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATg                                                                                                                                                <  4:286557/149‑1 (MQ=255)
aaCATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCAATTTTAGTGTCAAAACCATGCCTGATg                                                                                                                                                >  4:201116/1‑149 (MQ=255)
    ttGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAgg                                                                                                                                            >  2:260262/1‑149 (MQ=255)
       ttttCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCaa                                                                                                                                         <  3:11577/149‑1 (MQ=255)
            tcACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTTGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCatggttt                                                                                                                                    >  1:307085/1‑143 (MQ=255)
            tcACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCatggttt                                                                                                                                    <  2:307085/149‑7 (MQ=255)
                ccACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGAc                                                                                                                                >  4:141110/1‑149 (MQ=255)
                ccACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGAc                                                                                                                                <  2:68384/149‑1 (MQ=255)
                          aTGCCGACCCTCACCATCGATATTTTTAGAAAAATCTAGTCGTTACAAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCAc                                                                                                                      <  2:88071/149‑1 (MQ=255)
                          aTGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCAc                                                                                                                      <  1:149593/149‑1 (MQ=255)
                             ccAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCaa                                                                                                                   >  4:84324/1‑149 (MQ=255)
                                aCCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAatat                                                                                                                <  3:416088/149‑1 (MQ=255)
                                  ccTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGc                                                                                                              <  3:356066/149‑1 (MQ=255)
                                  ccTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGc                                                                                                              >  2:76503/1‑149 (MQ=255)
                                  ccTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGc                                                                                                              >  2:50502/1‑149 (MQ=255)
                                    tCACCATCGATATTTTTAGAAAGATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACGAATATGCGt                                                                                                            <  1:408882/149‑1 (MQ=255)
                                      aCCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCaa                                                                                                                                         >  1:17720/1‑118 (MQ=255)
                                      aCCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCaa                                                                                                                                         <  2:17720/118‑1 (MQ=255)
                                                      gAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCat                                                                                                                                         >  4:291641/1‑101 (MQ=255)
                                                      gAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGCATCAGGCat                                                                                                                                         <  3:291641/102‑2 (MQ=255)
                                                      gAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCAATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGga                                                                                          <  3:201116/149‑1 (MQ=255)
                                                       aaaaaTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGgag                                                                                         >  2:300702/1‑149 (MQ=255)
                                                       aaaaaTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGgag                                                                                         >  2:220479/1‑149 (MQ=255)
                                                         aaaTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGAt                                                                                       >  4:243753/1‑149 (MQ=255)
                                                           aTCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCa                                                                                     >  2:69305/1‑149 (MQ=255)
                                                                   gTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGACTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCg                                                                             >  1:271706/1‑149 (MQ=255)
                                                                    ttACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGt                                                                            >  2:421436/1‑149 (MQ=255)
                                                                        gAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGaa                                                                        >  3:440144/1‑149 (MQ=255)
                                                                                 gcgcAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAAt                                                               <  3:24116/149‑1 (MQ=255)
                                                                                    cAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCa                                                            <  4:280280/149‑1 (MQ=255)
                                                                                            ccACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAaataat                                                    >  1:212938/1‑149 (MQ=255)
                                                                                             cACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAAtt                                                   <  1:69305/149‑1 (MQ=255)
                                                                                               cTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAAtttt                                                 <  1:300702/149‑1 (MQ=255)
                                                                                                tCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTa                                                >  2:221234/1‑149 (MQ=255)
                                                                                                  cATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGa                                              >  3:157281/1‑149 (MQ=255)
                                                                                                  cATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGa                                              <  1:260262/149‑1 (MQ=255)
                                                                                                       tGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGGCCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTg                                         <  2:446521/149‑1 (MQ=255)
                                                                                                        ggCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGa                                        <  3:243753/149‑1 (MQ=255)
                                                                                                            tCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACt                                    >  3:177830/1‑149 (MQ=255)
                                                                                                                          aTTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttat                      <  1:421436/149‑1 (MQ=255)
                                                                                                                            tttAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTttatta                    <  1:50502/149‑1 (MQ=255)
                                                                                                                                   tCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGaa             <  4:177830/149‑1 (MQ=255)
                                                                                                                                     aaaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAgg           >  3:2108/1‑149 (MQ=255)
                                                                                                                                      aaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTAt                                                                                              <  2:358756/65‑1 (MQ=255)
                                                                                                                                      aaaCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTAt                                                                                              >  1:358756/1‑65 (MQ=255)
                                                                                                                                        aCCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGaaa        <  4:2108/149‑1 (MQ=255)
                                                                                                                                            tGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGta    <  4:157281/149‑1 (MQ=255)
                                                                                                                                              ccTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGtata  <  2:374529/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AACATTGTTTTCTCACCCACGGTGACATGCCAACCCTCACCATCGATATTTTTAGAAAAATCTGGTCGTTACGAAATGATGGCGCAATACTGCCACTCCATAATGGCATCCCCGGCCATCCGATTTTAGTGTCAAAACCATGCCTGATGCAGGCAATCCAACGACCCAATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATA  >  NC_000913/3015437‑3015727

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