Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,888,023 (A)8→9 intergenic (‑32/+38) rnd ← / ← fadD RNase D/fatty acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,888,0151.A98.1% 178.8 / ‑3.7 52intergenic (‑24/+46)rnd/fadDRNase D/fatty acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (31/20);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAC  >  NC_000913/1887877‑1888149
                                                                                                                                           |                                                                                                                                       
ggTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑a                                                                                                                                       >  2:1085522/1‑140 (MQ=255)
ggTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑a                                                                                                                                       >  1:1229135/1‑140 (MQ=255)
          gcgcGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑aaaaaaaaacg                                                                                                                             <  2:1229135/140‑4 (MQ=255)
          gcgcGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑aaaaaaaaacg                                                                                                                             <  1:485179/140‑4 (MQ=255)
           cgcgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑aaaaaaaaacgc                                                                                                                            >  2:4953/1‑136 (MQ=255)
           cgcgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑aaaaaaaaacgc                                                                                                                            >  2:1118495/1‑136 (MQ=255)
                    aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAAc                                                                                                                    <  2:462786/140‑1 (MQ=255)
                                tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTgacg                                                                                                        >  2:1154004/1‑140 (MQ=255)
                                tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑AAAAAAAACGCCGGATTAACCGGCGTCTgacga                                                                                                       <  1:1121118/140‑1 (MQ=255)
                                 ccAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTgacga                                                                                                       <  1:4953/140‑1 (MQ=255)
                                   aGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACgact                                                                                                     >  1:387384/1‑140 (MQ=255)
                                   aGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACgact                                                                                                     >  2:118554/1‑140 (MQ=255)
                                    gggCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACgactg                                                                                                    <  1:1118495/140‑1 (MQ=255)
                                     ggCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACgactga                                                                                                   >  1:1203741/1‑140 (MQ=255)
                                         aTCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTa                                                                                               <  1:1154004/140‑1 (MQ=255)
                                          tCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAGTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACtt                                                                                                >  2:1394326/1‑138 (MQ=255)
                                          tCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACtt                                                                                                <  1:1394326/138‑1 (MQ=255)
                                          tCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTaa                                                                                              >  1:863797/1‑140 (MQ=255)
                                               ggAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTc                                                                                         >  2:1079062/1‑140 (MQ=255)
                                                gAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAg                                                                                       >  1:776482/1‑140 (MQ=255)
                                                      cACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCttt                                                                                  >  1:1135442/1‑140 (MQ=255)
                                                               cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCAc                                                                         >  2:846422/1‑140 (MQ=255)
                                                               cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCAc                                                                         >  1:876884/1‑140 (MQ=255)
                                                                     acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGc                                                                    >  2:802880/1‑140 (MQ=255)
                                                                         agaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgc                                                               >  1:187087/1‑140 (MQ=255)
                                                                          gaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                              <  1:1079062/140‑1 (MQ=255)
                                                                            agCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                            >  2:105472/1‑140 (MQ=255)
                                                                            agCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                            <  1:256036/140‑1 (MQ=255)
                                                                             gCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCtt                                                           >  1:942956/1‑140 (MQ=255)
                                                                               cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                         <  2:187087/140‑1 (MQ=255)
                                                                               cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                         <  2:863797/140‑1 (MQ=255)
                                                                                 gcgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGTTTTATTGTCCTCTTTGCcgcgc                                                               <  2:762711/132‑1 (MQ=255)
                                                                                 gcgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                              >  2:1149122/1‑133 (MQ=255)
                                                                                 gcgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                              <  1:1149122/133‑1 (MQ=255)
                                                                                 gcgcATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgc                                                               >  1:762711/1‑132 (MQ=255)
                                                                                       cGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                 >  1:1314171/1‑140 (MQ=255)
                                                                                            gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATtctc                                            <  1:846422/140‑1 (MQ=255)
                                                                                                tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgtcg                                        <  1:118554/140‑1 (MQ=255)
                                                                                                    cATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCaaa                                    >  2:1524123/1‑140 (MQ=255)
                                                                                                       ttGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAtt                                 >  2:425459/1‑140 (MQ=255)
                                                                                                          gTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                              >  1:365076/1‑140 (MQ=255)
                                                                                                           tAATTCAAATAGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTc                             <  2:942956/140‑1 (MQ=255)
                                                                                                            aaTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                            <  1:425459/140‑1 (MQ=255)
                                                                                                            aaTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCAGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTcc                            >  2:1364326/1‑140 (MQ=255)
                                                                                                             aTTCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCg                           >  1:731747/1‑140 (MQ=255)
                                                                                                              ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGa                          >  2:1157789/1‑140 (MQ=255)
                                                                                                              ttCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGa                          >  2:1071781/1‑140 (MQ=255)
                                                                                                               tCAAATTGTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGAc                         >  1:544469/1‑140 (MQ=255)
                                                                                                                      gTTTTCTCTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAt                  <  1:105472/140‑1 (MQ=255)
                                                                                                                          tctcTTTAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCg              <  2:776482/140‑1 (MQ=255)
                                                                                                                              tttAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTaa          <  2:731747/140‑1 (MQ=255)
                                                                                                                               ttAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAAc         >  2:1281594/1‑140 (MQ=255)
                                                                                                                                tAGTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt        <  1:1085522/140‑1 (MQ=255)
                                                                                                                                  gTGGGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtca      <  2:876884/140‑1 (MQ=255)
                                                                                                                                      gCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAt                  >  1:1394/1‑124 (MQ=255)
                                                                                                                                      gCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAt                  <  2:1394/124‑1 (MQ=255)
                                                                                                                                      gCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc  >  1:962414/1‑140 (MQ=255)
                                                                                                                                      gCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc  >  2:1085487/1‑140 (MQ=255)
                                                                                                                                      gCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAc  >  2:899423/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCAC  >  NC_000913/1887877‑1888149

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