Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,888,058 C→T 7.2% intergenic (‑67/+3) rnd ← / ← fadD ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,888,0580CT 7.2% 186.2 / 4.4 57intergenic (‑67/+3)rnd/fadDribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base C (32/20);  new base T (2/2);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

TCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAT  >  NC_000913/1887919‑1888178
                                                                                                                                           |                                                                                                                        
tCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAc                                                                                                                          >  2:715102/1‑140 (MQ=255)
tCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAAc                                                                                                                          >  1:284866/1‑140 (MQ=255)
         aGGCACGGACGTCTTCACACAAAGAAGCCTGCGCATCTTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCTGGCt                                                                                                                 <  2:459233/140‑1 (MQ=255)
         aGGCACGGACGGCTACACACAAACAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTAAGCGGGCGACAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCt                                                                                                                 <  2:1161942/140‑1 (MQ=255)
            cACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTa                                                                                                              <  2:159632/140‑1 (MQ=255)
                 acggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATGTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTc                                                                                                         >  2:1332839/1‑140 (MQ=255)
                 acggCTTCACACAAAGAAGCCAGCGCAGCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTGTAGTGTc                                                                                                         >  2:1332505/1‑140 (MQ=255)
                     cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACt                                                                                                     >  1:842963/1‑140 (MQ=255)
                     cTTCACACAAAGAAGCCAGCGCAACGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAATCGCCGGATAAACCGGCGTCTGACTACTGACATAACGCTCATGCTTTATTGTCCACt                                                                                                     >  1:588566/1‑140 (MQ=255)
                                    ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                          >  2:260047/1‑120 (MQ=255)
                                    ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGt                                                                                                          <  1:260047/120‑1 (MQ=255)
                                    ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                      <  1:1384899/140‑1 (MQ=255)
                                            tCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTACTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                              >  2:1385357/1‑140 (MQ=255)
                                             cGTCCGTGGTAATCATTTGTTAATTCAAATTGTTTTCTCTTTAGTGGTCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCATGCTATATTGTCCACTTAGCCGCGCGCTTCGTCACGTAAt                                                                             >  2:903735/1‑140 (MQ=255)
                                               tCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtc                                                                           >  2:801549/1‑140 (MQ=255)
                                               tCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtc                                                                           >  2:287439/1‑140 (MQ=255)
                                                        aTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCaa                                                                  >  2:1664693/1‑140 (MQ=255)
                                                         tCATTTGGTAATTCAAATTGTTTTCTGTTTAGTGGGGGTGAAAATATACGCCGGATTAACCGGCGTCTGACGTCTGAGTTAACGCTCAGGGTTTATTGTCCACTTTGGCGCGCGCTTCGTCTCGTAATTCTCGTCGCaaa                                                                 >  2:144675/1‑140 (MQ=255)
                                                             ttGGTAATTGATTTTGTTTTCTCTTTAGTGGGCGTCTTATTTTTCGGCGGATTATCGGGCGTCTGACGACTGTCTTATCGCTCAGGGTTTATTGTCGACTTTGGCGCGCGGTTCGTCTCGTAATTGTCGTCGCTATAttt                                                             >  2:157510/1‑140 (MQ=255)
                                                             ttGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttt                                                             >  1:1542332/1‑140 (MQ=255)
                                                             ttGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttt                                                             >  1:7271/1‑140 (MQ=255)
                                                             ttGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttt                                                             >  1:873204/1‑140 (MQ=255)
                                                                    ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAc                                                      >  2:1376925/1‑140 (MQ=255)
                                                                    ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAc                                                      >  1:1149105/1‑140 (MQ=255)
                                                                    ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAc                                                      >  1:1325672/1‑140 (MQ=255)
                                                                    ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAc                                                      >  1:459703/1‑140 (MQ=255)
                                                                    ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCGTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAc                                                      >  1:1496626/1‑140 (MQ=255)
                                                                    ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACAGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATACTCGTCGCAAAATTTTTCCGAc                                                      >  1:968848/1‑140 (MQ=255)
                                                                    ttCAAATAGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGACTGACGACTGACTTAAAGCAAAGGCTTTATTGTCCACTTTGCCGCGCTCTTCGTCACGAAATTCTCGTCGCAAAATTTTTCCGAc                                                      >  2:533527/1‑140 (MQ=255)
                                                                     tCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                     >  1:1202653/1‑140 (MQ=255)
                                                                           tGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAt                                               <  1:277477/140‑1 (MQ=255)
                                                                                 ctctTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                         <  2:663651/140‑1 (MQ=255)
                                                                                       aGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtca                                   >  2:373749/1‑140 (MQ=255)
                                                                                       aGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtca                                   >  2:373650/1‑140 (MQ=255)
                                                                                       aGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtca                                   >  2:373544/1‑140 (MQ=255)
                                                                                             cGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                             >  2:1117538/1‑140 (MQ=255)
                                                                                             cGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                             >  2:1182905/1‑140 (MQ=255)
                                                                                             cGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                             <  2:238919/140‑1 (MQ=255)
                                                                                              gTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTa                                                                                                              >  2:43908/1‑58 (MQ=255)
                                                                                              gTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTa                                                                                                              <  1:43908/58‑1 (MQ=255)
                                                                                              gTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                            >  1:781042/1‑140 (MQ=255)
                                                                                              gTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                            <  1:290911/140‑1 (MQ=255)
                                                                                              gTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                            >  1:780836/1‑140 (MQ=255)
                                                                                                      aaaCGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                    <  1:1546215/140‑1 (MQ=255)
                                                                                                         cGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGc                 <  2:193605/140‑1 (MQ=255)
                                                                                                         cGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATATTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGc                 <  2:796078/140‑1 (MQ=255)
                                                                                                             ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg             <  1:535174/140‑1 (MQ=255)
                                                                                                             ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg             <  1:533873/140‑1 (MQ=255)
                                                                                                             ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg             <  1:194395/140‑1 (MQ=255)
                                                                                                             ggATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg             <  1:715102/140‑1 (MQ=255)
                                                                                                               aTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt           <  2:1290934/140‑1 (MQ=255)
                                                                                                               aTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt           <  1:345381/140‑1 (MQ=255)
                                                                                                               aTTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt           <  1:345320/140‑1 (MQ=255)
                                                                                                                        cGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAt  >  2:1273651/1‑140 (MQ=255)
                                                                                                                        cGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAt  >  2:1273645/1‑140 (MQ=255)
                                                                                                                                        tAATGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt               >  1:556439/1‑111 (MQ=255)
                                                                                                                                        tAATGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt               <  2:556378/111‑1 (MQ=255)
                                                                                                                                        tAATGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt               <  2:556439/111‑1 (MQ=255)
                                                                                                                                        tAATGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt               >  1:556378/1‑111 (MQ=255)
                                                                                                                                           |                                                                                                                        
TCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAT  >  NC_000913/1887919‑1888178

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